Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 164,458 C→G 7.5% A333A (GCG→GCC pntB ← pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947164,4580CG 7.5% 337.7 / 14.8 107A333A (GCG→GCCpntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (53/46);  new base G (6/2);  total (59/48)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 110 total aligned reads displayed.

GTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGA  >  NC_002947/164314‑164606
                                                                                                                                                |                                                                                                                                                    
gTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAg                                                                                                                                                 <  1:81121/150‑1 (MQ=255)
gTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTCCAGCTCGTGTTGTGCGCGGGCTACCGCCAg                                                                                                                                                 <  1:1270676/150‑1 (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)
                  ggCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACgat                                                                                                                               <  2:311510/150‑1 (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)
                                    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)
                                              aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTAcc                                                                                                                                                      >  1:1595293/1‑99 (MQ=255)
                                              aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTAcc                                                                                                                                                      <  2:1595293/99‑1 (MQ=255)
                                              aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCAGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAgg                                                                                                   >  1:2283741/1‑150 (MQ=255)
                                                 cGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAg                                                                                                <  1:470584/150‑1 (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                                                                                           <  2:357285/149‑1 (MQ=255)
                                                       gCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGc                                                                                           >  1:357285/1‑149 (MQ=255)
                                                         catcgtgTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATc                                                                                        <  2:1167981/144‑1 (MQ=255)
                                                         cACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATc                                                                                        <  1:728062/150‑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)
                                                               gTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGACCACGCTGTCGGCGTTGCTTAGCAGGAAGGTGGCATCGTCaac                                                                                  >  2:125015/1‑148 (MQ=255)
                                                                   aTCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGaa                                                                              >  1:807786/1‑150 (MQ=255)
                                                                    tCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAAc                                                                             >  1:256271/1‑150 (MQ=255)
                                                                      gCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGc                                                                                                                <  1:574129/113‑1 (MQ=255)
                                                                      gCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGc                                                                                                                >  2:574129/1‑113 (MQ=255)
                                                                        aTACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                         <  1:1682088/150‑1 (MQ=255)
                                                                        aTACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                         <  1:2092916/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)
                                                                                 ggTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTg                                                                                             >  2:1547080/1‑121 (MQ=255)
                                                                                  gTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGc                                                               >  2:860887/1‑150 (MQ=255)
                                                                                   tCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGcc                                                              >  1:1690917/1‑150 (MQ=255)
                                                                                   tCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGcc                                                              >  1:857777/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)
                                                                                              tgGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGc                                                   >  2:2002094/1‑150 (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                                                                                                                               <  1:1525101/63‑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)
                                                                                                                  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)
                                                                                                                          gcgcGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGc                                                                                                           <  1:1993220/66‑1 (MQ=255)
                                                                                                                            gcgTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCg                     >  1:361969/1‑150 (MQ=255)
                                                                                                                            gcgTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCg                     >  2:1065599/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             <  2:818139/150‑1 (MQ=255)
                                                                                                                                    tgCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGaa             <  1:933951/150‑1 (MQ=255)
                                                                                                                                     ttccgGGCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAg            >  1:924092/4‑150 (MQ=255)
                                                                                                                                         gggCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCa                  >  1:31910/1‑140 (MQ=255)
                                                                                                                                          ggCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAgccgcc       >  1:320858/1‑150 (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)
                                                                                                                                                |                                                                                                                                                    
GTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGA  >  NC_002947/164314‑164606

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: