Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,729,466 T→A 6.5% 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 6.5% 293.5 / 4.9 87intergenic (‑22/+150)gltW/rrsGtRNA‑Glu/16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base T (44/37);  new base A (4/2);  total (48/39)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAG  >  NC_000913/2729318‑2729614
                                                                                                                                                    |                                                                                                                                                    
aGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGACCCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGt                                                                                                                                                      <  4:1329028/149‑1 (MQ=25)
 gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGtt                                                                                                                                                     <  3:306517/149‑1 (MQ=35)
 gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGtt                                                                                                                                                     >  3:845076/1‑149 (MQ=35)
 gTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGAGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGtt                                                                                                                                                     <  3:286478/149‑1 (MQ=25)
  tCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTg                                                                                                                                                    >  1:408443/1‑149 (MQ=35)
   cTGTTTTGAGAGATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGc                                                                                                                                                   <  2:967285/149‑1 (MQ=18)
    tGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCg                                                                                                                                                  <  3:68224/149‑1 (MQ=32)
     gTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGt                                                                                                                                                 <  3:1367468/149‑1 (MQ=32)
      ttttGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTg                                                                                                                                                >  3:541005/1‑149 (MQ=32)
          gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGatat                                                                                                                                            >  1:951328/1‑149 (MQ=32)
          gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGatat                                                                                                                                            >  2:184079/1‑149 (MQ=32)
          gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGatat                                                                                                                                            <  1:1377507/149‑1 (MQ=32)
          gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGataa                                                                                                                                            <  2:244443/149‑2 (MQ=12)
           agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATAtt                                                                                                                                           >  2:562057/1‑149 (MQ=32)
           agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGat                                                                                                                                              <  1:1019419/146‑1 (MQ=12)
           agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGat                                                                                                                                              >  2:1019419/1‑146 (MQ=12)
            gaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCAACGCGTGataaga                                                                                                                                          >  2:258732/1‑146 (MQ=255)
              tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                        <  3:1044008/149‑1 (MQ=32)
              cagtgctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                        >  3:847982‑M1/112‑149 (MQ=255)
               atTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                       >  4:1372994/1‑149 (MQ=32)
               atTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                       >  2:985776/1‑149 (MQ=32)
               atTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                       >  2:1041460/1‑149 (MQ=32)
               atTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCg                                                                                                                                       >  2:1412841/1‑149 (MQ=32)
                ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgt                                                                                                                                      <  2:996979/149‑1 (MQ=32)
                 tgctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCgtg                                                                                                                                     >  3:498883‑M1/109‑149 (MQ=255)
                  gctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                        <  3:77859‑M1/38‑1 (MQ=255)
                  gctcacacagattgtctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTc                                                                                                                                        >  4:77859‑M1/108‑145 (MQ=255)
                   ctcacacagattgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcacccgACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTa                                                                                                                                   >  3:766030‑M1/107‑149 (MQ=255)
                             ttgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTa                                                                                                                                   >  1:179804‑M1/97‑139 (MQ=255)
                             ttgtctgatagaaagtgaaaagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTa                                                                                                                                   <  2:179804‑M1/43‑1 (MQ=255)
                                    atgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcacccgACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTaa                                                                                                                  >  1:1329517‑M1/90‑149 (MQ=255)
                                        aaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCtt                                                                                                                         <  2:1380807‑M1/53‑1 (MQ=255)
                                        aaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCtt                                                                                                                         >  1:1380807‑M1/86‑138 (MQ=255)
                                                                                                 cGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                                                     <  1:928406/149‑1 (MQ=32)
                                                                                                 cGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                                                     <  1:656022/149‑1 (MQ=32)
                                                                                                     gtCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                                                 >  3:435083/1‑149 (MQ=32)
                                                                                                     gtCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                                                 <  1:562057/149‑1 (MQ=32)
                                                                                                      tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                                >  1:655557/1‑149 (MQ=32)
                                                                                                      tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                                >  2:266953/1‑149 (MQ=32)
                                                                                                      tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                                >  4:218474/1‑149 (MQ=32)
                                                                                                        cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                              >  2:1372724/1‑149 (MQ=32)
                                                                                                        cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                              >  4:88477/1‑149 (MQ=33)
                                                                                                        cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                              >  4:98497/1‑149 (MQ=33)
                                                                                                        cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttct                                              >  1:1082639/1‑147 (MQ=17)
                                                                                                         tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                             >  2:1123688/1‑149 (MQ=33)
                                                                                                         tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                             <  3:908047/149‑1 (MQ=33)
                                                                                                         tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                             >  2:536667/1‑149 (MQ=34)
                                                                                                          gggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca                                            >  2:308680/1‑149 (MQ=34)
                                                                                                           ggCCTCTAGACGAAGGGGACACGAAAGTTGCTTATCACGCGTTGCGTGATAGTTTCGTGTAGGGTGAGCTTGCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                           >  3:1139485/1‑149 (MQ=255)
                                                                                                           ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                           <  3:98497/149‑1 (MQ=35)
                                                                                                           ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                           <  2:1369690/149‑1 (MQ=35)
                                                                                                           ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                           <  1:1372724/149‑1 (MQ=35)
                                                                                                           ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                           <  2:28694/149‑1 (MQ=35)
                                                                                                            gCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc                                          >  3:1191012/1‑149 (MQ=35)
                                                                                                            gCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc                                          <  2:149273/149‑1 (MQ=35)
                                                                                                            gCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc                                          >  1:1363862/1‑149 (MQ=35)
                                                                                                               tctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                       >  3:592709/1‑149 (MQ=37)
                                                                                                                   gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                   >  3:658233/1‑149 (MQ=38)
                                                                                                                     cGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt                                 <  1:1457221/149‑1 (MQ=37)
                                                                                                                      gAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg                                >  3:597541/1‑149 (MQ=38)
                                                                                                                       aaGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                               <  3:194739/149‑1 (MQ=38)
                                                                                                                       aaGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                               <  4:271789/149‑1 (MQ=38)
                                                                                                                        aGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtg                              <  3:319027/149‑1 (MQ=38)
                                                                                                                         ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                             >  4:199184/1‑149 (MQ=38)
                                                                                                                         ggggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                             <  4:106243/149‑1 (MQ=39)
                                                                                                                          gggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgtg                            >  4:902523/1‑149 (MQ=37)
                                                                                                                           ggACACGAAAATTGCTTATCACTCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                           <  4:789074/149‑1 (MQ=37)
                                                                                                                           ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                           <  4:847778/149‑1 (MQ=38)
                                                                                                                           ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                           <  1:690860/149‑1 (MQ=38)
                                                                                                                           ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                           <  1:359357/149‑1 (MQ=38)
                                                                                                                           ggACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                           <  1:1239304/149‑1 (MQ=39)
                                                                                                                            gACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAg                          >  2:996769/1‑149 (MQ=39)
                                                                                                                             ACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGc                         <  2:979141‑M1/149‑1 (MQ=255)
                                                                                                                               ACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                       >  1:184759‑M1/1‑149 (MQ=255)
                                                                                                                               ACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                       <  2:706580‑M1/149‑1 (MQ=255)
                                                                                                                                     ATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                 >  2:1071518‑M1/1‑149 (MQ=255)
                                                                                                                                     ATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                 >  3:199184‑M1/1‑149 (MQ=255)
                                                                                                                                      TTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAg                <  3:899767‑M1/149‑1 (MQ=255)
                                                                                                                                        gCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTa              <  1:996769/149‑1 (MQ=255)
                                                                                                                                         tatATCACGCCACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCATCACTCCCAAAGCGTAAACGCCTTGCTATTCACTTTTCATAAGACAATCTGTGTGAGCACTGCAAAGTAc             >  2:73060/3‑149 (MQ=14)
                                                                                                                                         cTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTAc             <  1:767517/149‑1 (MQ=40)
                                                                                                                                         cTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTAc             <  3:1372994/149‑1 (MQ=40)
                                                                                                                                         aatATCACGCAACGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTAc             >  2:891261/3‑149 (MQ=255)
                                                                                                                                          ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACg            >  4:594697/1‑149 (MQ=39)
                                                                                                                                          ttATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACg            >  3:1394390/1‑149 (MQ=255)
                                                                                                                                           tATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGc           <  4:1036809/149‑1 (MQ=40)
                                                                                                                                               aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct       >  1:43350/1‑149 (MQ=255)
                                                                                                                                               aCGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttct       >  2:465148/1‑149 (MQ=40)
                                                                                                                                                cgcgTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt      <  3:1447463/149‑1 (MQ=39)
                                                                                                                                                 gcgTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCttt     >  4:887431/1‑149 (MQ=255)
                                                                                                                                                   gTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaa   <  1:1401225/149‑1 (MQ=40)
                                                                                                                                                    ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaag  <  3:902523/149‑1 (MQ=40)
                                                                                                                                                    ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaag  <  3:757877/149‑1 (MQ=39)
                                                                                                                                                    ttGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaag  <  1:1267294/149‑1 (MQ=40)
                                                                                                                                                    |                                                                                                                                                    
AGTCTGTTTTGAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAG  >  NC_000913/2729318‑2729614

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