Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,206 G→T 6.7% R713L (CGT→CTT)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,2060GT 6.7% 154.9 / 2.1 60R713L (CGT→CTT) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (33/23);  new base T (2/2);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATG  >  NC_000913/1660058‑1660348
                                                                                                                                                    |                                                                                                                                              
aTAAGCCAAATGCCAGTTTATGCCTCTTCCTTTGAAGGCTGGATATCACCTGATCGTAGAACCTTCCCACTGCAACTGTGTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGACg                                                                                                                                                <  2:807873/149‑1 (MQ=255)
    gCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGATCCTTCCCACTGCAACTTGTTGGTTTCCATTACAAATCCCGTACTCACTCGACCAACGGCAATATTGATCTCCAGAAGGCTGCCTGCCGTCag                                                                                                                                            >  1:515326/1‑149 (MQ=255)
    ccccATAGCCAGTTTAAGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTACCATTACAACACCCGACCTCACTCGACCTACGGCAATATTGATCTCCTGACGGCAGGCAGCCTTCag                                                                                                                                            <  1:819378/148‑1 (MQ=255)
     cccATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCagg                                                                                                                                           >  1:39515/1‑149 (MQ=255)
      ccATTGCCAGTTTATTCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACTTAGAACCTTACCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTTCGGCAATATTGATCTCCTGACGGCAGGCAGCCTTCagga                                                                                                                                          <  2:102958/149‑1 (MQ=255)
       ccttGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGACCCTTCCCACTGCAACTGTTTGGTTTCCATTACAACTCCCGTACTCGCTCTCCCTCCGGCAAGATTGCTCTCCTGACGGCTGCCTGCCGTCcggag                                                                                                                                         >  2:213836/3‑149 (MQ=255)
           gCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGgtgt                                                                                                                                     <  2:150306/149‑1 (MQ=255)
            ccAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGACGGCAGGCAGCCTTCaggagatcaa                                                                                                                                    >  1:288404/1‑144 (MQ=255)
             cAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACATACCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTgg                                                                                                                                   >  1:270605/1‑149 (MQ=255)
                    tGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAAGCTTCCCTCTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCTCTCGACCTACTGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCATcc                                                                                                                            <  2:1143740/149‑1 (MQ=255)
                    tGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAcc                                                                                                                            <  2:1159987/149‑1 (MQ=255)
                    cgCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAcc                                                                                                                            <  2:1143739/148‑1 (MQ=255)
                     gCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAccc                                                                                                                           <  1:501092/149‑1 (MQ=255)
                        tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAGCCCtat                                                                                                                        >  1:32075/1‑149 (MQ=255)
                        tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCtat                                                                                                                        >  1:274864/1‑149 (MQ=255)
                        tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCtat                                                                                                                        <  2:418485/149‑1 (MQ=255)
                        tCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCtat                                                                                                                        >  1:1109831/1‑149 (MQ=255)
                                tGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATgcgc                                                                                                                <  1:37166/149‑1 (MQ=255)
                                tGAAGGCTGGAACTCCCCTGAACGTAGAACCATCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATgcgc                                                                                                                <  2:873556/149‑1 (MQ=255)
                                        ggAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGt                                                                                                        >  1:518393/1‑149 (MQ=255)
                                             tCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATATCCTGAAGGCTGCCTGCCGTCAGTAGGTGTGGATCAACACTATAGATGCGCAGAAACGTGGGAt                                                                                                   >  1:510064/1‑149 (MQ=255)
                                              gcccTGAACTTAGAACCATCGATCTTCTACTGTTTGCTGTCAAAAACACTTCCCGTACTCACTCGACCTACGGCACTATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCATATAGATGCGCAGAAACGTGGGAtt                                                                                                  <  2:1096942/148‑1 (MQ=255)
                                              ccccTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGAtt                                                                                                  <  2:88854/149‑1 (MQ=255)
                                               ccttGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTAAAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCTGAAACGTGGGTTTg                                                                                                 <  2:52714/146‑1 (MQ=255)
                                               cccTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATATCCTGAAGGCTGCCTGCCGGCAGGAGGTGTGGAGCAACCCTATAGATGCGCAGAAACGTGGGATTg                                                                                                 >  1:378415/1‑149 (MQ=255)
                                                 cTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGcc                                                                                               >  1:734799/1‑149 (MQ=255)
                                                 cTGAACCTAGAACCTTCCCACTGCAACTGTATGGTTTCCATTACAAATAACGTACTCACTCGCCCTACGGCAATATTGATCTCCTGAAGGCTGCCTCCCGTCAGGAGGTGTCGATCACCCCTATAGATGCGCAGAAACGTGGGATTTcc                                                                                               >  2:1117998/1‑149 (MQ=255)
                                                   gAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTCCTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCaa                                                                                             <  2:553610/149‑1 (MQ=255)
                                                      cGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCAGTACTCACTCGACCTACGGCAATATTGATCTGCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAAAGTGGGATTGCCAACgg                                                                                          <  2:1023438/149‑1 (MQ=255)
                                                        gaGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTCCTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGCTGCGCAGAAACGTGGGATTGCCAACGGCg                                                                                        <  2:337365/148‑1 (MQ=255)
                                                             ccTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAAGCCTATAGCTGCGCAGAAACGTGGTAATGCCAACGGCGATATg                                                                                   >  2:1045758/1‑149 (MQ=255)
                                                             ccTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACCTGGGATTGCCAACGGCGATATg                                                                                   >  1:100593/1‑149 (MQ=255)
                                                                 cccACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGc                                                                               >  2:820212/1‑149 (MQ=255)
                                                                    aCTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCggg                                                                            >  1:335560/1‑149 (MQ=255)
                                                                        cAACTGTTTGGTTTCCATTACAAATCCCGTCCGCCATCGACCTCCGGCAATATTGATCTCCTGCAGGCTGCCTGCCGTCAGGGGGGGTGGATCCCCCCTATAGATGCGCCGCAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                        >  1:893555/1‑149 (MQ=255)
                                                                        cAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGtt                                                                        <  2:270605/149‑1 (MQ=255)
                                                                            tGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAAc                                                                    >  1:304651/1‑149 (MQ=255)
                                                                                         ttACAAATCCCGTACTCACTCCACCTACGGCAATATTGATCTCCTGAAGGCTGCCTCCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGCCGATATTGTGCGGGTCTTTAACCATCGTGGCGAAg                                                       >  1:121434/1‑149 (MQ=255)
                                                                                               aTCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGACGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGc                                                 <  1:131717/149‑1 (MQ=255)
                                                                                               aTCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGc                                                 <  1:217101/149‑1 (MQ=255)
                                                                                                tCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACACTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGGGCGGGTGTTTACCCATCGTGGCGAAGTTCGGCt                                                >  1:796084/1‑149 (MQ=255)
                                                                                                  ccGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCCACCCTATAGATGCGCAGAACCGTGGGAATGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTAc                                              >  2:687637/1‑149 (MQ=255)
                                                                                                    gaagtcactcGTCCTACGGCAATTTTGATCTCCTGATGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGCTTTTGTCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCa                                            <  2:438861/145‑1 (MQ=255)
                                                                                                    gTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCa                                            >  1:152291/1‑149 (MQ=255)
                                                                                                     tACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAg                                           <  1:870504/149‑1 (MQ=255)
                                                                                                          actcGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCaaa                                      >  1:1024060/1‑149 (MQ=255)
                                                                                                           ctcGACCTACGGCAATAATGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAg                                     <  2:39515/149‑1 (MQ=255)
                                                                                                              gACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTaa                                  <  1:820212/149‑1 (MQ=255)
                                                                                                               aCCTACCGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTCCCAACGGCGATATGGTGCGGGTGTTTAAACATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAac                                 <  2:510064/149‑1 (MQ=255)
                                                                                                                  tACGGCAATATTGATCTCCTGAAGGCTGCATGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGACAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAAcacc                              >  2:71024/1‑149 (MQ=255)
                                                                                                                        aaTATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTAt                        >  2:137870/1‑149 (MQ=255)
                                                                                                                         atatTGATCTCCTGAAGTCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGGGTTTAACCATCGTTGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTAtt                       <  1:298173/149‑1 (MQ=255)
                                                                                                                             tGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTcc                   >  2:1133794/1‑149 (MQ=255)
                                                                                                                                 ctcCTGACGGCAGGCAGCCTTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGAATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTgg               >  1:282754/1‑149 (MQ=255)
                                                                                                                                 ctcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAAATGAGCATAAACGTGGGATTGCCAACGTCCATATGGTGCGGGTGTTTAACAATCGAGGCCAAGTTCGGCTACCAGACAAATGAACACCACGTATTCTCCCtag               >  2:495903/1‑147 (MQ=255)
                                                                                                                                  tcCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGa              >  2:15319/1‑149 (MQ=255)
                                                                                                                                    cTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGt            >  1:416876/1‑149 (MQ=255)
                                                                                                                                    cTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGt            >  1:1105143/1‑149 (MQ=255)
                                                                                                                                      gAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGCAGttt          >  2:231090/1‑148 (MQ=255)
                                                                                                                                             gcctgccGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTAt   >  1:837527/1‑149 (MQ=255)
                                                                                                                                              cctgccGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCCCGTATTCTCCCTGGAGTTAGCGCTATg  >  2:855857/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATAAGCCCATTGCCAGTTTATGCCTCAACCTTTGAAGGCTGGAACTCCCCTGAACGTAGAACCTTCCCACTGCAACTGTTTGGTTTCCATTACAAATCCCGTACTCACTCGACCTACGGCAATATTGATCTCCTGAAGGCTGCCTGCCGTCAGGAGGTGTGGATCAACCCTATAGATGCGCAGAAACGTGGGATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATG  >  NC_000913/1660058‑1660348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑