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,7080CT97.8% 467.4 / ‑0.4 134intergenic (‑475/‑2199)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (71/60);  total (73/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 134 total aligned reads displayed.

GCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTC  >  NC_000913_3_bme_pgi/4233564‑4233850
                                                                                                                                                |                                                                                                                                              
gCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTAc                                                                                                                                           >  1:527762/1‑150 (MQ=255)
   cttGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAAt                                                                                                                                        >  4:839681/2‑150 (MQ=255)
       tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                                    <  1:1045268/150‑1 (MQ=255)
        cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTcc                                                                                                                                   >  4:448538/1‑150 (MQ=255)
         aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCa                                                                                                                                  <  3:71656/150‑1 (MQ=255)
            gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                                               >  2:385685/1‑150 (MQ=255)
            gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                                               >  2:1396811/1‑150 (MQ=255)
            gATGGGGGCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                                               >  4:403413/1‑150 (MQ=255)
                gggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTcac                                                                                                                           <  2:932895/150‑1 (MQ=255)
                 ggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTcaca                                                                                                                          >  1:626390/1‑150 (MQ=255)
                  gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACaa                                                                                                                         >  2:69885/1‑150 (MQ=255)
                   tCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAt                                                                                                                        >  1:1112064/1‑150 (MQ=255)
                    cATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAAtt                                                                                                                       <  1:69885/150‑1 (MQ=255)
                       gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                    <  2:1459435/150‑1 (MQ=255)
                       gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATtctc                                                                                                                    >  2:453309/1‑150 (MQ=255)
                        cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                                       >  2:1352926/1‑130 (MQ=255)
                        cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATc                                                                                                                                       <  1:1352926/130‑1 (MQ=255)
                         ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                                  >  3:851803/1‑150 (MQ=255)
                         ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                                  >  1:1593200/1‑150 (MQ=255)
                          gATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                                                                                 <  2:1593200/150‑1 (MQ=255)
                              tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                             >  4:549723/1‑150 (MQ=255)
                               ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                                            >  1:455564/1‑150 (MQ=255)
                                tCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCaga                                                                                                           <  1:264909/150‑1 (MQ=255)
                                   tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                        <  4:1018140/150‑1 (MQ=255)
                                   tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                                        >  3:647277/1‑150 (MQ=255)
                                          cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                 <  2:304563/150‑1 (MQ=255)
                                          cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                 >  1:1541561/1‑150 (MQ=255)
                                          cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                 >  4:204501/1‑150 (MQ=255)
                                           tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCt                                                                                                <  3:344195/150‑1 (MQ=255)
                                              cGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAt                                                                                             <  3:49649/150‑1 (MQ=255)
                                              cGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAt                                                                                             <  2:1512025/150‑1 (MQ=255)
                                                 tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                                          >  1:485974/1‑150 (MQ=255)
                                                   agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTa                                                                                                                                            <  1:250538/98‑1 (MQ=255)
                                                   agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                        >  3:328424/1‑150 (MQ=255)
                                                      agtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAagag                                                                                     <  2:485974/150‑1 (MQ=255)
                                                       gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAgaga                                                                                    <  1:1247300/150‑1 (MQ=255)
                                                           aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAAAATTACGCTAACGGAACTAAAACCATTACATTTTTCTGTGAATGGTGCTACAATCTTCCAAAGTCACAATTCTCAAACTCAGACGACTATTGCAAATGACCACAGAcgct                                                                                >  4:874605/1‑149 (MQ=255)
                                                            gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                                               >  1:232638/1‑148 (MQ=255)
                                                             cTAATAAATGCTTCACTGCGCTAAGGGACTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGGATTGCTAATGGCAAGAGACGctac                                                                              >  1:1104860/1‑150 (MQ=255)
                                                              taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGACGAGTATTGCTAATGGCAAGAGACGctacg                                                                             >  1:774133/1‑150 (MQ=255)
                                                              taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                                             >  3:754375/1‑150 (MQ=255)
                                                              taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                                             >  3:332576/1‑150 (MQ=255)
                                                                ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTTTTTTTTTTTTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGGGTATTGCTAATGGCAAGAGACGCTACGaa                                                                           >  2:1041657/1‑150 (MQ=255)
                                                                ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGaa                                                                           <  2:1310406/150‑1 (MQ=255)
                                                                agaaaTGCTTCACTGCGCGAAGGGTTTACTCTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACCATCTTCCAAAGTCACAATTCTCAAAAGCAGTAGAGTATTGCTAATGGCAAGAGACGCTACGaa                                                                           <  1:1436926/148‑1 (MQ=255)
                                                                    aTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCtt                                                                       >  1:116656/1‑150 (MQ=255)
                                                                       cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                    <  2:527762/150‑1 (MQ=255)
                                                                       cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                    <  1:518526/150‑1 (MQ=255)
                                                                        ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                   <  2:1112064/150‑1 (MQ=255)
                                                                        ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                   >  3:600930/1‑150 (MQ=255)
                                                                                  tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTa                                                                                                                                            >  2:838345/1‑67 (MQ=255)
                                                                                  tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTa                                                                                                                                            <  1:838345/67‑1 (MQ=255)
                                                                                  tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                         >  2:466516/1‑150 (MQ=255)
                                                                                  tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                         >  1:1175400/1‑150 (MQ=255)
                                                                                     gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGc                                                                                          >  1:1327401/1‑114 (MQ=255)
                                                                                      ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                     <  1:1155084/150‑1 (MQ=255)
                                                                                       gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                                    >  2:1027921/1‑150 (MQ=255)
                                                                                          tACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                                       >  1:715186/1‑96 (MQ=255)
                                                                                          tACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAagag                                                                                                       <  2:715186/96‑1 (MQ=255)
                                                                                              cTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                             >  2:1377389/1‑150 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                                    >  2:1095780/1‑62 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                                    <  1:1095780/62‑1 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                            >  3:16057/1‑150 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                            >  2:952738/1‑150 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                            >  1:952298/1‑150 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                            >  1:106510/1‑150 (MQ=255)
                                                                                               tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                            >  4:44735/1‑150 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                                              >  2:1037603/1‑82 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                                              <  1:1037603/82‑1 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                                               >  1:1404255/1‑111 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                                               <  2:1404255/113‑3 (MQ=255)
                                                                                                 aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                          >  4:361486/1‑150 (MQ=255)
                                                                                                    aTTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                                       <  1:1282582/150‑1 (MQ=255)
                                                                                                        cGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                                   <  1:457024/150‑1 (MQ=255)
                                                                                                         gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAag                                                                                       >  3:908103/1‑97 (MQ=255)
                                                                                                         gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                                  <  2:467186/150‑1 (MQ=255)
                                                                                                             aCGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                              <  3:397914/150‑1 (MQ=255)
                                                                                                             aCGGCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                              <  4:229868/150‑1 (MQ=255)
                                                                                                               ggCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                 >  4:553006/1‑81 (MQ=255)
                                                                                                               ggCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                 <  3:553006/81‑1 (MQ=255)
                                                                                                               ggCACTAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                            <  1:130989/150‑1 (MQ=255)
                                                                                                                    tAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                        >  3:109162/1‑85 (MQ=255)
                                                                                                                    tAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                                        <  4:109162/85‑1 (MQ=255)
                                                                                                                    tAAAACCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcg                       >  3:281588/1‑150 (MQ=255)
                                                                                                                      aaaCCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                          >  2:1492702/1‑145 (MQ=255)
                                                                                                                      aaaCCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                          <  1:1492702/145‑1 (MQ=255)
                                                                                                                       aaCCGTAACATTTTCTTGTGACTGGTGCTATAATCTTGCAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTAAGAAGCTTGAAGAAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGTCAcgcgc                    <  4:62146/150‑1 (MQ=255)
                                                                                                                        aCCATTACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgcg                   >  4:229167/1‑150 (MQ=255)
                                                                                                                            ttacaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCGAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGatat               >  2:1009764/3‑150 (MQ=255)
                                                                                                                            ttacaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGatat               <  2:1359787/148‑1 (MQ=255)
                                                                                                                                 tttttCTGTGACCGGTCCTACAATCTTCCACAGCCACAATCCTCAAAATCAGAAGAGTATTGCTAATGGCAACAGACGCTACTAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcg          <  4:332576/150‑1 (MQ=255)
                                                                                                                                   tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg        <  1:1377389/150‑1 (MQ=255)
                                                                                                                                   tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg        <  4:341311/150‑1 (MQ=255)
                                                                                                                                   tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg        <  1:79118/150‑1 (MQ=255)
                                                                                                                                      cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCgggg     <  2:106510/150‑3 (MQ=255)
                                                                                                                                       tgtgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAt    >  1:693669/1‑150 (MQ=255)
                                                                                                                                       tgtgACTGGTGCTACAATCTTCCAAAGTCACAACTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAACTGAAAAAGCATAGGGCGGAACCGGCCACGCTCGATATGGGAGCGGCAt    >  2:230534/1‑150 (MQ=255)
                                                                                                                                         tgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                                              <  2:866710/42‑1 (MQ=255)
                                                                                                                                         tgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCGAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTc  <  1:1594456/150‑1 (MQ=255)
                                                                                                                                           aCTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAg                                                                          >  2:1633356/1‑76 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTC  >  NC_000913_3_bme_pgi/4233564‑4233850

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