Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,528,705 T→C 8.4% V349A (GTA→GCA)  ftsA → cell division protein FtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,528,7050TC 8.4% 200.2 / 6.7 71V349A (GTA→GCA) ftsAcell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base T (29/36);  new base C (3/3);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GACGAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAC  >  NC_002947/1528563‑1528841
                                                                                                                                              |                                                                                                                                        
gacgaGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCt                                                                                                                                    >  6:57404/1‑149 (MQ=255)
   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                                                                                                                         >  1:121038/1‑145 (MQ=255)
              ccTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGCACGCCCAGGCGTACCgg                                                                                                                         <  2:121038/146‑2 (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                                                                                                                      >  6:32736/1‑139 (MQ=255)
                        gCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGcaca                                                                                                                      <  5:32736/139‑1 (MQ=255)
                          cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                    >  7:273720/1‑139 (MQ=255)
                          cGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGc                                                                                                                    <  8:273720/139‑1 (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                                       <  6:55947/131‑1 (MQ=255)
                                                                                                               gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGaa                                       >  5:55947/1‑131 (MQ=255)
                                                                                                               gAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGc                     >  8:104669/1‑149 (MQ=255)
                                                                                                                 aCTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGccc                   <  7:104669/149‑1 (MQ=255)
                                                                                                                     gCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGTAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGa                <  3:308111/148‑1 (MQ=255)
                                                                                                                       cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg             >  1:311953/1‑149 (MQ=255)
                                                                                                                       cGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCg             >  4:149490/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  >  1:97568/6‑145 (MQ=255)
                                                                                                                                      ctgtgcGGCACGCCCAGGCGTACCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTaacagcaac  <  2:97568/140‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
GACGAGCTCTTCACCCTTATCCAGGCCGAACTGCGCCGCAGTGGCTACGAAGACCTGGTGCCCGCCGGCATCGTCCTCACCGGCGGCACCGCGAAGATGGAAGGCGCCGTCGAACTGGCCGAGGAAATCTTCCACATGCCGGTACGCCTGGGCGTGCCGCACAGCGTTCGGGGGCTGAGCGACGTGGTGCGCAACCCGATCTATTCCACTGGTGTGGGCTTGCTCACCTACGGCCTGCTGAAGCAGTCCGAGGACTTGGCCCTGACCGGTAACAGCAAC  >  NC_002947/1528563‑1528841

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: ‑