Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,708 C→T intergenic (‑475/‑2199) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,233,7080CT95.8% 630.4 / 10.0 191intergenic (‑475/‑2199)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (106/77);  total (110/81)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 192 total aligned reads displayed.

GGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCG  >  NC_000913_3_bme_pgi/4233561‑4233854
                                                                                                                                                   |                                                                                                                                                  
   gCGGTTGTCAACGATGGGGGCATGCGGATTTTTCCTCCACACCTGGCGGTCTGTAGTTCAGATAATAAATGCTTCACTGCGCTTAGGGTTTACACTTAACATTACGCTAACGGCACTAATACCATCAGATTCTTTTGTGACTGGTGCTAc                                                                                                                                               >  3:618875/1‑150 (MQ=255)
      gTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAAt                                                                                                                                            <  1:167587/150‑1 (MQ=255)
          tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTa                                                                                                                                                >  2:1849399/1‑142 (MQ=255)
               gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                                                   <  1:1390599/150‑1 (MQ=255)
                    ggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                                                    <  1:808114/144‑1 (MQ=255)
                    ggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCaaa                                                                                                                                    >  2:808114/1‑144 (MQ=255)
                        aTGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtc                                                                                                                          <  1:57170/150‑1 (MQ=255)
                          gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtcac                                                                                                                        >  1:1445698/1‑148 (MQ=255)
                           cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCa                                                                                                                       >  2:1297887/1‑150 (MQ=255)
                            ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCa                                                                                                                                      <  2:681768/134‑1 (MQ=255)
                                 tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                                 <  2:1183573/150‑1 (MQ=255)
                                      tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                                                  >  1:102437/1‑144 (MQ=255)
                                      tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                                                  <  2:102437/144‑1 (MQ=255)
                                      tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCACAAAATCagaaga                                                                                                            <  2:1445698/150‑1 (MQ=255)
                                      tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                            >  3:1214040/1‑150 (MQ=255)
                                       ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                                           <  3:997169/150‑1 (MQ=255)
                                         aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                        >  1:109024/1‑135 (MQ=255)
                                          cTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTGCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAt                                                                                                        <  3:1262011/150‑1 (MQ=255)
                                           tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                                                       >  3:635862/1‑150 (MQ=255)
                                           tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                                                       >  1:845990/1‑150 (MQ=255)
                                             cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                     <  2:887800/150‑1 (MQ=255)
                                                    tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                                              >  2:1227770/1‑150 (MQ=255)
                                                      agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                            <  3:314636/150‑1 (MQ=255)
                                                          gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACaa                                                                                                                             <  4:1203304/113‑1 (MQ=255)
                                                          gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                                                        >  4:577356/1‑150 (MQ=255)
                                                               gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATgg                                                                                               >  1:2540066/1‑138 (MQ=255)
                                                               gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATgg                                                                                               <  2:2540066/138‑1 (MQ=255)
                                                               gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                                                   <  1:471517/150‑3 (MQ=255)
                                                                cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctac                                                                                  >  1:2198397/1‑150 (MQ=255)
                                                                  aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTTTTTTTTTTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacga                                                                                >  4:1272713/1‑150 (MQ=255)
                                                                  aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAagag                                                                                         <  3:143857/141‑1 (MQ=255)
                                                                  aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacga                                                                                >  4:541195/1‑150 (MQ=255)
                                                                   ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                                              >  2:871510/1‑119 (MQ=255)
                                                                   ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGaa                                                                               >  2:1103205/1‑150 (MQ=255)
                                                                   ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGaa                                                                               <  2:549064/150‑1 (MQ=255)
                                                                     aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTa                                                                                                                                                >  1:877796/1‑83 (MQ=255)
                                                                       aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAAt                                                                                                                   >  4:494493/1‑110 (MQ=255)
                                                                       aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAAt                                                                                                                   <  3:494493/110‑1 (MQ=255)
                                                                       aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCtt                                                                           >  2:927329/1‑150 (MQ=255)
                                                                       aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCtt                                                                           >  1:1886244/1‑150 (MQ=255)
                                                                        tGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTg                                                                          <  4:621620/150‑1 (MQ=255)
                                                                        tGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATAACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTg                                                                          <  2:845990/150‑1 (MQ=255)
                                                                         gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGa                                                                         >  2:524944/1‑150 (MQ=255)
                                                                           ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTCTTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                       >  1:1854952/1‑150 (MQ=255)
                                                                           ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                                                       >  1:2333328/1‑118 (MQ=255)
                                                                           ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTaa                                                                                                  >  4:451306/1‑123 (MQ=255)
                                                                           ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                       >  2:723454/1‑150 (MQ=255)
                                                                           ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                       <  1:1945454/150‑1 (MQ=255)
                                                                           ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                       >  1:325918/1‑150 (MQ=255)
                                                                                tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTAGGAAGCTTGAAGCAACg                                                                  >  4:1366807/1‑150 (MQ=255)
                                                                                tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACg                                                                  >  2:1942066/1‑150 (MQ=255)
                                                                                     tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                                              <  3:543840/101‑1 (MQ=255)
                                                                                       aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                                                <  1:2275585/97‑1 (MQ=255)
                                                                                        gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                        >  4:1460384/1‑136 (MQ=255)
                                                                                        gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCaaa                                                          >  1:1964880/1‑150 (MQ=255)
                                                                                         ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                         >  3:769729/1‑150 (MQ=255)
                                                                                         ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTCATGACAAGCTATGTTTCGATGCTTGAAGCAACGGTCGCAAAg                                                         <  4:1166105/150‑1 (MQ=255)
                                                                                           tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                                                   >  2:1594700/1‑138 (MQ=255)
                                                                                           tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                                                   <  1:1594700/138‑1 (MQ=255)
                                                                                           tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                       >  1:644697/1‑150 (MQ=255)
                                                                                            ttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAt                                                                                                                            >  3:1299661/1‑80 (MQ=255)
                                                                                            ttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAt                                                                                                                            <  4:1300284/80‑1 (MQ=255)
                                                                                               cacTCAACATTACGCTAACGGCACTAAAACCATCACATATTTCTGTGACTGGTGCTACAATCTCCCAAAGTCACAATTCTCAACATCAGAAGAGAATTGCTAATGGCACGAGACGCTACGAAGCTTGAAAGCACGGTCGCAAAGCTGaaa                                                   >  2:2509758/1‑150 (MQ=255)
                                                                                                acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAACGTCACAATACTCAAAATCAGAAGAGTATGGCTAATGGCAAGAGACGCTACGAAGCGTGAAGCAACGGTCGCAAAGCTGaaaa                                                  >  3:692690/1‑150 (MQ=255)
                                                                                                acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACACTTATAAAAATACAAAGAGTCAAGCAACTAGGAAGAGACGCAAAGAAGATTGAAGAAAAGGTCGCAAAGCTGaaaa                                                  >  1:373752/1‑150 (MQ=255)
                                                                                                 cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                                              >  2:641198/1‑89 (MQ=255)
                                                                                                 cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                 >  2:265354/1‑150 (MQ=255)
                                                                                                  tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGt                                                                                                                                  >  3:1173589/1‑68 (MQ=255)
                                                                                                  tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                                         <  4:16595/93‑1 (MQ=255)
                                                                                                  tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                                                                                         >  3:16595/1‑93 (MQ=255)
                                                                                                  tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                                >  4:1244281/1‑150 (MQ=255)
                                                                                                    aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                                              >  1:546520/1‑86 (MQ=255)
                                                                                                      cATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                            <  1:1288149/150‑1 (MQ=255)
                                                                                                        ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAt                                                                                                 <  1:306314/95‑1 (MQ=255)
                                                                                                         tACGCTAACGGCACTAAAACCATCACATTTTTTTTTTTTTTTTTCTACAATCTTCCCAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCACGAGACTCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                         >  2:451397/1‑150 (MQ=255)
                                                                                                         tACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                         <  2:311918/150‑1 (MQ=255)
                                                                                                          aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                                        <  1:2420601/150‑1 (MQ=255)
                                                                                                           cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                                       >  1:2031645/1‑150 (MQ=255)
                                                                                                                  ggCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                >  3:367211/1‑150 (MQ=255)
                                                                                                                   gCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                               >  1:1100414/1‑150 (MQ=255)
                                                                                                                   gCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                               >  1:701518/1‑150 (MQ=255)
                                                                                                                    cACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAc                                                                                       >  2:1370174/1‑93 (MQ=255)
                                                                                                                    cACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                              >  2:190209/1‑150 (MQ=255)
                                                                                                                     aCTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTCCAAAAGTACCAATTCACACACTCAGAAGAGAAACGCCTATGGAAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGccc                             <  3:146061/150‑1 (MQ=255)
                                                                                                                       tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg                           >  2:475014/1‑150 (MQ=255)
                                                                                                                         aaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgctac                                                                                  >  2:2279700/1‑90 (MQ=255)
                                                                                                                          aaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                       >  1:1778761/1‑119 (MQ=255)
                                                                                                                          aaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                        >  1:1214287/1‑150 (MQ=255)
                                                                                                                            ccATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                      <  2:325918/150‑1 (MQ=255)
                                                                                                                             cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGat                     >  3:452275/1‑150 (MQ=255)
                                                                                                                             cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGat                     <  3:1244281/150‑1 (MQ=255)
                                                                                                                              aTCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata                    <  1:906035/150‑1 (MQ=255)
                                                                                                                                     ttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc             >  2:2339467/1‑150 (MQ=255)
                                                                                                                                     ttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCACGCGATATgcgcgc             >  4:1228884/1‑150 (MQ=255)
                                                                                                                                         cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCGCGCGCATATCGCGCggggcc       <  1:265354/150‑5 (MQ=255)
                                                                                                                                           gtgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt       >  2:2477952/1‑150 (MQ=255)
                                                                                                                                            tgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCaa                                                                    <  1:2243544/88‑1 (MQ=255)
                                                                                                                                            tgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCaa                                                                    >  2:2243544/1‑88 (MQ=255)
                                                                                                                                               cTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc   >  4:863684/1‑150 (MQ=255)
                                                                                                                                                tGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg  >  1:312692/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCG  >  NC_000913_3_bme_pgi/4233561‑4233854

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