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)

CTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGC  >  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)
                                                              ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTgg                                                                              >  1:787472/1‑150 (MQ=255)
                                                              ttACGCTAACGGCAATAAAACCACCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGAGGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTgg                                                                              <  2:970508/150‑1 (MQ=255)
                                                               tACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                                                                    >  2:989927/1‑112 (MQ=255)
                                                               tACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                                                                    <  1:989927/112‑1 (MQ=255)
                                                               tACGCTAACGGCACTAAAACCATCACAGTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGGATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                             <  2:221570/150‑1 (MQ=255)
                                                                aCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGc                                                                            <  1:68194/150‑1 (MQ=255)
                                                                        ggCACTAAAACCATCACATTTTTCTGTGCCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                    <  1:285108/150‑1 (MQ=255)
                                                                        ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                    <  1:1221057/150‑1 (MQ=255)
                                                                        ggCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCgg                                                                    <  2:815803/150‑1 (MQ=255)
                                                                         gCACTAAAACCACCACATTTTTCTTTGTGGGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGc                                                                   <  1:1360578/150‑1 (MQ=255)
                                                                                aaCGATCACATTTTTCTGGGACTGGGGCTACAATCTTCCAAAGTCAAAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                                                            <  2:182128/150‑1 (MQ=255)
                                                                                 aCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgcg                                                           >  2:684780/1‑150 (MQ=255)
                                                                                       acaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCTACGGTCGCAAAGCTGAAAGAGCATTGGGCGGAATCGGCCCCGCGCGATATgg                                                    >  1:421492/1‑149 (MQ=255)
                                                                                          tttttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcg                                                 >  1:852503/1‑150 (MQ=255)
                                                                                              tCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGCATATCGCGCggg                                             <  1:2610/150‑2 (MQ=255)
                                                                                                   gACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGACAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCa                                        <  2:1379280/150‑1 (MQ=255)
                                                                                                     cTGGCGCATCAATCTTCCAAAGGCACAATCCTCAAAAGCAGAGCAGTATTGCTAATGGCAAAAGACGCGACGAAGCTTGAAGCAACCGTCCCAAAGCTCAAAAAGCATTGGGCGGAATCGGCCCCGCGCTATATCCGCGCGGCATTCAgc                                      <  2:609455/150‑1 (MQ=255)
                                                                                                      tGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                     >  2:1259500/1‑150 (MQ=255)
                                                                                                       ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcgc                                    <  2:781282/150‑1 (MQ=255)
                                                                                                        gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                   >  1:61779/1‑150 (MQ=255)
                                                                                                          gcTACAATCTTCCAAAGTCACAATTCCCAAAACCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGa                                 <  1:281349/150‑1 (MQ=255)
                                                                                                              cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAcccc                             >  2:1127180/1‑150 (MQ=255)
                                                                                                                aTCTTCCAAAGTCACAATTCGCCAAATCAGAGGAGGATTGCTAAGGGCAAGAGACGCGACGAAGCTTGAAGCCACGGTCGCAAAGCTGAAACAGCATTGGGCAGAATCGGCCCCGCGCGATATACGCACGCCAGTCAGCGCCGCACccgg                           >  2:616026/1‑150 (MQ=255)
                                                                                                                   ttCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccg                        >  2:1190602/1‑150 (MQ=255)
                                                                                                                      cAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGtt                     >  1:692255/1‑150 (MQ=255)
                                                                                                                       aaaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGTGCGGCAGTCAGCGCCGACCCCGGCCGGTTc                    <  2:1004464/150‑1 (MQ=255)
                                                                                                                       aaaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCACTCTCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTc                    >  1:960225/1‑150 (MQ=255)
                                                                                                                          gTCACAATTCCCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 <  1:296111/150‑1 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAGGCACCGGTCGCAAAGCGGAAAACCCATTGAGGCCAAATCAGGCCTCGCTAGAGGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                >  2:296111/1‑150 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 <  2:559355/149‑1 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 >  1:559355/1‑149 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                >  2:1211474/1‑150 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc                >  1:1325630/1‑150 (MQ=255)
                                                                                                                           tCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCCGTTCgggg                >  1:202826/1‑149 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 <  2:541631/147‑1 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg                 >  1:541631/1‑147 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCgt              >  2:1328279/1‑149 (MQ=255)
                                                                                                                             acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa              >  2:978238/1‑150 (MQ=255)
                                                                                                                                 ttCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                   >  2:1280710/1‑125 (MQ=255)
                                                                                                                                 ttCTCAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                                   <  1:1280710/125‑1 (MQ=255)
                                                                                                                                   ctcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCAAAGAAGCTTGCAGCCACGGTCGCAAAGCTGCAAAAGGATCGGGCGGAAGCGGCCCATCGCGGTATGCGGGCGGCAGGAAGCGCCGACCCCGGCCGGTTCGGGGCATATTcc        >  2:475154/1‑149 (MQ=255)
                                                                                                                                    tcAAAATCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGc       >  1:70861/1‑150 (MQ=255)
                                                                                                                                      aaaaGCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGctc     <  2:1301150/150‑1 (MQ=255)
                                                                                                                                        aaTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTg   <  1:1186019/150‑1 (MQ=255)
                                                                                                                                         aTCAGAGGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGc  <  1:1390412/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCTCTGC  >  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: