Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,729,467 T→C 7.0% intergenic (‑23/+149) gltW ← / ← rrsG tRNA‑Glu/16S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,729,4670TC 7.0% 228.8 / 7.8 86intergenic (‑23/+149)gltW/rrsGtRNA‑Glu/16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base C (3/3);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

TGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGG  >  NC_000913/2729322‑2729615
                                                                                                                                                 |                                                                                                                                                    
tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCGTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGt                                                                                                                                                  <  2:905068/150‑1 (MQ=17)
tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGt                                                                                                                                                  <  2:1379661/150‑1 (MQ=32)
tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGt                                                                                                                                                  <  1:1403230/150‑1 (MQ=32)
  ttttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGa                                                                                                                                                <  2:1007097/150‑1 (MQ=32)
    ttGATATATTGAGGTCGGCGACTTTCACTCACAAACCAGCTAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGTCGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGata                                                                                                                                              <  1:667576/150‑1 (MQ=12)
    ttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGata                                                                                                                                              <  1:1219552/150‑1 (MQ=12)
     tGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGatat                                                                                                                                             <  1:507194/150‑1 (MQ=32)
     ctttgtagtgctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  1:360210‑M1/117‑150 (MQ=255)
     ctttgcagtgctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  1:735262‑M1/117‑150 (MQ=255)
     ctttgcagggctcacaccgattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  2:1496206‑M1/117‑150 (MQ=255)
      gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGataag                                                                                                                                            <  2:206517/150‑3 (MQ=12)
       agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGataagc                                                                                                                                           <  1:672643/150‑4 (MQ=255)
          tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  2:893811/1‑150 (MQ=32)
          tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  1:1684660/1‑150 (MQ=32)
          tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  1:1595734/1‑150 (MQ=32)
          tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                        <  1:1238324/150‑1 (MQ=32)
            ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgtg                                                                                                                                      >  2:777276/1‑150 (MQ=32)
            gtgctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgtg                                                                                                                                      >  1:1128288‑M1/110‑150 (MQ=255)
             tgctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgtgt                                                                                                                                     >  2:613065‑M1/109‑150 (MQ=255)
             tGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgtgt                                                                                                                                     >  1:1068390/1‑150 (MQ=31)
              gAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTa                                                                                                                                    >  2:1194482/1‑150 (MQ=11)
               ctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAg                                                                                                                                   >  1:589672‑M1/107‑150 (MQ=255)
               ctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAg                                                                                                                                   >  1:1373563‑M1/107‑150 (MQ=255)
               ctcacacagatagtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAg                                                                                                                                   >  1:631182‑M1/107‑150 (MQ=255)
                 cacacagattggctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataacgccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGGGTAggg                                                                                                                                 >  2:1515091‑M1/105‑150 (MQ=255)
                      cgtcTTTCACTCACAAACCCTCAAGTGGCGTCACCTTGGTGATTCGAACCCCTGTTAACGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCTTTGCGTGATATTTTCGTGTAGGGTGAGc                                                                                                                            <  2:1702600/147‑1 (MQ=1)
                                                                                          gggCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                                                        >  1:1436891/1‑150 (MQ=11)
                                                                                           ggCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAt                                                       >  2:510231/1‑150 (MQ=31)
                                                                                            gCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                                                      >  2:403584/1‑150 (MQ=32)
                                                                                            gCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                                                      <  1:847277/150‑1 (MQ=32)
                                                                                             cGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTACACGCCTTGCTATTc                                                     >  2:696335/1‑150 (MQ=17)
                                                                                             cGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTc                                                     >  2:690295/1‑150 (MQ=32)
                                                                                              ggTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                                                    <  1:111669/150‑1 (MQ=32)
                                                                                                tgtCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                                                  <  2:1180818/150‑1 (MQ=32)
                                                                                                   ccTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                               >  2:720085/1‑150 (MQ=32)
                                                                                                     tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCAGTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca                                             >  1:907923/1‑150 (MQ=21)
                                                                                                      gggCCTCTAGACGATGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTAGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                            <  1:293178/150‑1 (MQ=25)
                                                                                                      gggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                            <  1:1444116/150‑1 (MQ=35)
                                                                                                      gggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                            <  2:513294/150‑1 (MQ=35)
                                                                                                       ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc                                           <  2:395767/150‑1 (MQ=35)
                                                                                                        gCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca                                          <  2:909442/150‑1 (MQ=37)
                                                                                                         ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg                                         >  2:1605948/1‑150 (MQ=37)
                                                                                                         ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg                                         >  1:8717/1‑150 (MQ=37)
                                                                                                          ctctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                        >  1:147990/1‑150 (MQ=37)
                                                                                                           tctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGAc                                       <  1:316395/150‑1 (MQ=37)
                                                                                                            ctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACa                                      >  2:1323227/1‑150 (MQ=37)
                                                                                                              aGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                    >  2:1092535/1‑150 (MQ=37)
                                                                                                              aGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                    >  1:683477/1‑150 (MQ=37)
                                                                                                              aGACCAACGGGACACGAAAATTCCTTATCACGCGTTGCGTGATATTTTCGTGTTGGGAGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACACGTAAACGCCTTGCAATTCACATTTGATCAGACAAt                                    >  2:66074/1‑150 (MQ=11)
                                                                                                                  gAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                                <  1:1192014/150‑1 (MQ=37)
                                                                                                                   aaGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtg                               <  1:159641/150‑1 (MQ=37)
                                                                                                                    aGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                              <  1:634366/150‑1 (MQ=37)
                                                                                                                    aGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                              >  1:483575/1‑150 (MQ=37)
                                                                                                                    aGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTATCTTCGGCCTCACTCCCAACGCGTAAACGCCTGGCGTTTCACGTTTCATCAGACAATCtgtgt                              >  2:55649/1‑150 (MQ=11)
                                                                                                                     ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgtg                             <  1:553741/150‑1 (MQ=37)
                                                                                                                     ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgtg                             <  2:1130017/150‑1 (MQ=37)
                                                                                                                       ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAg                           >  2:267389/1‑150 (MQ=37)
                                                                                                                        gACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGc                          >  1:958668/1‑150 (MQ=37)
                                                                                                                        gACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGc                          >  1:849219/1‑150 (MQ=37)
                                                                                                                          CACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                        <  1:1411263‑M1/150‑1 (MQ=255)
                                                                                                                          CACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                        <  1:1056313‑M1/150‑1 (MQ=255)
                                                                                                                          CACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCACAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                        <  2:432179‑M1/150‑1 (MQ=255)
                                                                                                                            CGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACTCGTAAACGCCATGCTATTCACTTTTCATCAGACAATCTCTGTGAGCACTg                      <  2:725290‑M1/150‑1 (MQ=255)
                                                                                                                            CGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTg                      <  2:303679‑M1/150‑1 (MQ=255)
                                                                                                                            CGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTg                      <  2:1486167‑M1/150‑1 (MQ=255)
                                                                                                                              AAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCa                    >  2:585011‑M1/1‑150 (MQ=255)
                                                                                                                                AATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                  >  2:869690‑M1/1‑150 (MQ=255)
                                                                                                                                AATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACACCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                  <  1:524124‑M1/150‑1 (MQ=255)
                                                                                                                                   aaaatATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTa               >  2:1660233/5‑150 (MQ=255)
                                                                                                                                    tattATCACGCGTTGCGTTATATTTTCGTGTAGGGTGAGCTTTCTTTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTTAGCACTGCAAAGTAc              <  2:1316240/148‑1 (MQ=255)
                                                                                                                                    gCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTAc              <  2:683477/150‑1 (MQ=39)
                                                                                                                                     cTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACg             >  1:965279/1‑150 (MQ=39)
                                                                                                                                      ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGc            <  2:1068390/150‑1 (MQ=39)
                                                                                                                                      atATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGc            >  1:521312/2‑150 (MQ=255)
                                                                                                                                      atATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGc            >  2:1458135/2‑150 (MQ=255)
                                                                                                                                       tATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgtg                             >  2:1302119/1‑132 (MQ=37)
                                                                                                                                       tATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgtg                             <  1:1302119/132‑1 (MQ=37)
                                                                                                                                       tATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGct           <  2:1336230/150‑1 (MQ=39)
                                                                                                                                        aTCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGctt          >  1:1251929/1‑150 (MQ=39)
                                                                                                                                          cACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct        <  2:259687/150‑1 (MQ=39)
                                                                                                                                          cACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct        <  2:507980/150‑1 (MQ=39)
                                                                                                                                           aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt       >  1:1399365/1‑150 (MQ=39)
                                                                                                                                           aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt       >  1:1637020/1‑150 (MQ=39)
                                                                                                                                            cgcgTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCttt      <  2:922706/150‑1 (MQ=39)
                                                                                                                                              cgTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaa    <  2:1411263/150‑1 (MQ=39)
                                                                                                                                               gTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaag   <  1:51658/150‑1 (MQ=39)
                                                                                                                                                ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaagg  <  1:510231/150‑1 (MQ=39)
                                                                                                                                                ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaagg  <  1:411815/150‑1 (MQ=39)
                                                                                                                                                 |                                                                                                                                                    
TGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGG  >  NC_000913/2729322‑2729615

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