Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,772 G→T intergenic (‑539/‑2135) 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,7720GT99.6% 866.0 / ‑4.2 235intergenic (‑539/‑2135)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (124/110);  total (124/111)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 238 total aligned reads displayed.

TAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGA  >  NC_000913_3_bme_pgi/4233626‑4233917
                                                                                                                                                  |                                                                                                                                                  
 aataaATGCTTCACTGCGCTAAGGGTTGCCGCTCAACATTTCGCTAACGGAACTAAAACCATGACATTTTTCTGGGACAGGCGCTACAATCTTGCAGAGTCACGATGATCGGGATCAGAAGAGTTTTGTTAATGGCAAGGGACGcgacga                                                                                                                                                <  2:95419/150‑1 (MQ=255)
     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                                                                                            <  3:95950/150‑1 (MQ=255)
       tGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTg                                                                                                                                          >  3:285628/1‑150 (MQ=255)
        gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGa                                                                                                                                         <  1:191362/150‑1 (MQ=255)
            cACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCa                                                                                                                                     >  3:704203/1‑150 (MQ=255)
                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCACAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                                                                                                              <  1:1571788/150‑1 (MQ=255)
                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                                                                                                              <  3:481059/150‑1 (MQ=255)
                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                                                                                                              <  3:357051/150‑1 (MQ=255)
                     aaGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCa                                                                                                                            >  3:908380/1‑150 (MQ=255)
                      aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCaa                                                                                                                           >  1:1019555/1‑150 (MQ=255)
                         gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                                                        >  1:2249124/1‑150 (MQ=255)
                         gTTTACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                                                        <  1:2283119/150‑1 (MQ=255)
                           ttACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTg                                                                                                                      >  2:1225114/1‑150 (MQ=255)
                           ttACACTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTg                                                                                                                      <  2:193507/150‑1 (MQ=255)
                               acTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                                                                                  >  1:2167657/1‑150 (MQ=255)
                                cTCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                                                 <  4:315057/150‑1 (MQ=255)
                                 tCAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                                                                                                >  2:328282/1‑150 (MQ=255)
                                  cAACATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                                                                                               <  2:1992481/150‑1 (MQ=255)
                                   aaCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                                                                                              >  2:103339/1‑150 (MQ=255)
                                    aCATTACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAt                                                                                                             >  4:746825/1‑150 (MQ=255)
                                       ttACGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTgg                                                                                                          <  1:258372/150‑1 (MQ=255)
                                          cGCTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                                                                                                       <  2:1422878/150‑1 (MQ=255)
                                           gCTAACGGCACTAAAACCATTACATTTTTCGGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCgg                                                                                                      >  2:1369761/1‑150 (MQ=255)
                                            cTAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                                                         >  3:1271562/1‑146 (MQ=255)
                                             tAACGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGaa                                                                                                    <  1:1423609/150‑1 (MQ=255)
                                               aCGGCACTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                  >  1:1253326/1‑150 (MQ=255)
                                                     cTAAAACCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAGGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcccc                                                                                            >  3:1516742/1‑150 (MQ=255)
                                                        aaaCCATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcg                                                                                         >  4:1581713/1‑150 (MQ=255)
                                                           ccATTACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                                                     <  3:116677/150‑1 (MQ=255)
                                                               tacaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATg                                                                                 >  1:281704/2‑150 (MQ=255)
                                                                  aTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgc                                                                              >  4:356696/1‑150 (MQ=255)
                                                                          gtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaa                                                                                                                    <  2:1913029/105‑1 (MQ=255)
                                                                          gtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt                                                                      >  2:1267480/1‑150 (MQ=255)
                                                                              cTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc                                                                  >  2:1137203/1‑150 (MQ=255)
                                                                               tGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                                                 <  2:1253326/150‑1 (MQ=255)
                                                                                ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcgc                                                                <  2:1539318/150‑1 (MQ=255)
                                                                                ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcgc                                                                <  1:87453/150‑1 (MQ=255)
                                                                                 gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                               >  1:1973641/1‑150 (MQ=255)
                                                                                  cgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGc                                                                         <  4:914999/139‑1 (MQ=255)
                                                                                   gcTACACTCTTCCAAAGTCCCAATTCTCAAAATCAGAAGAGTAATGCTAATGGCAAGAGACGATACGAAGCTTGAAGCCACGGTCGCAAAGATGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGAATTCAGCGCCGa                                                             >  2:467427/1‑150 (MQ=255)
                                                                                   gcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGcc                                                                                              >  2:165802/1‑118 (MQ=255)
                                                                                     tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGaa                                                                                                    <  4:1426706/110‑1 (MQ=255)
                                                                                     tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAcc                                                           >  2:2567090/1‑150 (MQ=255)
                                                                                        aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                        <  1:1216601/150‑1 (MQ=255)
                                                                                         aTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaa                                                                                                                   <  2:473085/91‑1 (MQ=255)
                                                                                         aTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaa                                                                                                                   >  1:473085/1‑91 (MQ=255)
                                                                                            ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCGCGCGCATATCGCGCggggccg                                                                     <  4:1384124/131‑6 (MQ=255)
                                                                                            ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCGCGCGCATATCGCGCggggccg                                                                     >  3:1384124/1‑126 (MQ=255)
                                                                                             tCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccgg                                                   >  3:956809/1‑150 (MQ=255)
                                                                                              ccAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                                                         >  2:599052/1‑96 (MQ=255)
                                                                                              ccAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                                                         <  1:599052/96‑1 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                                                     >  4:529160/1‑113 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                                                                                     <  3:529160/113‑1 (MQ=255)
                                                                                                  aGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCgg                                              <  1:2481223/150‑1 (MQ=255)
                                                                                                  aGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCgg                                              <  3:746825/150‑1 (MQ=255)
                                                                                                    tCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCgggg                                                                        >  3:165525/1‑120 (MQ=255)
                                                                                                    tCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCgggg                                                                        <  4:165525/122‑3 (MQ=255)
                                                                                                      acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                   <  4:1055706/94‑1 (MQ=255)
                                                                                                      acaATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                   >  3:1055706/1‑94 (MQ=255)
                                                                                                       caATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        <  4:1603167/56‑1 (MQ=255)
                                                                                                       caATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        >  3:1603167/1‑56 (MQ=255)
                                                                                                       caATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                                                                                             >  1:501551/1‑67 (MQ=255)
                                                                                                        aaTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata                                        >  1:239965/1‑150 (MQ=255)
                                                                                                        aaTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGata                                        >  2:763683/1‑150 (MQ=255)
                                                                                                             tcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCATATCGCGCgggg                                                                        <  3:1144468/111‑3 (MQ=255)
                                                                                                              cAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct                                  <  3:969966/150‑1 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc                                 >  4:1228103/1‑150 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc                                 >  2:2075162/1‑150 (MQ=255)
                                                                                                                  aTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGTCCGGTTCGGGCGAGAGTCGCTCTGc                              >  2:916630/1‑150 (MQ=255)
                                                                                                                     agaagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                           >  1:954330/1‑150 (MQ=255)
                                                                                                                     agaagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                           <  1:328282/150‑1 (MQ=255)
                                                                                                                     agaagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTg                           <  2:2464195/150‑1 (MQ=255)
                                                                                                                        agagTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgat                        <  4:956809/150‑1 (MQ=255)
                                                                                                                         gagTATTGGTTATGGAAAGCGGCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGGGCGGTATGCGCGCGGCGTTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgatg                       <  2:691744/150‑1 (MQ=255)
                                                                                                                         gagTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCggg                                                                         <  3:301230/100‑2 (MQ=255)
                                                                                                                          agTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt                                                                      >  3:1450271/1‑102 (MQ=255)
                                                                                                                          agTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGgatga                      <  1:1995164/150‑1 (MQ=255)
                                                                                                                           gTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAAt                                                                                                   <  2:2128440/73‑1 (MQ=255)
                                                                                                                            tATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCgg                                              >  1:843239/1‑124 (MQ=255)
                                                                                                                            tATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtc                    <  1:1677826/150‑1 (MQ=255)
                                                                                                                             aTTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCAtt                                                                      <  1:1307372/99‑1 (MQ=255)
                                                                                                                                gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtcttct                <  2:553374/150‑1 (MQ=255)
                                                                                                                                 cTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                                                  >  1:1712452/1‑68 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCa                                                                        >  1:1757307/1‑92 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCa                                                                        <  2:1757307/92‑1 (MQ=255)
                                                                                                                                   aaTGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                                <  3:857496/83‑1 (MQ=255)
                                                                                                                                   aaTGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGtt             >  2:747112/1‑150 (MQ=255)
                                                                                                                                      ggCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCg                                                                                                       <  2:2171239/58‑1 (MQ=255)
                                                                                                                                       gCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGAc         >  1:1399526/1‑150 (MQ=255)
                                                                                                                                       gCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGAc         <  4:887378/150‑1 (MQ=255)
                                                                                                                                           gagaCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCg                                               <  1:1074837/108‑1 (MQ=255)
                                                                                                                                           gagaCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCg                                               >  2:1074837/1‑108 (MQ=255)
                                                                                                                                           gagaCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGt     <  4:1195698/150‑1 (MQ=255)
                                                                                                                                           gagaCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGt     <  2:2249124/150‑1 (MQ=255)
                                                                                                                                              aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                                                         <  4:649260/48‑1 (MQ=255)
                                                                                                                                              aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                                                         >  3:649260/1‑48 (MQ=255)
                                                                                                                                              aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                                                                                                        >  1:688554/1‑49 (MQ=255)
                                                                                                                                              aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGa  <  4:293469/150‑1 (MQ=255)
                                                                                                                                              aCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGa  >  2:146124/1‑150 (MQ=255)
                                                                                                                                               cgctACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                                                >  4:856220/1‑56 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGA  >  NC_000913_3_bme_pgi/4233626‑4233917

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: