Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 164,462 A→T 6.2% L332Q (CTG→CAG)  pntB ← pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947164,4620AT 6.2% 249.9 / 5.1 81L332Q (CTG→CAG) pntBpyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (33/42);  new base T (3/2);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGAT  >  NC_002947/164314‑164607
                                                                                                                                                    |                                                                                                                                                 
gTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAg                                                                                                                                                  <  1:1327194/150‑1 (MQ=255)
    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                                                                                                                                          <  2:565816/149‑1 (MQ=255)
         ttCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGc                                                                                                                                          >  1:565816/1‑149 (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                                                                                                                   >  1:2092385/1‑150 (MQ=255)
                               ttCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCg                                                                                                                   <  2:1865082/150‑1 (MQ=255)
                                  aTGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACga                                                                                                                                 <  2:51155/133‑1 (MQ=255)
                                  aTGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACga                                                                                                                                 >  1:51155/1‑133 (MQ=255)
                                      gCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGc                                                                                                            <  1:1912976/150‑1 (MQ=255)
                                       cccGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                           >  2:957225/1‑150 (MQ=255)
                                       cccGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                           >  1:210495/1‑150 (MQ=255)
                                         cGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGa                                                                                                         <  2:2277734/150‑1 (MQ=255)
                                         cGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGa                                                                                                         <  1:2106675/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                                                                                                        <  1:53757/146‑2 (MQ=255)
                                              aTACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGaa                                                                                                        >  2:53757/1‑145 (MQ=255)
                                               tACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGa                                                                                                   <  1:1480636/150‑1 (MQ=255)
                                                  ggCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGtt                                                                                                              <  1:194569/136‑1 (MQ=255)
                                                  ggCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGtt                                                                                                              >  2:194569/1‑136 (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:177586/150‑1 (MQ=255)
                                                                               aCGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCAcc                                                                   <  2:2137604/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                                                                                                           <  1:70927/75‑1 (MQ=255)
                                                                                                                  ctctTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCt                                                                                                           >  2:70927/1‑75 (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                                           >  1:711968/1‑112 (MQ=255)
                                                                                                                                             tACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTgc                                           <  2:711968/112‑1 (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)
                                                                                                                                                    |                                                                                                                                                 
GTACGGCACTTCAGCCTCGGCCAGCAGTACGTTCATGTGCCCGGGCATACGGCCCGCCACCGGGTGAATCGCATACTTCACGGTCACGCCGTTGTGGGTCAGCTTCTCGGTCAGCTCTTTCAGCGCGTGTTGTGCGCGGGCTACCGCCAGGCCGTAGCCCGGAACGATGATCACGCTGTCGGCGTTGCTGAGCAGGAAGGTGGCATCGTCAGCCGAACCGGACTTCACCGGCCGCTGTTCTTGCGAACCTTGCGCCGCGCCAGTGTCGGTATCGCCACCGAAGCCGCCAAGGAT  >  NC_002947/164314‑164607

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