Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,764,437 T→C 7.3% D203G (GAC→GGC)  PP_2418 ← cobalamin ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,764,4370TC 7.3% 302.3 / 8.8 105D203G (GAC→GGC) PP_2418cobalamin ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (45/52);  new base C (6/2);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

Only 100 of 110 total aligned reads displayed.

CGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAGCTCATCGAACTG  >  NC_002947/2764288‑2764585
                                                                                                                                                     |                                                                                                                                                    
cggccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                                      >  1:470656/1‑150 (MQ=255)
cggccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                                      >  2:527706/1‑150 (MQ=255)
cggccACACAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                                      >  2:2377794/1‑150 (MQ=255)
 ggccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTc                                                                                                                                                     <  1:840081/150‑1 (MQ=255)
   ccACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGa                                                                                                                                                   >  2:1623262/1‑150 (MQ=255)
      caTAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGaatt                                                                                                                                                <  1:1963619/150‑3 (MQ=37)
           gggTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCg                                                                                                                                           <  1:1011948/150‑1 (MQ=255)
           gggTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCg                                                                                                                                           <  2:633787/150‑1 (MQ=255)
             gTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGc                                                                                                                                         <  1:2196754/150‑1 (MQ=255)
             gTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGc                                                                                                                                         >  2:715583/1‑150 (MQ=255)
              tCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGcc                                                                                                                                        >  2:10163/1‑150 (MQ=255)
                 aTGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGt                                                                                                                                     <  2:257183/150‑1 (MQ=255)
                  tGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGtt                                                                                                                                    <  2:765931/150‑1 (MQ=255)
                        ggTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGa                                                                                                                              <  1:249357/150‑1 (MQ=255)
                              cgcTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGctct                                                                                                                        <  2:2212348/150‑1 (MQ=255)
                               gcTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCtt                                                                                                                       >  2:307692/1‑150 (MQ=255)
                               gcTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCtt                                                                                                                       <  1:227371/150‑1 (MQ=255)
                                      aaCCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCa                                                                                                                                       <  1:1921628/127‑1 (MQ=255)
                                         cAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTg                                                                                                             <  1:1784667/150‑1 (MQ=255)
                                         cAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTg                                                                                                             >  1:276908/1‑150 (MQ=255)
                                         cAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGaa                                                                                                                                                  <  2:1859374/113‑1 (MQ=25)
                                         cAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGaa                                                                                                                                                  >  1:1859374/1‑113 (MQ=25)
                                          aGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGt                                                                                                            >  1:396555/1‑150 (MQ=255)
                                            aaCAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCTCGCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCt                                                                                                          <  2:1796036/150‑1 (MQ=255)
                                            aaCAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCt                                                                                                          <  2:1795854/150‑1 (MQ=255)
                                              cAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCTGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTAc                                                                                                        <  1:1623262/150‑1 (MQ=255)
                                                 cgcCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGAGTTCAGTGCCGCGCCGCGCTCCTGAACGCCCAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTg                                                                                                     >  1:2379543/1‑150 (MQ=255)
                                                  gcCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTg                                                                                                     >  2:1571805/1‑149 (MQ=255)
                                                  gcCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTg                                                                                                     <  1:1571805/149‑1 (MQ=255)
                                                   ccGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTg                                                                                                                      <  2:1910554/131‑1 (MQ=255)
                                                   ccGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTg                                                                                                                      >  1:1910554/1‑131 (MQ=255)
                                                         ttGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCtt                                                                                             >  2:2139283/1‑150 (MQ=255)
                                                          tGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTg                                                                                            <  2:7234/150‑1 (MQ=255)
                                                          tGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTg                                                                                            <  1:10163/150‑1 (MQ=255)
                                                          tGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCCGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTg                                                                                            <  1:1535874/150‑1 (MQ=255)
                                                           gCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGc                                                                                           >  2:289059/1‑150 (MQ=255)
                                                             cGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAg                                                                                         <  2:470656/150‑1 (MQ=255)
                                                                cAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCg                                                                                      >  1:2038927/1‑150 (MQ=255)
                                                                 aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     <  1:1956082/150‑1 (MQ=255)
                                                                 aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     <  1:2139283/150‑1 (MQ=255)
                                                                 aGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     <  2:412935/150‑1 (MQ=255)
                                                                   gTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                                      >  2:1690658/1‑83 (MQ=255)
                                                                   gTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGt                                                                                                                                                      <  1:1690658/83‑1 (MQ=255)
                                                                   gTCCTTGCCGGCCAGCTGTGCCCTGGCGCAGTCGAAACATCGCTGGCACTCTGCGGTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCTCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGt                                                                                   <  1:2377794/150‑1 (MQ=255)
                                                                      cTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGGTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGcc                                                                                >  2:336143/1‑150 (MQ=255)
                                                                       ttGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCg                                                                               >  2:917781/1‑150 (MQ=255)
                                                                            ggCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGc                                                                          >  1:1496412/1‑150 (MQ=255)
                                                                            ggCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCGCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTACCGTCTGc                                                                          >  2:1309600/1‑150 (MQ=255)
                                                                             gCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCtt                                                                                                                       <  2:1943249/104‑1 (MQ=255)
                                                                             gCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCtt                                                                                                                       >  1:1943249/1‑104 (MQ=255)
                                                                              ccAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGACGGCACCCGCCTGCaa                                                                        <  1:606015/150‑2 (MQ=255)
                                                                               cAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCAcc                                                                       <  2:1337643/150‑1 (MQ=255)
                                                                                  cTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGAc                                                                    <  2:808015/150‑1 (MQ=255)
                                                                                         cTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGTTCATGGACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTg                                                             >  1:1356877/1‑150 (MQ=255)
                                                                                         cTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     <  2:1020945/126‑1 (MQ=255)
                                                                                         cTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     >  1:1020945/1‑126 (MQ=255)
                                                                                          tGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGt                                                            >  1:1762863/1‑150 (MQ=255)
                                                                                            gcgcTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTcc                                                          <  2:2000480/150‑1 (MQ=255)
                                                                                             cgcTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGATGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCt                                                         <  2:1378142/150‑1 (MQ=255)
                                                                                               cTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTg                                                       <  2:396555/150‑1 (MQ=255)
                                                                                                tGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAg                                                                 >  1:1508332/1‑139 (MQ=255)
                                                                                                         gcCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     >  2:2385770/1‑110 (MQ=255)
                                                                                                         gcCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCgg                                                                                     <  1:2385770/110‑1 (MQ=255)
                                                                                                           ccTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTAc                                                                                                        >  2:160172/1‑89 (MQ=255)
                                                                                                           ccTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTAc                                                                                                        <  1:160172/89‑1 (MQ=255)
                                                                                                           ccTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTcc                                           >  2:922606/1‑150 (MQ=255)
                                                                                                           ccTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGGCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTcc                                           >  1:247666/1‑150 (MQ=255)
                                                                                                              tGCAATTCGGCTTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGAt                                        <  2:1356877/150‑1 (MQ=255)
                                                                                                                cAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATgg                                      <  2:2038927/150‑1 (MQ=255)
                                                                                                                atcttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATgg                                      >  1:1469132/4‑150 (MQ=37)
                                                                                                                 aaTTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGgc                                     >  1:1503912/1‑150 (MQ=255)
                                                                                                                  aTTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGgca                                    <  1:546368/150‑1 (MQ=255)
                                                                                                                   ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTg                                                                           <  2:2237642/110‑1 (MQ=37)
                                                                                                                   ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTg                                                                           >  1:2237642/1‑110 (MQ=37)
                                                                                                                   ttCGACGTTCAGGACCGCGGCGCGGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGgcag                                   >  2:1780742/1‑150 (MQ=38)
                                                                                                                       gCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACg                               >  2:930974/1‑150 (MQ=255)
                                                                                                                         gTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTa                             >  1:2145338/1‑150 (MQ=255)
                                                                                                                          ttcaggaccgcggcgcgGCACTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAg                            >  2:352967/8‑150 (MQ=255)
                                                                                                                          ttcaggaccgcggcgcgGCACTGAACGCCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAg                            >  2:868921/8‑150 (MQ=255)
                                                                                                                          ttCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAg                            >  2:2188104/1‑150 (MQ=255)
                                                                                                                          ttCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCACGATTGCCTTGCCCTCGCAGTCCGATGGCAGCACGTAg                            <  1:2410877/150‑1 (MQ=255)
                                                                                                                            cAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTTGCCCCTTACCAAGCGGAGGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCGCGTAGGt                          >  1:1476424/1‑150 (MQ=255)
                                                                                                                             aGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTg                         <  1:922606/150‑1 (MQ=255)
                                                                                                                             aGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTg                         <  1:336143/150‑1 (MQ=255)
                                                                                                                               tGCCGCGCCGCGGTCGTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTAGCAGTCCGATGGCAGCACGTAGGTggg                       >  2:867826/1‑150 (MQ=255)
                                                                                                                                 ccgcgccgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGa                     <  1:715583/150‑1 (MQ=255)
                                                                                                                                  cgcgccgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGAt                    <  1:191862/150‑1 (MQ=255)
                                                                                                                                   gcgccgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCCCGTAGGTGGGGATc                   >  2:1149399/1‑150 (MQ=255)
                                                                                                                                       cgcgGTCCTGAACGTCGAAGATCTCGGCCAGTTAGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCGGCACGTAGGTGGGGATCTTcg               >  2:1179498/1‑149 (MQ=255)
                                                                                                                                        gcgGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTTTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAg              <  1:917781/150‑1 (MQ=255)
                                                                                                                                              cTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCtt                                                        >  2:1066859/1‑102 (MQ=255)
                                                                                                                                              cTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCtt                                                        <  1:1066859/102‑1 (MQ=255)
                                                                                                                                               tGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAg                                                                 >  2:1926276/1‑92 (MQ=255)
                                                                                                                                               tGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAg                                                                 <  1:1926276/92‑1 (MQ=255)
                                                                                                                                                gAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAGCTCATCGa      >  1:1473657/1‑150 (MQ=255)
                                                                                                                                                 aaCGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAGCTCATCGaa     <  1:527706/150‑1 (MQ=255)
                                                                                                                                                   cGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAGCTCATCGAACt   >  2:1656474/1‑150 (MQ=255)
                                                                                                                                                    gTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGc                                                      <  1:418504/98‑1 (MQ=255)
                                                                                                                                                    gTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGc                                                      >  2:418504/1‑98 (MQ=255)
                                                                                                                                                    gTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGAAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAGCTCATCGAACTg  >  1:1772519/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGCCACATAGGGGTCGATGTCCAGGTCGGCGCTGGAGAACCAGAACAGCGCCGAAGTTGCCGACAGGTCCTTGCCGGCCAGCTGTGCCCTGGCGCTGTCGAGGCGCCCTTGCAATTCGGCGTTCAGTGCCGCGCCGCGGTCCTGAACGTCGAAGATCTCGGCCAGTTCGCTGACGCTCTTGTAGATGCTGTCTACCTGGAACGCTTGCAGGCGGGTGCCGTCTGCACCGACCAGGTTGTCCTTGCCCTCGCAGTCCGATGGCAGCACGTAGGTGGGGATCTTCAGCTCATCGAACTG  >  NC_002947/2764288‑2764585

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: