Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 164,466 C→G 5.9% 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 5.9% 269.0 / 6.4 86G331R (GGC→CGC) pntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (35/45);  new base G (3/2);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

GGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGTT  >  NC_002947/164318‑164613
                                                                                                                                                    |                                                                                                                                                   
ggCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCg                                                                                                                                                    >  1:2219184/1‑150 (MQ=255)
 gCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTACTCTGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGt                                                                                                                                                   >  1:1730483/1‑150 (MQ=255)
  cACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTa                                                                                                                                                  <  2:2161176/150‑1 (MQ=255)
    cTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGc                                                                                                                                                <  2:2362900/150‑1 (MQ=255)
     ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGcc                                                                                                                                                     >  2:1525030/1‑144 (MQ=255)
     ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGcc                                                                                                                                                     <  1:1525030/144‑1 (MQ=255)
     ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGc                                                                                                                                                >  1:565816/1‑149 (MQ=255)
     ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGc                                                                                                                                                <  2:565816/149‑1 (MQ=255)
       cAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGCCCCGCCACCGGCTGCAACGCATCCTTCACGGACCCGGCGGTGTGTGTCTGCTTCTCGGTCTGCTCGTTCATGGCGGGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCg                                                                                                                                             <  2:1464842/150‑1 (MQ=255)
          ccTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGTGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGaa                                                                                                                                          >  2:483676/1‑150 (MQ=255)
          ccTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGaa                                                                                                                                          >  2:340619/1‑150 (MQ=255)
          ccTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGaa                                                                                                                                          <  1:2229052/150‑1 (MQ=255)
           cTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAAc                                                                                                                                         >  1:1404493/1‑150 (MQ=255)
             cGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACga                                                                                                                                       <  2:641226/150‑1 (MQ=255)
                ccAGCAGTACGGTCATGTGCCCGGGCATACGGCCCGCCACCGCGTGCCCCGCATGCTGCTCTGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACgatga                                                                                                                                    <  2:1024200/150‑1 (MQ=255)
                    cagTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCAc                                                                                                                                >  1:1865082/1‑150 (MQ=255)
                      gTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGc                                                                                                                              >  1:2277734/1‑150 (MQ=255)
                      gTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGCCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGc                                                                                                                              <  2:2103585/150‑1 (MQ=255)
                          gTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTc                                                                                                                          >  2:68570/1‑150 (MQ=255)
                           ttCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCg                                                                                                                         <  2:1865082/150‑1 (MQ=255)
                           ttCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCg                                                                                                                         >  1:2092385/1‑150 (MQ=255)
                              aTGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACga                                                                                                                                       >  1:51155/1‑133 (MQ=255)
                              aTGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACga                                                                                                                                       <  2:51155/133‑1 (MQ=255)
                                  gCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGc                                                                                                                  <  1:1912976/150‑1 (MQ=255)
                                   cccGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                                 >  2:957225/1‑150 (MQ=255)
                                   cccGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                                 >  1:210495/1‑150 (MQ=255)
                                     cGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGa                                                                                                               <  1:2106675/150‑1 (MQ=255)
                                     cGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGa                                                                                                               <  2:2277734/150‑1 (MQ=255)
                                     cGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGa                                                                                                               >  2:162840/1‑150 (MQ=255)
                                       ggCCTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCCCGCTGTCGGCGTTGCTGAGc                                                                                                             <  2:1404493/150‑1 (MQ=255)
                                        gCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCt                                                                                                                             <  1:570445/133‑1 (MQ=255)
                                        gCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCt                                                                                                                             >  2:570445/1‑133 (MQ=255)
                                        gCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCa                                                                                                            <  1:957225/150‑1 (MQ=255)
                                          aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGaa                                                                                                              >  2:53757/1‑145 (MQ=255)
                                          aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGaa                                                                                                              <  1:53757/146‑2 (MQ=255)
                                           tACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGa                                                                                                         <  1:1480636/150‑1 (MQ=255)
                                              ggCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGtt                                                                                                                    >  2:194569/1‑136 (MQ=255)
                                              ggCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGtt                                                                                                                    <  1:194569/136‑1 (MQ=255)
                                               gCCCGCCACCGGGTGAGTCGCATGCTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTCGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGt                                                                                                     <  2:799448/150‑1 (MQ=255)
                                                    ccACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACGGCCTGGCGGTAGCCCgc                                                                                                                                            >  2:1078827/1‑105 (MQ=255)
                                                    ccACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACGGCCTGGCGGTAGCCCgc                                                                                                                                            <  1:1078827/106‑2 (MQ=255)
                                                      aCCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCg                                                                                              >  1:2295420/1‑150 (MQ=255)
                                                         gggTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCa                                                                                           <  2:1333664/150‑1 (MQ=255)
                                                          ggTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAg                                                                                          >  1:2327661/1‑150 (MQ=255)
                                                           gTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCCGGCCGTACCCCGGAACGACGACCACGCTGCCGGCGTTCCTGCCCACTAAGGTTGCATCGTCAGc                                                                                         <  2:2365267/150‑1 (MQ=255)
                                                           gTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGc                                                                                         <  2:284726/150‑1 (MQ=255)
                                                                   cATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCgg                                                                                 <  1:1009757/150‑1 (MQ=255)
                                                                    aTACTTCACGGTCACGCCGTTGTGGGTCAGCTTCCCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                <  2:375518/150‑1 (MQ=255)
                                                                        ttCACGGTCACGCCGTTGTGGGTCCGCCTCTCTGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTCGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTc                                                                            >  1:1242660/1‑150 (MQ=255)
                                                                           gcGGTCGCCCCCTTCTGCGTCAGCTGCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCAcc                                                                         <  2:1957740/149‑1 (MQ=255)
                                                                           aCGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCAcc                                                                         <  2:2137604/150‑1 (MQ=255)
                                                                           aCGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCAcc                                                                         <  2:177586/150‑1 (MQ=255)
                                                                              gTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGc                                                                      <  1:483676/150‑1 (MQ=255)
                                                                                   gCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTg                                                                 >  2:289549/1‑150 (MQ=255)
                                                                                          tgGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGc                                                          <  1:939751/150‑1 (MQ=255)
                                                                                            ggTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGa                                                        <  2:1724172/150‑1 (MQ=255)
                                                                                              tCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCAGTAGCCCGGAACGATGATCAGGCTGGCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAAc                                                      >  1:885837/1‑150 (MQ=255)
                                                                                                        ggTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgcgccgc                                            <  1:1652686/150‑1 (MQ=255)
                                                                                                          tCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgcgccgcgc                                          >  2:1148505/1‑150 (MQ=255)
                                                                                                           cAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgcgccgcgcc                                         >  1:15766/1‑150 (MQ=255)
                                                                                                             gCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                <  1:480272/109‑1 (MQ=255)
                                                                                                             gCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                >  2:480272/1‑109 (MQ=255)
                                                                                                              ctctTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                                 >  2:70927/1‑75 (MQ=255)
                                                                                                              ctctTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                                 <  1:70927/75‑1 (MQ=255)
                                                                                                              ctctTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAgt                                      >  2:2142387/1‑150 (MQ=255)
                                                                                                                   tCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCgg                                 <  2:15766/150‑1 (MQ=255)
                                                                                                                    cAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGt                                <  1:2142387/150‑1 (MQ=255)
                                                                                                                     aGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTc                                   >  2:1919166/1‑146 (MQ=255)
                                                                                                                     aGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTc                                   <  1:1919166/146‑1 (MQ=255)
                                                                                                                      gcgcGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAAcc                                                     <  2:1216633/127‑1 (MQ=255)
                                                                                                                      gcgcGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAAcc                                                     >  1:1216633/1‑127 (MQ=255)
                                                                                                                        gcgTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCg                            >  2:2344724/1‑150 (MQ=255)
                                                                                                                                    cgGGCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAgccg                >  2:883241/1‑150 (MQ=255)
                                                                                                                                      ggCTACGGCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAgccgcc              >  1:1933231/1‑150 (MQ=255)
                                                                                                                                         tACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgc                                                 <  2:711968/112‑1 (MQ=255)
                                                                                                                                         tACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgc                                                 >  1:711968/1‑112 (MQ=255)
                                                                                                                                         tACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAg           >  2:487157/1‑150 (MQ=255)
                                                                                                                                          acggCCTGGCGGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAgg          <  1:883241/147‑1 (MQ=255)
                                                                                                                                            cGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGAt        <  1:68570/150‑1 (MQ=255)
                                                                                                                                            cGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTATTGCCAACCGTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGAt        <  1:1239807/150‑1 (MQ=255)
                                                                                                                                                 ggCCGTAGCCCGGACCGCTGATCCCGCTGTCCTCTTTGCCTACCAGGAAGTTTTCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGt   <  2:2092385/150‑1 (MQ=255)
                                                                                                                                                 ggCCGTAGCCCGGAACGATGATCTCGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGt   <  1:1461861/150‑1 (MQ=255)
                                                                                                                                                 ggCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGt   >  1:633113/1‑150 (MQ=255)
                                                                                                                                                  gCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCAt                                                                                                <  1:1624255/56‑1 (MQ=255)
                                                                                                                                                  gCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCAt                                                                                                >  2:1624255/1‑56 (MQ=255)
                                                                                                                                                  gCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                >  1:2023288/1‑72 (MQ=255)
                                                                                                                                                  gCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGa                                                                                <  2:2023288/72‑1 (MQ=255)
                                                                                                                                                  gCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTAGAATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGtt  >  1:169782/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGATGACGTT  >  NC_002947/164318‑164613

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