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,7720GT98.9% 345.3 / NA 95intergenic (‑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 (50/43);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCG  >  NC_000913_3_bme_pgi/4233626‑4233916
                                                                                                                                                  |                                                                                                                                                  
taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGctacg                                                                                                                                                 >  1:770648/1‑150 (MQ=255)
     aaTGCTTCACTGCGCTAATGGTTTACGCTCAACGTGACGCTAACTTCACTATTTTCAGCACATTTGTCTGGGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                                                                                            <  3:228732/150‑1 (MQ=255)
     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                                                                                            >  4:192720/1‑150 (MQ=255)
     aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCt                                                                                                                                            >  2:696275/1‑150 (MQ=255)
        gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGa                                                                                                                                         <  3:661084/150‑1 (MQ=255)
         cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        <  3:209731/150‑1 (MQ=255)
         cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        >  4:209731/1‑150 (MQ=255)
          ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAg                                                                                                                                       <  4:393303/150‑1 (MQ=255)
              cTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                                                                                                                   <  4:473550/150‑1 (MQ=255)
               tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACg                                                                                                                                  >  2:340071/1‑150 (MQ=255)
                 cgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGt                                                                                                                                >  1:1020812/1‑150 (MQ=255)
                   cTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCg                                                                                                                              <  3:308218/150‑1 (MQ=255)
                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        >  2:725474/1‑139 (MQ=255)
                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGaa                                                                                                                                        <  1:725474/139‑1 (MQ=255)
                    tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                                                                                             >  2:718151/1‑150 (MQ=255)
                     aaGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCa                                                                                                                            >  2:60033/1‑150 (MQ=255)
                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                                                         >  2:733070/1‑150 (MQ=255)
                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                                                         >  4:586594/1‑150 (MQ=255)
                        ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                                                         <  2:289449/150‑1 (MQ=255)
                               acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                                                                                  >  1:242687/1‑150 (MQ=255)
                               acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                                                                                  >  1:242688/1‑150 (MQ=255)
                               acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                                                                                  >  1:259096/1‑150 (MQ=255)
                                 tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaaa                                                                                                                >  2:439068/1‑149 (MQ=255)
                                  cAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaaac                                                                                                               >  3:643999/1‑148 (MQ=255)
                                         aCGCTNACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGc                                                                                                        <  3:197680/150‑1 (MQ=255)
                                         aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGc                                                                                                        <  2:903838/150‑1 (MQ=255)
                                          cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCg                                                                                                       <  2:242687/150‑1 (MQ=255)
                                            cTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGa                                                                                                     <  1:60033/150‑1 (MQ=255)
                                                cGGCACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCg                                                                                                 >  4:226365/1‑150 (MQ=255)
                                                   cACTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGcc                                                                                              >  3:437061/1‑150 (MQ=255)
                                                     cTAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACACTTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGcccc                                                                                            >  1:727028/1‑150 (MQ=255)
                                                      tacaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCCAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCcg                                                                                           >  2:1012133/4‑150 (MQ=255)
                                                      tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATc                                                                                                  >  3:118805/1‑143 (MQ=255)
                                                      tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATc                                                                                                  <  4:118805/143‑1 (MQ=255)
                                                      tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCcg                                                                                           >  3:41368/1‑150 (MQ=255)
                                                         aaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                                                                                        >  4:300348/1‑150 (MQ=255)
                                                          aCCATCACATTTCTCTTTGACTGGTGCTACAATTTTCCAAAGGCACGATTCACAAAATCAGAAGAGTGTTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCCACGGCCGCAAAGCTGACAAAACATTGGGCGGAATCGGCCCcgcgcg                                                                                       <  2:267660/150‑1 (MQ=255)
                                                           ccATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGa                                                                                     <  3:192720/150‑1 (MQ=255)
                                                           ccATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGa                                                                                     <  2:770648/150‑1 (MQ=255)
                                                            cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGat                                                                                    >  2:583125/1‑150 (MQ=255)
                                                                acaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                               >  4:385796/1‑150 (MQ=255)
                                                                acaTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATAGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATgc                                                                               >  3:179406/1‑150 (MQ=255)
                                                                 caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                                                                                             >  1:1444/1‑105 (MQ=255)
                                                                 caTTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGc                                                                                                                             <  2:1444/105‑1 (MQ=255)
                                                                    ttttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc                                                                           >  3:414801/1‑150 (MQ=255)
                                                                     tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg                                                                          <  3:456605/150‑1 (MQ=255)
                                                                      ttCTGTGACTGGTGCTACAAGCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCgg                                                                         <  2:615171/150‑1 (MQ=255)
                                                                        cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCGCGCGCATATCGCGCggagcc                                                                     <  2:1020812/150‑5 (MQ=255)
                                                                         tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCggggc                                                                      <  2:504694/150‑4 (MQ=255)
                                                                          gtgACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCGCGCGCATATCGCGCggggccga                                                                   <  1:718151/150‑7 (MQ=255)
                                                                                 gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaa                                                                                                                    >  2:91850/3‑98 (MQ=255)
                                                                                 gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGaa                                                                                                                    <  1:91850/96‑1 (MQ=255)
                                                                                 gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                              <  1:583125/148‑1 (MQ=255)
                                                                                 gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                              >  3:414663/3‑150 (MQ=255)
                                                                                     tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACgg                                                                                                                                 <  1:583173/81‑1 (MQ=255)
                                                                                     tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACgg                                                                                                                                 <  1:583239/81‑1 (MQ=255)
                                                                                     tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACgg                                                                                                                                 >  2:583239/1‑81 (MQ=255)
                                                                                     tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACgg                                                                                                                                 >  2:583173/1‑81 (MQ=255)
                                                                                        aaTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGTGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGGCGACCccg                                                       >  2:881023/1‑150 (MQ=255)
                                                                                         aTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                              <  3:143121/142‑1 (MQ=255)
                                                                                         aTCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                              >  4:143121/1‑142 (MQ=255)
                                                                                            ttCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccg                                                   >  3:264172/1‑150 (MQ=255)
                                                                                              ccAAAGGCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATAGGCGCGCGGCGCTCAGCGCCGACCCCGGCCGGt                                                 >  2:322198/1‑150 (MQ=255)
                                                                                               cAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGtt                                                <  3:385796/150‑1 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                               <  4:55764/150‑1 (MQ=255)
                                                                                                aaaGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                                               >  3:55764/1‑150 (MQ=255)
                                                                                                  aGTCACAATTCTCAAAATCAGAAGAGTATTGCGAAGGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGCCATTCAGCGCCGACCCCGGCCGGTTCgg                                             <  2:991147/150‑1 (MQ=255)
                                                                                                    acacaATTATCAAAATGAAAAGAGTATTGATAATGGCAAGAGACGCAACGAAGCTTGAAGCAACTGTCGAAAAACAGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                                           <  4:643999/149‑1 (MQ=255)
                                                                                                          ttCTCAAAATCAGAAGAGTAGTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATAtt                                     <  2:415126/150‑1 (MQ=255)
                                                                                                             tcAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc                                  >  1:1081942/1‑150 (MQ=255)
                                                                                                              cAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGTGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct                                 <  1:439068/150‑1 (MQ=255)
                                                                                                              cAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct                                 >  1:323530/1‑150 (MQ=255)
                                                                                                              cAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct                                 >  3:569869/1‑150 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc                                >  2:527641/1‑150 (MQ=255)
                                                                                                               aaaaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc                                >  4:356847/1‑150 (MQ=255)
                                                                                                                 aaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                                                                                                                   <  3:105062/51‑1 (MQ=255)
                                                                                                                 aaTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAAc                                                                                                                                   >  4:105062/1‑51 (MQ=255)
                                                                                                                  aTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc                             <  3:441949/150‑1 (MQ=255)
                                                                                                                       aagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGga                        >  1:1079723/1‑150 (MQ=255)
                                                                                                                             aTTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtct                  <  2:74735/150‑1 (MQ=255)
                                                                                                                                gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                       <  3:659489/110‑1 (MQ=255)
                                                                                                                                gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                       >  4:659489/1‑110 (MQ=255)
                                                                                                                                gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                         <  2:88931/124‑1 (MQ=255)
                                                                                                                                gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                         >  1:88931/1‑124 (MQ=255)
                                                                                                                                gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCACAAAGCTGAAAAAACCTTGGGCGGACTCGGCCCCGCGCGATACGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtcttct               <  1:322198/150‑1 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCcgcg                                                                                         <  2:755219/76‑1 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCcgcg                                                                                         >  1:755219/1‑76 (MQ=255)
                                                                                                                                  tAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGt             <  3:586594/150‑1 (MQ=255)
                                                                                                                                      ggCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                       <  2:807021/104‑1 (MQ=255)
                                                                                                                                      ggCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                       >  1:807021/1‑104 (MQ=255)
                                                                                                                                       gCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGAc        >  2:214627/1‑150 (MQ=255)
                                                                                                                                       gCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGAc        <  2:1079723/150‑1 (MQ=255)
                                                                                                                                       gCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGAc        <  1:733070/150‑1 (MQ=255)
                                                                                                                                       gCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGACATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGAc        >  4:205716/1‑150 (MQ=255)
                                                                                                                                        cAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCGGGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGa         >  1:175654/1‑150 (MQ=255)
                                                                                                                                             gaCGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCg  <  1:1012133/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCG  >  NC_000913_3_bme_pgi/4233626‑4233916

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: