Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,528,712 G→A 8.6% L351L (CTG→CTA) ‡ ftsA → cell division protein FtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,528,7120GA 8.6% 192.6 / 4.5 69L351L (CTG→CTA) ‡ftsAcell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base G (28/35);  new base A (3/3);  total (31/38)
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.59e-01

GAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTA  >  NC_002947/1528566‑1528858
                                                                                                                                                  |                                                                                                                                                  
gaGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTggg                                                                                                                                                  <  3:136341/149‑1 (MQ=255)
 aGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGc                                                                                                                                                 <  5:57404/149‑1 (MQ=255)
 aGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGc                                                                                                                                                 <  7:216325/149‑1 (MQ=255)
  gCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCg                                                                                                                                                >  5:210023/1‑149 (MQ=255)
    tctTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTg                                                                                                                                              >  7:181243/1‑149 (MQ=255)
    tctTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTCCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTg                                                                                                                                              >  7:4663/1‑149 (MQ=255)
         aCCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGca                                                                                                                                         >  5:16559/1‑149 (MQ=255)
           ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTgg                                                                                                                                                   >  8:59935/1‑137 (MQ=255)
           ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTgg                                                                                                                                                   <  7:59935/137‑1 (MQ=255)
           ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                       <  3:276843/149‑1 (MQ=255)
           ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCgg                                                                                                                                          <  2:121038/146‑2 (MQ=255)
           ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCgg                                                                                                                                          >  1:121038/1‑145 (MQ=255)
           cattATCCAGGACGAACTGCGCCGCAGTGGCTACGAAGAAATGGTGCCCGCCGGCATCGTCATCACCGGCGGCACCGCGAAGATGGAAGGCCCCGTCGAACTGCACGAGGAAATCTTCAACATGCCGCTACGCCTGCGAGTGCCGca                                                                                                                                         >  4:166110/3‑147 (MQ=255)
             ttATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                     >  4:226223/1‑149 (MQ=255)
                   aGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCggc                                                                                                                                         <  4:82630/139‑3 (MQ=255)
                   aGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCggc                                                                                                                                         >  3:82630/1‑137 (MQ=255)
                    ggCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCggg                                                                                                                              >  8:255135/1‑149 (MQ=255)
                    ggCCGAACTGCGCCGCAGTGGCGACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGGGCCGCACAGCGTTCggg                                                                                                                              <  6:324993/149‑1 (MQ=255)
                     gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                       <  5:32736/139‑1 (MQ=255)
                     gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                                       >  6:32736/1‑139 (MQ=255)
                       cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                     <  8:273720/139‑1 (MQ=255)
                       cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                                     >  7:273720/1‑139 (MQ=255)
                                ccgcAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACg                                                                                                                  >  7:284141/1‑149 (MQ=255)
                                     gTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtggtg                                                                                                             <  8:181243/149‑1 (MQ=255)
                                      tGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                            <  8:112741/149‑1 (MQ=255)
                                       ggCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgc                                                                                                            <  3:106971/148‑1 (MQ=255)
                                          tACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGCCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCaa                                                                                                        >  6:168449/1‑149 (MQ=255)
                                                gACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGAt                                                                                                  <  3:2663/149‑1 (MQ=255)
                                                    tGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACgtg                                                                                                                >  2:28036/1‑131 (MQ=255)
                                                    tGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACgtg                                                                                                                <  1:28036/131‑1 (MQ=255)
                                                      gTGCCCGCCGGCATCGGCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTAt                                                                                              <  5:41287/147‑1 (MQ=255)
                                                        gCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCa                                                                                          <  8:284141/149‑1 (MQ=255)
                                                           cgccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTg                                                                                       >  5:82444/1‑149 (MQ=255)
                                                           cgccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTg                                                                                       >  1:247152/1‑149 (MQ=255)
                                                             ccgGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGgt                                                                                     <  6:210023/149‑1 (MQ=255)
                                                                gCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGtgtg                                                                                  <  1:293789/149‑1 (MQ=255)
                                                                gCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGtgtg                                                                                  <  3:226223/149‑1 (MQ=255)
                                                                   tCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgcg                                                                                                           <  8:80761/121‑1 (MQ=255)
                                                                   tCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTgcg                                                                                                           >  7:80761/1‑121 (MQ=255)
                                                                     gTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCtt                                                                             >  8:15779/1‑149 (MQ=255)
                                                                         tccccGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTc                                                                         >  5:265366/4‑149 (MQ=255)
                                                                           aCCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTc                                                                         <  4:103994/147‑1 (MQ=255)
                                                                                 ggcACCGCGAAGATGGAAGGCGCCGTCGAACTGGCTGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACgg                                                                 <  8:96095/149‑1 (MQ=255)
                                                                                    aCCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCct                                                              <  6:82444/149‑1 (MQ=255)
                                                                                     gcgcgAAGATGGAATGCTCCGTCGAACTGGCCGAGGAAATCTTGCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGTGTGAGCGACGTGGTGCGCATCCCGATGTATTCAACTGGTTTGGGCTTGCTCACCTACGGCctg                                                             <  3:166110/148‑1 (MQ=255)
                                                                                      cgcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgc                                                            >  5:158672/1‑149 (MQ=255)
                                                                                       gcgAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGTAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgct                                                           >  4:308111/1‑149 (MQ=255)
                                                                                         gAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCctgctg                                                          >  3:30677/1‑148 (MQ=255)
                                                                                         gAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGa                                                         <  4:30677/149‑1 (MQ=255)
                                                                                          aaGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                        <  7:255135/149‑1 (MQ=255)
                                                                                                gAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGt                                                   <  2:221569/148‑1 (MQ=255)
                                                                                                  aGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTcc                                                 >  1:308932/1‑148 (MQ=255)
                                                                                                    gcgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAg                                              <  6:16559/149‑1 (MQ=255)
                                                                                                     cgcCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTTCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAgg                                             <  3:200335/149‑1 (MQ=255)
                                                                                                       ccGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGAc                                           <  7:15779/149‑1 (MQ=255)
                                                                                                            gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                        >  5:55947/1‑131 (MQ=255)
                                                                                                            gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                                        <  6:55947/131‑1 (MQ=255)
                                                                                                            gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGc                                      >  8:104669/1‑149 (MQ=255)
                                                                                                              aCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGccc                                    <  7:104669/149‑1 (MQ=255)
                                                                                                                  gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGTAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGa                                 <  3:308111/148‑1 (MQ=255)
                                                                                                                    cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                              >  4:149490/1‑149 (MQ=255)
                                                                                                                    cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg                              >  1:311953/1‑149 (MQ=255)
                                                                                                                     gAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                             <  2:308932/149‑1 (MQ=255)
                                                                                                                            tCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacag                       >  7:29208/1‑148 (MQ=255)
                                                                                                                                ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagc                      <  5:168449/145‑1 (MQ=255)
                                                                                                                                ccACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagc                      <  6:158672/145‑1 (MQ=255)
                                                                                                                                  acaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                             <  6:151056/136‑1 (MQ=255)
                                                                                                                                  acaTGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCgg                             >  5:151056/1‑136 (MQ=255)
                                                                                                                                   ctgtgcGGCACGCCCAGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac                   <  2:97568/140‑1 (MQ=255)
                                                                                                                                   ctgtgcGGCACGCCCAGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac                   >  1:97568/6‑145 (MQ=255)
                                                                                                                                             cGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACa      <  2:311953/148‑1 (MQ=255)
                                                                                                                                                cTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTa  >  3:188393/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAACAGCAGCAGCAACAGCTA  >  NC_002947/1528566‑1528858

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