Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,689 C→T intergenic (‑456/‑2218) 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,6890CT99.5% 773.8 / ‑5.1 212intergenic (‑456/‑2218)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (110/100);  total (111/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 213 total aligned reads displayed.

TTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGA  >  NC_000913_3_bme_pgi/4233542‑4233834
                                                                                                                                                   |                                                                                                                                                 
     aataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAtttt                                                                                                                                            >  4:289862/1‑150 (MQ=255)
           gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCAGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCtgtg                                                                                                                                      <  3:429431/150‑1 (MQ=255)
            ggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAttttt                                                                                                                                           >  4:1427671/1‑144 (MQ=255)
            ggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGa                                                                                                                                     <  2:239122/150‑1 (MQ=255)
             gTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAAGCAGTGCATTTTTCTGTGAc                                                                                                                                    >  1:691744/1‑150 (MQ=255)
                  cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGgcg                                                                                                                               >  4:357051/1‑150 (MQ=255)
                      gCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTAc                                                                                                                           <  2:392574/150‑1 (MQ=255)
                        ggTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACaa                                                                                                                         <  2:110992/150‑1 (MQ=255)
                          ttGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATc                                                                                                                       <  3:212624/150‑1 (MQ=255)
                           tGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCt                                                                                                                      <  1:2069300/150‑1 (MQ=255)
                            gTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGTCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCtt                                                                                                                     <  1:2049729/150‑1 (MQ=255)
                             tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTc                                                                                                                    >  2:1381213/1‑150 (MQ=255)
                              cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCTGTCATTATTTCATCTAATAAAGGCGTCACTGCTCGAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTcc                                                                                                                   >  1:1050514/1‑150 (MQ=255)
                              cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACAtt                                                                                                                                              >  4:375612/1‑123 (MQ=255)
                              cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTcc                                                                                                                   <  3:925048/150‑1 (MQ=255)
                               aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTAc                                                                                                                           >  2:1205881/1‑141 (MQ=255)
                               aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTAc                                                                                                                           <  1:1205881/141‑1 (MQ=255)
                               aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                                                                  >  1:2603773/1‑150 (MQ=255)
                                         tCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAAt                                                                                                        <  2:2514599/150‑1 (MQ=255)
                                           aTGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtc                                                                                                      <  4:661482/150‑1 (MQ=255)
                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                                    <  4:1262438/150‑1 (MQ=255)
                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACt                                                                                                                                   <  3:170418/117‑1 (MQ=255)
                                                         tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        >  1:2449783/1‑150 (MQ=255)
                                                         tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        <  1:723275/150‑1 (MQ=255)
                                                         tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                          <  2:1973093/148‑1 (MQ=255)
                                                         tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                          >  1:1973093/1‑148 (MQ=255)
                                                          ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                                                                            <  2:1720272/145‑1 (MQ=255)
                                                              tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                                                                   <  2:1433777/150‑1 (MQ=255)
                                                               ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTg                                                                                  <  2:61587/150‑1 (MQ=255)
                                                                 tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                                                                 >  1:1473454/1‑133 (MQ=255)
                                                                    cGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAAt                                                                             <  1:1384993/150‑1 (MQ=255)
                                                                       tCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                  >  1:1915369/1‑126 (MQ=255)
                                                                         agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                        >  1:371643/1‑150 (MQ=255)
                                                                         agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                        >  1:553374/1‑150 (MQ=255)
                                                                          gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAag                                                                       >  4:1102647/1‑150 (MQ=255)
                                                                          gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAag                                                                       <  3:1225403/150‑1 (MQ=255)
                                                                           tagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTgg                                                                                                                                 <  3:1264335/91‑1 (MQ=255)
                                                                             gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtct                                                                                                     <  3:1194523/117‑1 (MQ=255)
                                                                                cAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAATTCAGAAGAGTATTGCTAATGGCAAGAGAcgc                                                                 <  1:2579638/150‑1 (MQ=255)
                                                                                  gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                                <  4:709770/101‑1 (MQ=255)
                                                                                  gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                                >  3:709770/1‑101 (MQ=255)
                                                                                  gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTAAAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgcta                                                               >  1:1894807/1‑148 (MQ=255)
                                                                                   caaatcaATGCTTCACAGCGCACAGCGTTTCCACTAAACATGAAGCTAACGGCACTAAAACCATCACCTTTTTCCGTGACTAGCGCTACAATCTTCCAAAGCCACCATCCTCAAAATCAGAAGAGCATCGCTAATGGCACGAGAcgcccc                                                              >  2:2481223/3‑147 (MQ=255)
                                                                                   cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        <  2:202714/124‑1 (MQ=255)
                                                                                   cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAcgctac                                                              >  4:116677/1‑147 (MQ=255)
                                                                                    taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATgg                                                                           <  2:72201/136‑1 (MQ=255)
                                                                                    taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                             >  3:1362602/1‑150 (MQ=255)
                                                                                     aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                                    <  1:2414826/110‑1 (MQ=255)
                                                                                     aataaATGCTTCACTGCGCTAAGGGTTGCCGCTCAACATTTCGCTAACGGAACTAAAACCATGACATTTTTCTGGGACAGGCGCTACAATCTTGCAGAGTCACGATGATCGGGATCAGAAGAGTTTTGTTAATGGCAAGGGACGcgacga                                                            <  2:95419/150‑1 (MQ=255)
                                                                                        aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCtg                                                                                                                                        >  2:1921522/1‑71 (MQ=255)
                                                                                        aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCtg                                                                                                                                        <  1:1921522/71‑1 (MQ=255)
                                                                                         aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                        <  1:891435/150‑1 (MQ=255)
                                                                                             cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                    <  4:892450/150‑1 (MQ=255)
                                                                                             cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                    <  2:1894807/150‑1 (MQ=255)
                                                                                              ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCaa                                                                        <  1:1008746/129‑1 (MQ=255)
                                                                                              ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                   <  2:472937/150‑1 (MQ=255)
                                                                                               tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGc                                                  >  3:1074010/1‑150 (MQ=255)
                                                                                               tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGc                                                  >  3:1459398/1‑150 (MQ=255)
                                                                                                cACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGa                                                                                                                                     >  2:2062097/1‑66 (MQ=255)
                                                                                                cACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGa                                                                                                                                     <  1:2062176/66‑1 (MQ=255)
                                                                                                cACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCa                                                 >  3:704203/1‑150 (MQ=255)
                                                                                                 aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                  <  2:2158814/100‑1 (MQ=255)
                                                                                                 aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                                                                  >  1:2158814/1‑100 (MQ=255)
                                                                                                 aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCaa                                                <  4:208284/150‑1 (MQ=255)
                                                                                                  cTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                               <  1:2485539/150‑1 (MQ=255)
                                                                                                   tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                                                                            <  2:2432806/88‑1 (MQ=255)
                                                                                                   tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                                    <  3:278201/96‑1 (MQ=255)
                                                                                                       cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCACAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                          <  1:1571788/150‑1 (MQ=255)
                                                                                                       cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                          <  3:481059/150‑1 (MQ=255)
                                                                                                       cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                          <  3:357051/150‑1 (MQ=255)
                                                                                                        tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        >  2:2318842/1‑103 (MQ=255)
                                                                                                        tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        <  1:2318842/103‑1 (MQ=255)
                                                                                                        tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                         >  2:2283119/1‑150 (MQ=255)
                                                                                                          aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCaa                                       >  3:794373/1‑150 (MQ=255)
                                                                                                             gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                                                              >  3:195989/1‑92 (MQ=255)
                                                                                                             gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                    >  1:2249124/1‑150 (MQ=255)
                                                                                                             gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                    <  1:2283119/150‑1 (MQ=255)
                                                                                                               ttACACTCAACATTACGCTAACGGCACTAAAACCATTACAt                                                                                                                                               <  1:1247877/41‑1 (MQ=255)
                                                                                                               ttACACTCAACATTACGCTAACGGCACTAAAACCATTACAt                                                                                                                                               >  2:1247877/1‑41 (MQ=255)
                                                                                                                 acacTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaa                                <  1:450366/150‑1 (MQ=255)
                                                                                                                   acTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                              >  1:1422878/1‑150 (MQ=255)
                                                                                                                     tccaCATTACGCTAACGACACTAAACCCATTACATTTTTCTGTGACTGGCGCTACCATCTTCCAACGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTAAAGCAACGGTCGCAAAGCTGAAAAAg                            >  2:102409/4‑150 (MQ=255)
                                                                                                                     tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                            >  2:328282/1‑150 (MQ=255)
                                                                                                                      cAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                           <  2:1992481/150‑1 (MQ=255)
                                                                                                                          aTTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTATCAAAATCAGAAGAGTATTGCTAATGGCAAGAGAAGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                       >  2:2118381/1‑150 (MQ=255)
                                                                                                                             aCGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                                    <  2:2352337/70‑1 (MQ=255)
                                                                                                                             aCGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                                                                    >  1:2352337/1‑70 (MQ=255)
                                                                                                                             aCGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                    >  4:188460/1‑150 (MQ=255)
                                                                                                                              cGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                   <  2:1422878/150‑1 (MQ=255)
                                                                                                                               gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                          <  2:1339528/78‑1 (MQ=255)
                                                                                                                               gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaa                                                                                          >  1:1339528/1‑78 (MQ=255)
                                                                                                                               gCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                  >  3:293469/1‑150 (MQ=255)
                                                                                                                                cTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                     <  4:1271562/146‑1 (MQ=255)
                                                                                                                                 tAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGaa                >  1:1755091/1‑150 (MQ=255)
                                                                                                                                       cACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        <  2:1641437/72‑1 (MQ=255)
                                                                                                                                       cACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                                                                        >  1:1641437/1‑72 (MQ=255)
                                                                                                                                         cTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAGGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc        >  3:1516742/1‑150 (MQ=255)
                                                                                                                                            aaaCCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg     >  4:1581713/1‑150 (MQ=255)
                                                                                                                                            aaaCCATTACATTTTTCTGTGACTGGCGCTACAATCTTCAAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg     >  4:959018/1‑150 (MQ=255)
                                                                                                                                              aCCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgcg   <  1:103339/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGA  >  NC_000913_3_bme_pgi/4233542‑4233834

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: