Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 474,195 Δ1 bp 8.8% coding (793/1983 nt) dnaG ← DNA primase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947474,1950C. 8.8% 384.5 / 32.9 125coding (793/1983 nt)dnaGDNA primase
Reads supporting (aligned to +/- strand):  ref base C (55/57);  new base . (6/5);  total (62/63)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 131 total aligned reads displayed.

CCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGC  >  NC_002947/474053‑474343
                                                                                                                                              |                                                                                                                                                     
ccGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGAtc                                                                                                                                                >  2:3064377/1‑150 (MQ=255)
 cGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGAtct                                                                                                                                              >  1:2431013/1‑150 (MQ=255)
     aaaaCAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCtcgtc                                                                                                                                          >  1:2402502/1‑150 (MQ=255)
       aacagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGa                                                                                                                                        <  2:2725662/150‑1 (MQ=255)
         cagaacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATTATCTCGTCGAgg                                                                                                                                      <  1:660803/150‑1 (MQ=255)
          agaaacgCACGCGGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGt                                                                                                                                     >  2:1654023/7‑150 (MQ=255)
            aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTGGAGGTTg                                                                                                                                   >  1:942948/1‑150 (MQ=255)
            aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTg                                                                                                                                   <  2:2070613/150‑1 (MQ=255)
            aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGAGATCATCGTCGTCGa                                                                                                                                        >  1:2168046/1‑145 (MQ=255)
            aacagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGAGATCATCGTCGTCGa                                                                                                                                        <  2:2168046/145‑1 (MQ=255)
              cagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCg                                                                                                                                 <  1:2032609/150‑1 (MQ=255)
              cagCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCg                                                                                                                                 <  1:539327/150‑1 (MQ=255)
                 cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGtt                                                                                                                              <  1:1074288/150‑1 (MQ=255)
                 cACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGtt                                                                                                                              <  2:1221934/150‑1 (MQ=255)
                    gCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTtgtg                                                                                                                           >  2:1197074/1‑150 (MQ=255)
                      tGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATTACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCt                                                                                                                         >  1:896609/1‑150 (MQ=255)
                        ggCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCgg                                                                                                                                >  2:2459771/1‑141 (MQ=255)
                             cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                                  >  2:1285471/1‑150 (MQ=255)
                             cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                                  >  2:224694/1‑150 (MQ=255)
                             cacGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTgcgcgc                                                                                                                  >  1:1245424/1‑150 (MQ=255)
                                     aCAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCCTTGCGCTGCTGTTGATACTGCAGGGCACTGACGTCCATGTCGCCCTCGCCGCAGATGATCTCGTCGATGTTGCGGTTGTGCTTGCGCGCCTCGTACa                                                                                                          <  1:701599/150‑1 (MQ=255)
                                       aGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGc                                                                                                        <  1:2086306/150‑1 (MQ=255)
                                       aGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGAGATCATCGTCGTCGAgg                                                                                                                                      >  2:504224/1‑120 (MQ=255)
                                        ggCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGcc                                                                                                       <  1:224694/150‑1 (MQ=255)
                                         gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAg                                                                                                                                       >  1:943138/1‑117 (MQ=255)
                                         gcgcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                      <  1:1654023/150‑1 (MQ=255)
                                           gcTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGt                                                                                                    >  1:1327697/1‑150 (MQ=255)
                                            cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                      <  1:1143501/147‑1 (MQ=255)
                                            cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGccc                                                                                                      >  2:1143501/1‑147 (MQ=255)
                                            cTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTa                                                                                                   <  1:1720274/150‑1 (MQ=255)
                                                cggTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGAAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGc                                                                                               >  1:78992/2‑150 (MQ=255)
                                                  gTGCTCTTCGCTGGTGGCAGTGCCAAGGGTTGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTc                                                                                             >  1:70361/1‑150 (MQ=255)
                                                        ttCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCACGCCTCGTACAGCCGCAACAGCTCCTGccc                                                                                       >  2:177494/1‑150 (MQ=255)
                                                         tCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGcccc                                                                                      <  2:1245424/150‑1 (MQ=255)
                                                          cGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGCTGATCTCGTCGCGGTTGCGGTTGTGCTTGCGCCCCTCGTACAGCCCGTCCCGCTCCTGCCCCt                                                                                     >  2:1068118/1‑150 (MQ=255)
                                                             tggtggCGGTGCCAAGGGTGTCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  <  1:157816/150‑1 (MQ=255)
                                                             tggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  <  1:469467/150‑1 (MQ=255)
                                                             tggtggCGGTGCCAAGGGTGGCCACGGCATCGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTGGCGCGCCTCTTAGTGCCGGTACAGCTTCACCACCTtgt                                                                                  <  2:896609/150‑1 (MQ=255)
                                                              ggtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgtg                                                                                 <  1:2953270/150‑1 (MQ=255)
                                                               gtggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgt                                                                                  <  1:1713469/148‑1 (MQ=255)
                                                                tggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGtt                                                                                                                              <  2:1730037/103‑1 (MQ=255)
                                                                tggCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGtt                                                                                                                              >  1:1730037/1‑103 (MQ=255)
                                                                   cGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCCGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGaaaa                                                                            <  2:64048/150‑1 (MQ=255)
                                                                       gCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                        <  1:1842674/150‑1 (MQ=255)
                                                                        ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  2:1310926/150‑1 (MQ=255)
                                                                        ccAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCcgg                                                                       <  1:3064377/150‑1 (MQ=255)
                                                                         cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgtg                                                                                 >  1:1316696/1‑139 (MQ=255)
                                                                         cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGAGATCATCGTCGTCGAgg                                                                                                                                      <  1:2792945/86‑1 (MQ=38)
                                                                         cAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGAGATCATCGTCGTCGAgg                                                                                                                                      >  2:2792945/1‑86 (MQ=38)
                                                                                ggCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCgg                                                               <  1:1599083/150‑1 (MQ=255)
                                                                                 gCCACGGCATTGCGCCGGCCTAGCTGGGGCAAGGTAGCCACGTCTATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCggg                                                              <  2:1568726/150‑1 (MQ=255)
                                                                                   cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACATCTCCTGCCCCTtg                                                                                   <  2:755922/127‑1 (MQ=255)
                                                                                   cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtg                                                                                   >  1:755922/1‑127 (MQ=255)
                                                                                   cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                            <  1:1245652/150‑1 (MQ=255)
                                                                                   cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                            <  2:70361/150‑1 (MQ=255)
                                                                                   cACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTATGCTCCAGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGa                                                            >  2:1951175/1‑150 (MQ=255)
                                                                                         aTTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGGACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAg                                                      >  1:966924/1‑150 (MQ=255)
                                                                                              gcAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACt                                                 >  1:2578103/1‑150 (MQ=255)
                                                                                                aGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTcc                                                                                            <  1:1265514/105‑1 (MQ=255)
                                                                                                aGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTcc                                                                                            >  2:1265514/1‑105 (MQ=255)
                                                                                                aGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTg                                               >  2:1898798/1‑150 (MQ=255)
                                                                                                 ggCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTgg                                              >  2:2527309/1‑150 (MQ=255)
                                                                                                 ggCCTTGCTGGGCCAAGGCAATCACGTCCATGTAGCCTCCGACGACCATGATCTCCTCGAGGCTGCGGTTGTGCTTTCGCGCCTCGTACAGCCCGTACAGCTCCTGCCACTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTgg                                              <  2:942948/150‑1 (MQ=255)
                                                                                                  gCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGa                                                                                                                                        <  2:1832191/59‑1 (MQ=255)
                                                                                                  gCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCg                                                                                                     <  1:2295697/94‑1 (MQ=255)
                                                                                                  gCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCg                                                                                                     >  2:2295697/1‑94 (MQ=255)
                                                                                                        cTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGc                                                                                                        <  2:2631149/85‑1 (MQ=255)
                                                                                                        cTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGc                                                                                                        >  1:2631149/1‑85 (MQ=255)
                                                                                                            gCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGtcgtcg                                   >  2:850865/1‑150 (MQ=255)
                                                                                                             ccAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGc                                  <  1:1197074/150‑1 (MQ=255)
                                                                                                                gggCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGa                               <  2:2450516/150‑1 (MQ=255)
                                                                                                                gggCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGa                               <  2:2402502/150‑1 (MQ=255)
                                                                                                                gggCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGAGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGa                               <  1:177494/150‑1 (MQ=255)
                                                                                                                    aaTGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCtt                                                                                    >  2:38886/1‑93 (MQ=255)
                                                                                                                    aaTGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCAc                           >  1:2470013/1‑150 (MQ=255)
                                                                                                                     aTGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCACTTGTGGAAAACCGGCGGTTCCGCGTGGTTCAGGTACTTGGGCTTGTCGTCGCCGAgagat                          >  2:727885/1‑146 (MQ=255)
                                                                                                                         cGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGc                      >  1:470625/1‑150 (MQ=255)
                                                                                                                          gTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTgtc                                       <  1:2369790/132‑1 (MQ=255)
                                                                                                                          gTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTgtc                                       >  2:2369790/1‑132 (MQ=255)
                                                                                                                           tCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       >  1:2404519/1‑99 (MQ=255)
                                                                                                                           tCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCggg                                                                       <  2:2404519/99‑1 (MQ=255)
                                                                                                                              aTGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCg                 >  2:1943127/1‑150 (MQ=255)
                                                                                                                               tGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                     >  1:1992271/1‑97 (MQ=255)
                                                                                                                               tGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTCTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGt                                                                     <  2:1992271/97‑1 (MQ=255)
                                                                                                                                    cccTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                        <  1:1060807/137‑1 (MQ=255)
                                                                                                                                    cccTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCg                        >  2:1060807/1‑137 (MQ=255)
                                                                                                                                     ccTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCga          >  2:1590935/1‑150 (MQ=255)
                                                                                                                                      cTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                        >  1:1239379/1‑87 (MQ=255)
                                                                                                                                      cTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCgg                                                                        <  2:1239379/87‑1 (MQ=255)
                                                                                                                                      cTCGACGAGATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAg                                                      <  1:307596/105‑1 (MQ=255)
                                                                                                                                       tCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg        >  2:338458/1‑150 (MQ=255)
                                                                                                                                       tCGACGAGATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgtg                                                                                 <  1:2041122/77‑1 (MQ=37)
                                                                                                                                       tCGACGAGATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTtgtg                                                                                 >  2:2041122/1‑77 (MQ=37)
                                                                                                                                       tCGACGAGATCATCGTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatg        >  2:18854/1‑150 (MQ=255)
                                                                                                                                         gacgacgaTGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCgatgat      <  1:194633/150‑1 (MQ=255)
                                                                                                                                          acgacgaTGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGatgatg     >  1:423145/1‑150 (MQ=255)
                                                                                                                                          acgacgaTGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGCGAGTTCAGGTACTTCGGATTGTCCTCGCCGATGACCCGGCCACCGAAGGCGatgatg     <  2:507261/150‑1 (MQ=255)
                                                                                                                                           cgacgaTGATCTCGGCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATgc    <  1:2527309/150‑1 (MQ=255)
                                                                                                                                             acgaTGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATgcgc  <  2:423145/150‑1 (MQ=255)
                                                                                                                                             acgaTGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATgcgc  <  1:338458/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
CCGTCAAAACAGAACAGCACGCTGGGCACCACGCGGAACAGGCGCTTGAGGTGCTCTTCGCTGGTGGCGGTGCCAAGGGTGGCCACGGCATTGCGCAGGCCTTGCTGGGCCAGGGCAATGACGTCCATGTAGCCCTCGACGACGATGATCTCGTCGAGGTTGCGGTTGTGCTTGCGCGCCTCGTACAGCCCGTACAGCTCCTGCCCCTTGTGGAAAACCGGGGTTTCCGGGGAGTTCAGGTACTTGGGCTTGTCGTCGCCGAGCACCCGGCCACCGAAGGCGATGATGCGC  >  NC_002947/474053‑474343

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