Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 164,466 C→G 8.7% G331R (GGC→CGC)  pntB ← pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947164,4660CG 8.7% 330.6 / 16.0 108G331R (GGC→CGC) pntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (53/45);  new base G (7/3);  total (60/48)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGTTGA  >  NC_002947/164323‑164615
                                                                                                                                               |                                                                                                                                                     
ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGACCGTAGcc                                                                                                                                                 >  2:2067664/1‑150 (MQ=255)
ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGAGGGCTACGGCCTGGCGGTAGcc                                                                                                                                                 >  1:125381/1‑150 (MQ=255)
 tCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACGGCCTGGCGGTAGccc                                                                                                                                                <  2:1927373/150‑1 (MQ=255)
 tCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACGGCCTGGCGGTAGccc                                                                                                                                                <  2:125381/150‑1 (MQ=255)
 tCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGccc                                                                                                                                                <  2:1491709/150‑1 (MQ=255)
     ccTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGaa                                                                                                                                            <  1:1890602/150‑1 (MQ=255)
     ccTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGaa                                                                                                                                            >  2:1890602/1‑150 (MQ=255)
               cagTACGTTCATGTTCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCATGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCAc                                                                                                                                  <  2:174423/150‑1 (MQ=255)
               cagTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCAc                                                                                                                                  <  2:2270302/150‑1 (MQ=255)
                agTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACg                                                                                                                                 >  2:1272600/1‑150 (MQ=255)
                 gTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGc                                                                                                                                >  2:728062/1‑150 (MQ=255)
                     gTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTc                                                                                                                            >  2:2196522/1‑150 (MQ=255)
                     gTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTGGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTCCCGCCAGTCCGTAGCCCGGAACGATGATCACGCTGTc                                                                                                                            >  2:541399/1‑150 (MQ=255)
                         aTGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCg                                                                                                                        >  2:1682088/1‑150 (MQ=255)
                           gtgCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGtt                                                                                                                      >  1:1724092/1‑150 (MQ=255)
                           gtgCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGGGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGtt                                                                                                                      <  1:429401/150‑1 (MQ=255)
                              cccGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                                   <  1:348651/150‑1 (MQ=255)
                                cGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGa                                                                                                                 <  2:32183/150‑1 (MQ=255)
                                   gCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCATTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCt                                                                                                                               <  1:283433/133‑1 (MQ=255)
                                   gCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCATTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCt                                                                                                                               >  2:283433/1‑133 (MQ=255)
                                     aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCAGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAgg                                                                                                            >  1:2283741/1‑150 (MQ=255)
                                        cGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAg                                                                                                         <  1:470584/150‑1 (MQ=255)
                                         ggCCCGCCACCGGGTGAATCGCATACTTCACGGCCACTCCGATGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTATTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAgg                                                                                                        >  1:1398664/1‑150 (MQ=255)
                                           cccGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTg                                                                                                      <  2:1995626/150‑1 (MQ=255)
                                           cccGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTg                                                                                                      <  1:1272600/150‑1 (MQ=255)
                                           cccGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCAGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTg                                                                                                      <  2:2283741/150‑1 (MQ=255)
                                           cccGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGGGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGCTGACCACGCTGAAAAAACTGCTGAGCAGGAAGGTg                                                                                                      >  1:1329974/1‑150 (MQ=255)
                                             cGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGc                                                                                                    >  1:2062826/1‑150 (MQ=255)
                                              gCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGc                                                                                                    >  1:357285/1‑149 (MQ=255)
                                              gCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGc                                                                                                    <  2:357285/149‑1 (MQ=255)
                                                catcgtgTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATc                                                                                                 <  2:1167981/144‑1 (MQ=255)
                                                cACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATc                                                                                                 <  2:1075118/150‑1 (MQ=255)
                                                 aCCGGGTGAATCGCATACTTCACGGTCGACCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCg                                                                                                >  1:1642707/1‑150 (MQ=255)
                                                     ggTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAg                                                                                            <  2:2252973/150‑1 (MQ=255)
                                                      gTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTc                                                                                                                            >  2:2345970/1‑117 (MQ=255)
                                                      gTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTc                                                                                                                            <  1:2345970/117‑1 (MQ=255)
                                                          aTCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGaa                                                                                       >  1:807786/1‑150 (MQ=255)
                                                           tCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAAc                                                                                      >  1:256271/1‑150 (MQ=255)
                                                           tCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTATTGCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAAc                                                                                      <  2:1642707/150‑1 (MQ=255)
                                                             gCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGc                                                                                                                         >  2:574129/1‑113 (MQ=255)
                                                             gCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGc                                                                                                                         <  1:574129/113‑1 (MQ=255)
                                                              cATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTTCCTGGAGGCGGTAGTCCGGAACGATGATCACGCTGTCGGCGTTGCGGAGCAGGAAGGTGTGATCGTCAGCCGAACCgg                                                                                   >  1:899811/1‑150 (MQ=255)
                                                               aTACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                  <  1:2092916/150‑1 (MQ=255)
                                                               aTACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                  <  1:1682088/150‑1 (MQ=255)
                                                                tACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGAc                                                                                 <  1:1781874/150‑1 (MQ=255)
                                                                   ttCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTc                                                                              >  1:440643/1‑150 (MQ=255)
                                                                   ttCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTc                                                                              >  1:1503806/1‑150 (MQ=255)
                                                                        ggTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTg                                                                                                      <  1:1547080/121‑1 (MQ=255)
                                                                         gTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGc                                                                        >  2:860887/1‑150 (MQ=255)
                                                                          tCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGcc                                                                       >  1:857777/1‑150 (MQ=255)
                                                                          tCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGcc                                                                       >  1:1690917/1‑150 (MQ=255)
                                                                           cACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCg                                                                      <  2:1690917/150‑1 (MQ=255)
                                                                              gCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTg                                                                   >  1:818139/1‑150 (MQ=255)
                                                                               ccGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGt                                                                  >  2:1535921/1‑150 (MQ=255)
                                                                                   tgtgGGTCAGTTTCTCGGTCACCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCtt                                                              <  1:125015/150‑1 (MQ=255)
                                                                                   tgtgGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCtt                                                              <  1:175533/150‑1 (MQ=255)
                                                                                    gtgGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGCCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTg                                                             <  1:937958/150‑1 (MQ=255)
                                                                                       ggTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGa                                                          >  2:1985795/1‑150 (MQ=255)
                                                                                        gTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACgatga                                                                                                                                      <  1:1120663/73‑1 (MQ=255)
                                                                                        gTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACgatga                                                                                                                                      >  2:1120663/1‑73 (MQ=255)
                                                                                             cTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCa                                                                                                              >  2:142458/1‑92 (MQ=255)
                                                                                             cTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCa                                                                                                              <  1:142458/92‑1 (MQ=255)
                                                                                              ttCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAg                                                                                                         >  1:1609334/1‑96 (MQ=255)
                                                                                              ttCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgc                                                   >  2:180971/1‑150 (MQ=255)
                                                                                              ttCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTCTCGGAGTTGCTGTGCAGGAAg                                                                                                         <  2:1609334/96‑1 (MQ=255)
                                                                                                ctcGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACgat                                                                                                                                        >  2:1525101/1‑63 (MQ=255)
                                                                                                   ggTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgcgccgc                                              >  2:933951/1‑150 (MQ=255)
                                                                                                   ggTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgcgccgc                                              >  2:1539545/1‑150 (MQ=255)
                                                                                                     tCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTTCTGAGCAGGAAGGTGCCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgcgccgcgc                                            >  1:1145082/1‑150 (MQ=255)
                                                                                                       aGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTg                                                                                                                     <  2:497084/75‑1 (MQ=255)
                                                                                                       aGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTg                                                                                                                     >  1:497084/1‑75 (MQ=255)
                                                                                                       aGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAg                                                                                                         >  1:2220692/1‑87 (MQ=255)
                                                                                                       aGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAg                                                                                                         <  2:2220692/87‑1 (MQ=255)
                                                                                                       aGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCa                                          >  1:2311517/1‑150 (MQ=255)
                                                                                                         ctctTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAgt                                        <  1:2002094/150‑1 (MQ=255)
                                                                                                         ctctTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAgt                                        >  1:496023/1‑150 (MQ=255)
                                                                                                            tttCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTc                                     >  2:1474241/1‑150 (MQ=255)
                                                                                                             ttCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGAGGCAGGCCTGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCCGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCg                                    <  2:1398664/150‑1 (MQ=255)
                                                                                                              tCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCgg                                   >  1:1572869/1‑150 (MQ=255)
                                                                                                                 gcgcGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGc                                                                                                                    >  2:1993220/1‑66 (MQ=255)
                                                                                                                   gcgTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCg                              >  2:1065599/1‑150 (MQ=255)
                                                                                                                   gcgTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCg                              >  1:361969/1‑150 (MQ=255)
                                                                                                                     gtgtTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGGGTTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGcc                            >  1:1331652/1‑150 (MQ=255)
                                                                                                                       gttgtGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGCGCAGGAACGTGCCATCTTCAGCCTAACCTGACGTCACCGTCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCAc                          <  2:2079177/150‑1 (MQ=255)
                                                                                                                        atcgttccgGGCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCAcc                         >  1:1074600/8‑150 (MQ=255)
                                                                                                                           tgCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGaa                      <  1:933951/150‑1 (MQ=255)
                                                                                                                           tgCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGaa                      <  2:818139/150‑1 (MQ=255)
                                                                                                                            ttccgGGCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAg                     >  1:924092/4‑150 (MQ=255)
                                                                                                                                gggCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAgccgc                 >  1:469763/1‑150 (MQ=255)
                                                                                                                                gggCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAgccgc                 >  2:1290267/1‑150 (MQ=255)
                                                                                                                                gggCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCa                           >  1:31910/1‑140 (MQ=255)
                                                                                                                                gggCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCa                           <  2:31910/140‑1 (MQ=255)
                                                                                                                                  gCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCa               >  2:2121216/1‑150 (MQ=255)
                                                                                                                                    tACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAg             >  1:2080538/1‑150 (MQ=255)
                                                                                                                                    tACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGAGACCTTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAg             <  1:633866/150‑1 (MQ=255)
                                                                                                                                     aCCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAgg            <  2:361969/150‑1 (MQ=255)
                                                                                                                                      ccGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCTAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGa           <  1:1135192/150‑1 (MQ=255)
                                                                                                                                          cAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGAc       <  1:1539545/150‑1 (MQ=255)
                                                                                                                                           aGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACg      >  2:2195442/1‑150 (MQ=255)
                                                                                                                                               cGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGTTga  <  2:1329974/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGTTGA  >  NC_002947/164323‑164615

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