Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,808 G→A intergenic (‑575/‑2099) 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,8080GA93.1% 268.0 / 7.0 87intergenic (‑575/‑2099)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (39/42);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

ACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATGGCGCTGTTGA  >  NC_000913_3_bme_pgi/4233667‑4233953
                                                                                                                                             |                                                                                                                                                   
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                                                                                                                                       <  4:118805/143‑1 (MQ=255)
             tAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATc                                                                                                                                       >  3:118805/1‑143 (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)
                           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)
                                        gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                                                                   <  1:583125/148‑1 (MQ=255)
                                        gtgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                                                                                   >  3:414663/3‑150 (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)
                                                                         aTCAGAAGAGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc                                                                  <  3:441949/150‑1 (MQ=255)
                                                                              aagaGTATTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGga                                                             >  1:1079723/1‑150 (MQ=255)
                                                                                    aTTGCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGAtct                                                       <  2:74735/150‑1 (MQ=255)
                                                                                       gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                                                            >  4:659489/1‑110 (MQ=255)
                                                                                       gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccg                                                                                            <  3:659489/110‑1 (MQ=255)
                                                                                       gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                                                              >  1:88931/1‑124 (MQ=255)
                                                                                       gCTAATGGCAAGAGACGCTACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa                                                                              <  2:88931/124‑1 (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)
                                                                                                        ctacgaAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGt                                   <  1:696275/148‑1 (MQ=255)
                                                                                                         tacgaAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTg                                  <  4:264172/149‑1 (MQ=255)
                                                                                                             aaGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCTATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCgg                              <  2:727028/150‑1 (MQ=255)
                                                                                                             aaGCTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCgg                              >  1:324351/1‑150 (MQ=255)
                                                                                                                catGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTGCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGtt                           <  2:719352/148‑1 (MQ=255)
                                                                                                                cTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGtt                           <  4:179406/150‑1 (MQ=255)
                                                                                                                cTTGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGtt                           <  4:569869/150‑1 (MQ=255)
                                                                                                                 ttGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCgg                              >  4:85006/1‑146 (MQ=255)
                                                                                                                 ttGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCgg                              <  3:85006/146‑1 (MQ=255)
                                                                                                                  tGAAGCAACGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTaa                         >  2:132392/1‑150 (MQ=255)
                                                                                                                        aaCGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga                   >  3:44584/1‑150 (MQ=255)
                                                                                                                        aaCGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga                   >  3:2676/1‑150 (MQ=255)
                                                                                                                        aaCGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga                   <  1:553395/150‑1 (MQ=255)
                                                                                                                        aaCGGGCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAAcgacga                   <  3:226365/150‑1 (MQ=255)
                                                                                                                         aCGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACgag                  >  1:1019355/1‑150 (MQ=255)
                                                                                                                          cGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACgaga                 >  3:655817/1‑150 (MQ=255)
                                                                                                                          cGGTCGCAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACgaga                 <  4:2676/150‑1 (MQ=255)
                                                                                                                           ggTCGCAAAGTTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAGCGACGAGAc                >  3:450371/1‑150 (MQ=255)
                                                                                                                                cAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATgg           >  1:1012422/1‑150 (MQ=255)
                                                                                                                                cAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATgg           >  1:900724/1‑150 (MQ=255)
                                                                                                                                cAAAGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATgg           <  4:414663/150‑1 (MQ=255)
                                                                                                                                 aaaGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGATCACGACGAGACCATGgc          >  1:881105/1‑150 (MQ=255)
                                                                                                                                 aaaGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACgaga                 <  3:502950/143‑1 (MQ=255)
                                                                                                                                 aaaGCTGAAAAAACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAACCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACgaga                 >  4:502950/1‑143 (MQ=255)
                                                                                                                                        aaaaaaCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATGGCGCTGTTg   >  4:458096/1‑150 (MQ=255)
                                                                                                                                         aaaaACATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATGGCGCTGTTGa  >  2:1004554/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
ACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGCCTGGATGATCTTCTGTTCGACTGGTCGAAGTGCCGGGTTAACGACGAGACCATGGCGCTGTTGA  >  NC_000913_3_bme_pgi/4233667‑4233953

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