Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,904,907 C→T E231K (GAA→AAA)  bglF ← beta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,904,9070CT98.4% 208.3 / NA 63E231K (GAA→AAA) bglFbeta‑glucoside specific PTS enzyme II/BglG kinase/BglG phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (27/35);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGG  >  NC_000913/3904670‑3905149
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
ccAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGcac                                                                                                                                                                                                                                       <  2:171586/251‑1 (MQ=255)
ccAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGcac                                                                                                                                                                                                                                       <  1:522986/251‑1 (MQ=255)
   ggCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTTATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAAc                                                                                                                                                                                                                                    <  1:414550/251‑1 (MQ=255)
    gCCCCAGTGCCGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGATTCAGAATGCTGCACAAcc                                                                                                                                                                                                                                   >  2:84686/1‑251 (MQ=255)
     ccccAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCa                                                                                                                                                                                                                                  <  2:476778/251‑1 (MQ=255)
      cccAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAAcc                                                                                                                                                                                                                                   >  2:440844/1‑249 (MQ=255)
          gTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAg                                                                                                                                                                                                                             <  2:501787/251‑1 (MQ=255)
           tGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGc                                                                                                                                                                                                                               >  2:320687/1‑248 (MQ=255)
                  cGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAaata                                                                                                                                                                                                                     <  2:388614/251‑1 (MQ=255)
                   gAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAaataa                                                                                                                                                                                                                    <  1:84686/251‑1 (MQ=255)
                   gAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAaata                                                                                                                                                                                                                     >  2:195315/1‑250 (MQ=255)
                   gAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAAACGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGATTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTTAGCATGCTGCACAACCAGGCAGaaaaacta                                                                                                                                                                                                                     >  1:354912/1‑247 (MQ=255)
                     aCATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATg                                                                                                                                                                                                                  <  2:108410/251‑1 (MQ=255)
                      cATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATgg                                                                                                                                                                                                                 <  1:440844/251‑1 (MQ=255)
                       aTGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATggg                                                                                                                                                                                                                >  2:103882/1‑251 (MQ=255)
                                   tGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg                                                                                                                                                                                                     <  2:248343/250‑1 (MQ=255)
                                    gCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTTCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGAc                                                                                                                                                                                                   <  1:115751/251‑1 (MQ=255)
                                     ccAGAAGCCGCCCATTACCGCGCCCGCCAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACg                                                                                                                                                                                                  <  2:108809/251‑1 (MQ=255)
                                      cAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACg                                                                                                                                                                                                  >  1:13482/1‑250 (MQ=255)
                                       aGAAGCCGCCCATTTCCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTTATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAg                                                                                                                                                                                                <  1:487951/251‑1 (MQ=255)
                                        gAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGt                                                                                                                                                                                               <  1:151937/251‑1 (MQ=255)
                                        gAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGt                                                                                                                                                                                               <  1:103882/251‑1 (MQ=255)
                                        gAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGACAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAAGGCTGCACAACCAGGCAGAAAAAATAATGGGAATATCCGATGACGAGt                                                                                                                                                                                               <  1:451369/251‑1 (MQ=255)
                                            ccgccCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATGAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAAtt                                                                                                                                                                                           <  1:270231/251‑1 (MQ=255)
                                                           cccGCAAATGCAGGAACCGTCTGATAAGGCCAGCGATAACCGGCGGCAATTAGTTCGCGTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCAGCAGACATAGCAATGGTGTGAAGAAATTTTTGATGGCAGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAAGAACCGATGACGATTAATTCAACAATGTGACCgg                                                                                                                                                                            <  2:435397/251‑1 (MQ=255)
                                                           cccGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGGTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCgg                                                                                                                                                                            <  2:282261/251‑1 (MQ=255)
                                                                       ggAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATCCCCAGGaa                                                                                                                                                                <  2:325342/251‑1 (MQ=255)
                                                                           ccGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg                                                                                                                                                                                                     >  1:119709/1‑210 (MQ=255)
                                                                           ccGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAAGTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATg                                                                                                                                                                                                     <  2:119709/210‑1 (MQ=255)
                                                                                  aTAAAGCCAGAGATAACCGGCGCCAATCAGTTCGCTTATCCAGGTTGCTAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGATTGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGACATCCAGccccc                                                                                                                                                     >  2:122246/1‑250 (MQ=255)
                                                                                        ccAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAACGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTACGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCAt                                                                                                                                               >  2:158288/1‑251 (MQ=255)
                                                                                                       ggcAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGt                                                                                                                                <  2:432077/251‑1 (MQ=255)
                                                                                                             cAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaaca                                                                                                                                                                                      <  2:320471/191‑1 (MQ=255)
                                                                                                             cAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTcaaca                                                                                                                                                                                      >  1:320471/1‑191 (MQ=255)
                                                                                                              aGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGa                                                                                                                          >  1:577108/1‑250 (MQ=255)
                                                                                                               gTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGGTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGaaa                                                                                                                        <  1:364846/251‑1 (MQ=255)
                                                                                                                   gCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGACTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGc                                                                                                                      >  2:396695/1‑249 (MQ=255)
                                                                                                                       aTCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAAACATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTc                                                                                                                  >  1:491056/1‑249 (MQ=255)
                                                                                                                          cAGGTTGGTAGCGGCCCCAACAGCAGAAAGGTTACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTACTTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCACGTTCTCGACAGCAGTCAGAAt                                                                                                             >  2:162577/1‑251 (MQ=255)
                                                                                                                            ggTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAAtt                                                                                                            >  2:124936/1‑250 (MQ=255)
                                                                                                                                   aGCGGCCCCACCAGCAGAAAGGTGACTGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCTCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATg                                                                                                    <  2:454016/251‑1 (MQ=255)
                                                                                                                                         cccACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCACTGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTACTGGATGCACTa                                                                                               >  1:95351/1‑250 (MQ=255)
                                                                                                                                          ccaccaGCAGAAAGGTGCCGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGATGCACAACCAGGCAGAAAAAATACTGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGCTGCACTAAc                                                                                             >  2:205213/1‑251 (MQ=255)
                                                                                                                                          ccaccaGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAAc                                                                                             <  1:145898/251‑1 (MQ=255)
                                                                                                                                            accaGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGc                                                                                           >  1:384699/1‑251 (MQ=255)
                                                                                                                                               agcagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTcc                                                                                        <  2:103983/251‑1 (MQ=255)
                                                                                                                                                gcagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTcca                                                                                       <  1:320687/251‑1 (MQ=255)
                                                                                                                                                 cagAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTccac                                                                                      >  1:169737/1‑251 (MQ=255)
                                                                                                                                                   gAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTccacc                                                                                     >  1:394431/1‑250 (MQ=255)
                                                                                                                                                       ggTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATc                                                                                 >  1:318796/1‑250 (MQ=255)
                                                                                                                                                          gACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCAcc                                                                              <  2:579615/250‑1 (MQ=255)
                                                                                                                                                             gggTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAACTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCAGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATgg                                                                          >  2:2528/1‑251 (MQ=255)
                                                                                                                                                                   gATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTa                                                                     >  1:445039/1‑250 (MQ=255)
                                                                                                                                                                             aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCCCAACCAGGCAGCAAAAATAATGGGAATAACCGATGACGAGTACTTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCCGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAcc                                                          >  1:586303/1‑251 (MQ=255)
                                                                                                                                                                             aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAcc                                                          <  1:541569/251‑1 (MQ=255)
                                                                                                                                                                                     cAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCCCGAAAGCAGTCAGAATTAATGGAGGCGCTAACGCTCCACCAATCACCATGGCAGTAAAGGGATGACCGCCGAAGc                                                  <  2:577108/251‑1 (MQ=255)
                                                                                                                                                                                                        ttGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATcc                               >  1:339685/1‑251 (MQ=255)
                                                                                                                                                                                                               ccGACGGTAACCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATTCCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGGCAGAATTAATTGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGAtt                        <  1:489720/251‑1 (MQ=255)
                                                                                                                                                                                                                        aCCACGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAAGAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAAGGGATGCACTAACGCCCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTGCCCCGCGGTGTATCCCAGGATTATCGGgaag               <  2:229704/251‑1 (MQ=255)
                                                                                                                                                                                                                           aCGCATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAAc            <  2:394431/251‑1 (MQ=255)
                                                                                                                                                                                                                              cATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAACTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCa          >  1:264527/1‑250 (MQ=255)
                                                                                                                                                                                                                              cATTAAGTCGGCGTTTCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCaa         <  1:195315/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cGGCGTTCCAAAAGGCTGCACAACCAGGCACAAAAAATAATGGGAATAACCGATGACGAGTAATGCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGTGCATCCGCTTTTTGCCCGTTCGCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAACCGCCATGGCAGTAAATGGATTACCGCCGAAGCGTTGCCCCGCGGGGTATCCCAGGATTATCGGGAAGAACCAAAACAAgg  <  2:51514/250‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
CCAGGCCCCAGTGCAGTCCGAACATGACGAAGATTTGCCAGAAGCCGCCCATTACCGCGCCCGCAAATGCAGGAACCGCCTGATAAAGCCAGAGATAACCGGCGGCAATCAGTTCGCTTATCCAGGTTGATAGCGGCCCCACCAGCAGAAAGGTGACGGGTGTGATAACCATCAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAAGAACCAAAACAAGG  >  NC_000913/3904670‑3905149

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: