Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,746 A→G intergenic (‑513/‑2161) 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,7460AG100.0% 286.7 / NA 83intergenic (‑513/‑2161)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (52/31);  total (52/31)

CTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGC  >  NC_000913_3_bme_pgi/4233603‑4233889
                                                                                                                                               |                                                                                                                                                 
cTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTAt                                                                                                                                             <  2:25717/150‑1 (MQ=255)
      gCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTACAACCAACACAGTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCCCAATTCTCAAAATCAGAGGAGGATTGCTaa                                                                                                                                       >  2:854469/1‑150 (MQ=255)
                 ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGAc                                                                                                                            <  2:772322/150‑1 (MQ=255)
                       taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAGGTCACAATTCTCAAAATCAGAGGAGTATTGATAATGGCAAGAGACGcgacg                                                                                                                      >  1:555795/1‑150 (MQ=255)
                        aataaATGCTTCACTGCGCTAAGGGTTTACACCCAACATTACGCTAACGGCACTAAAACCATCACATATTTCCTTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGCAGTATCACAAATGGCAAGAGACGcgacga                                                                                                                     >  2:962685/1‑150 (MQ=255)
                            aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCt                                                                                                                 >  1:665822/1‑150 (MQ=255)
                              tGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTg                                                                                                               >  2:888246/1‑150 (MQ=255)
                               gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                                                                     <  2:1273394/127‑1 (MQ=255)
                               gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                                                                     >  1:1273394/1‑127 (MQ=255)
                                cTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                                                                             <  1:105608/150‑1 (MQ=255)
                                 ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAACCTTGaac                                                                                                            >  1:849515/1‑149 (MQ=255)
                                  tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGc                                                                                                           <  2:708670/150‑1 (MQ=255)
                                       gcgcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCACAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGGGACGAGGCTTGAAGCAACgg                                                                                                      >  1:371823/1‑150 (MQ=255)
                                           tAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGACCAGAGACGCGACGAAGCTTGAAGCAACGGTCGc                                                                                                  >  2:1041358/1‑150 (MQ=255)
                                            aaGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCa                                                                                                 >  1:111519/1‑150 (MQ=255)
                                             aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACGTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAACTCCGAGGAGTATTGCTAATGGCAAGAGACGGGACGAAGCTTGTAGCAAAGTTCGcca                                                                                                >  1:257516/1‑148 (MQ=255)
                                             aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCaa                                                                                                >  1:422112/1‑150 (MQ=255)
                                              gggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTATTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCACAATCAGAGGAGTATTGCTAATGTCAAGAGACGCGAATAAGCTTTGAGCAACGGTCGCaaa                                                                                               >  2:1162329/1‑150 (MQ=255)
                                               ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGGCGAAGGTTGCAGCAACGGTCGCAAAg                                                                                              >  1:1352376/1‑150 (MQ=255)
                                               ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                              >  2:1304425/1‑150 (MQ=255)
                                                gTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAAACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                             >  2:485994/1‑150 (MQ=255)
                                                 tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGAAAGAGAAGCGACGCAGCTTGAACCAACGGTTACAAAACt                                                                                            >  2:271813/1‑150 (MQ=255)
                                                 ttaacacTCAACATGAGTGTAACGGCACGAAAACCATCACATGTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                                                            <  2:216060/147‑1 (MQ=255)
                                                     cacTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaa                                                                                        <  2:1272599/150‑1 (MQ=255)
                                                      acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACCATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                                                       >  1:893378/1‑150 (MQ=255)
                                                      acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTCAAGCAACGGCCGCAAAGCTGaaaa                                                                                       <  1:1207316/150‑1 (MQ=255)
                                                      acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTAAAGCAGCGGTCGCAAAGCTGaaaa                                                                                       >  1:881527/1‑150 (MQ=255)
                                                       cTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                                                      >  1:781282/1‑150 (MQ=255)
                                                        tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCACGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGACAAAg                                                                                     >  1:102299/1‑150 (MQ=255)
                                                        tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                                                                     >  1:105031/1‑150 (MQ=255)
                                                        tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                                                                     >  1:739365/1‑150 (MQ=255)
                                                        tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAg                                                                                     >  2:1186019/1‑150 (MQ=255)
                                                          aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGACAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCa                                                                                   >  1:646001/1‑150 (MQ=255)
                                                             aatACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAACCTGAAGAAGCATTg                                                                                >  2:241252/3‑150 (MQ=255)
                                                             aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                             >  2:1238140/1‑137 (MQ=255)
                                                             aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                                                             <  1:1238140/137‑1 (MQ=255)
                                                             aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTGCAAGCAACGGACGCAAAGCTGAAAAAGCATTg                                                                                >  1:1122472/1‑150 (MQ=255)
                                                             aTTAAGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg                                                                                >  2:241229/1‑150 (MQ=255)
                                                              ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg‑g                                                                              >  1:787472/1‑150 (MQ=255)
                                                              ttACGCTAACGGCAATAAAACCACCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGAGGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg‑g                                                                              <  2:970508/150‑1 (MQ=255)
                                                               tACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                                                                    >  2:989927/1‑112 (MQ=255)
                                                               tACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                                                                    <  1:989927/112‑1 (MQ=255)
                                                               tACGCTAACGGCACTAAAACCATCACAGTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGGATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTg‑gg                                                                             <  2:221570/150‑1 (MQ=255)
                                                                aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGc                                                                            <  1:68194/150‑1 (MQ=255)
                                                                        ggCACTAAAACCATCACATTTTTCTGTGCCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCgg                                                                    <  1:285108/150‑1 (MQ=255)
                                                                        ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCgg                                                                    <  1:1221057/150‑1 (MQ=255)
                                                                        ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCgg                                                                    <  2:815803/150‑1 (MQ=255)
                                                                         gCACTAAAACCACCACATTTTTCTTTGTGGGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGc                                                                   <  1:1360578/150‑1 (MQ=255)
                                                                                aaCGATCACATTTTTCTGGGACTGGGGCTACAATCTTCCAAAGTCAAAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCcgcgc                                                            <  2:182128/150‑1 (MQ=255)
                                                                                 aCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCcgcgcg                                                           >  2:684780/1‑150 (MQ=255)
                                                                                       acaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCTACGGTCGCAAAGCTGAAAGAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATgg                                                    >  1:421492/1‑149 (MQ=255)
                                                                                          tttttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATgcgcg                                                 >  1:852503/1‑150 (MQ=255)
                                                                                              tCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCGCATAT‑CGCGCggg                                             <  1:2610/150‑2 (MQ=255)
                                                                                                   gACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGACAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCa                                        <  2:1379280/150‑1 (MQ=255)
                                                                                                     cTGGCGCATCAATCTTCCAAAGGCACAATCCTCAAAAGCAGAGCAGTATTGCTAATGGCAAAAGACGCGACGAAGCTTGAAGCAACCGTCCCAAAGCTCAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCT‑ATATCCGCGCGGCATTCAgc                                      <  2:609455/150‑1 (MQ=255)
                                                                                                      tGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAgcg                                     >  2:1259500/1‑150 (MQ=255)
                                                                                                       ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAgcgc                                    <  2:781282/150‑1 (MQ=255)
                                                                                                        gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGcc                                   >  1:61779/1‑150 (MQ=255)
                                                                                                          gcTACAATCTTCCAAAGTCACAATTCCCAAAACCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGa                                 <  1:281349/150‑1 (MQ=255)
                                                                                                              cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGAcccc                             >  2:1127180/1‑150 (MQ=255)
                                                                                                                aTCTTCCAAAGTCACAATTCGCCAAATCAGAGGAGGATTGCTAAGGGCAAGAGACGCGACGAAGCTTGAAGCCACGGTCGCAAAGCTGAAACAGCATTG‑GGCAGAATCGGCCCCGCGCG‑ATATACGCACGCCAGTCAGCGCCGCACccgg                           >  2:616026/1‑150 (MQ=255)
                                                                                                                   ttCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCccggccg                        >  2:1190602/1‑150 (MQ=255)
                                                                                                                      cAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGtt                     >  1:692255/1‑150 (MQ=255)
                                                                                                                       aaaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGTGCGGCAGTCAGCGCCGACCCCGGCCGGTTc                    <  2:1004464/150‑1 (MQ=255)
                                                                                                                       aaaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCACTCTCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                    >  1:960225/1‑150 (MQ=255)
                                                                                                                          gTCACAATTCCCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 <  1:296111/150‑1 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAGGCACCGGTCGCAAAGCGGAAAACCCATTGAGGCCAAATCAG‑GCCTCGCT‑AGAGGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                >  2:296111/1‑150 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 <  2:559355/149‑1 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 >  1:559355/1‑149 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                >  2:1211474/1‑150 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                >  1:1325630/1‑150 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCCGTTCgggg                >  1:202826/1‑149 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 <  2:541631/147‑1 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 >  1:541631/1‑147 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCgt              >  2:1328279/1‑149 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa              >  2:978238/1‑150 (MQ=255)
                                                                                                                                 ttCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGcc                                   >  2:1280710/1‑125 (MQ=255)
                                                                                                                                 ttCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGcc                                   <  1:1280710/125‑1 (MQ=255)
                                                                                                                                   ctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCAAAGAAGCTTGCAGCCACGGTCGCAAAGCTGCAAAAGGATCG‑GGCGGAAGCGGCCCATCGCG‑GTATGCGGGCGGCAGGAAGCGCCGACCCCGGCCGGTTCGGGGCATATTcc        >  2:475154/1‑149 (MQ=255)
                                                                                                                                    tcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc       >  1:70861/1‑150 (MQ=255)
                                                                                                                                      aaaaGCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc     <  2:1301150/150‑1 (MQ=255)
                                                                                                                                        aaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTg   <  1:1186019/150‑1 (MQ=255)
                                                                                                                                         aTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc  <  1:1390412/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTG‑GGCGGAATCGGCCCCGCGCG‑ATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGC  >  NC_000913_3_bme_pgi/4233603‑4233889

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