Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,729,466 T→A 7.0% intergenic (‑22/+150) 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,4660TA 7.0% 204.6 / 7.1 86intergenic (‑22/+150)gltW/rrsGtRNA‑Glu/16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (41/39);  new base A (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 = 7.04e-02

CAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGG  >  NC_000913/2729317‑2729615
                                                                                                                                                     |                                                                                                                                                     
cAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGTAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGt                                                                                                                                                       <  1:431583/150‑1 (MQ=37)
     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:1238324/150‑1 (MQ=32)
               tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  1:1595734/1‑150 (MQ=32)
               tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                        >  1:1684660/1‑150 (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                                                      <  1:847277/150‑1 (MQ=32)
                                                                                                 gCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                                                      >  2:403584/1‑150 (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                                            <  2:513294/150‑1 (MQ=35)
                                                                                                           gggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                            <  1:1444116/150‑1 (MQ=35)
                                                                                                            ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc                                           <  2:395767/150‑1 (MQ=35)
                                                                                                             gCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca                                          <  2:909442/150‑1 (MQ=37)
                                                                                                              ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg                                         >  1:8717/1‑150 (MQ=37)
                                                                                                              ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg                                         >  2:1605948/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                             <  2:1130017/150‑1 (MQ=37)
                                                                                                                          ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgtg                             <  1:553741/150‑1 (MQ=37)
                                                                                                                            ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAg                           >  2:267389/1‑150 (MQ=37)
                                                                                                                             gACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGc                          >  1:849219/1‑150 (MQ=37)
                                                                                                                             gACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGc                          >  1:958668/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:1486167‑M1/150‑1 (MQ=255)
                                                                                                                                 CGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTg                      <  2:303679‑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            >  2:1458135/2‑150 (MQ=255)
                                                                                                                                           atATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGc            >  1:521312/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:507980/150‑1 (MQ=39)
                                                                                                                                               cACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct        <  2:259687/150‑1 (MQ=39)
                                                                                                                                                aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt       >  1:1637020/1‑150 (MQ=39)
                                                                                                                                                aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt       >  1:1399365/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)
                                                                                                                                                     |                                                                                                                                                     
CAGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGG  >  NC_000913/2729317‑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: