Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,745 Δ2 bp intergenic (‑512/‑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,7450A.100.0% 294.8 / NA 61intergenic (‑512/‑2162)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (28/33);  total (28/33)
*NC_000913_3_bme_pgi4,233,7460A.100.0% 294.8 / NA 61intergenic (‑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 . (28/33);  total (28/33)

TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCT  >  NC_000913_3_bme_pgi/4233599‑4233885
                                                                                                                                                  ||                                                                                                                                           
tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCaggaga                                                                                                                                           >  4:14242/1‑150 (MQ=255)
     tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGc                                                                                                                                    <  4:362085/150‑1 (MQ=255)
      ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCt                                                                                                                                   <  3:180312/150‑1 (MQ=255)
       cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTa                                                                                                                                  <  2:354448/150‑1 (MQ=255)
            ggTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGc                                                                                                                             <  3:261409/150‑1 (MQ=255)
                  tagtTCAGCTAATAAATGCTTCCCTGCGCTAAGGGTTTACCCTCAGCATTACGTTAACGGCACTAAAACCAGCACAGTTTTCTGTGAGTGGCGCTACAATTATACGAAGTCACAATCCTCAAAATCAG‑‑GAGTATTGCTGATGGCAAgaga                                                                                                                       <  4:210334/150‑1 (MQ=255)
                   agtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGAc                                                                                                                      <  1:170655/150‑1 (MQ=255)
                     ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGAcgc                                                                                                                    <  1:200912/150‑1 (MQ=255)
                     ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGAcgc                                                                                                                    <  2:33037/150‑1 (MQ=255)
                      tCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGAcgcg                                                                                                                   >  2:474055/1‑150 (MQ=255)
                       cAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGcga                                                                                                                  >  1:350804/1‑150 (MQ=255)
                       cAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGcga                                                                                                                  >  4:165603/1‑150 (MQ=255)
                         gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGcgacg                                                                                                                >  1:474290/1‑150 (MQ=255)
                           taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGaa                                                                                                              >  1:63903/1‑150 (MQ=255)
                            aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAg                                                                                                             <  1:474055/150‑1 (MQ=255)
                            aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAg                                                                                                             <  2:350804/150‑1 (MQ=255)
                             ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGc                                                                                                            <  4:318101/150‑1 (MQ=255)
                               aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCtt                                                                                                          >  4:122077/1‑150 (MQ=255)
                               aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCtt                                                                                                          <  1:156791/150‑1 (MQ=255)
                                  tGCTTCACCGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGaa                                                                                                       <  2:369879/150‑1 (MQ=255)
                                      tCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCaa                                                                                                   >  1:212969/1‑150 (MQ=255)
                                                 aGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAg                                                                                        >  1:454673/1‑150 (MQ=255)
                                                   ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCt                                                                                      >  2:258778/1‑150 (MQ=255)
                                                     tttACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGa                                                                                    >  1:576897/1‑150 (MQ=255)
                                                           cTCAACATTACGCTAACGGCACTAAAACCATCACGTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                                                              <  4:216385/150‑1 (MQ=255)
                                                              aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                                                           >  1:360715/1‑150 (MQ=255)
                                                              aaCATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCAtt                                                                           >  3:29275/1‑150 (MQ=255)
                                                                 aTTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTggg                                                                        >  4:272174/1‑150 (MQ=255)
                                                                     cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGa                                                                    >  3:211888/1‑150 (MQ=255)
                                                                     cGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGa                                                                    >  1:329160/1‑150 (MQ=255)
                                                                      gCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGaa                                                                   >  1:470599/1‑150 (MQ=255)
                                                                      gCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGaa                                                                   <  3:122077/150‑1 (MQ=255)
                                                                        tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCaa                                                                                                   >  2:245844/1‑116 (MQ=255)
                                                                        tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCaa                                                                                                   <  1:245844/116‑1 (MQ=255)
                                                                        tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                 >  4:227838/1‑150 (MQ=255)
                                                                        tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                 <  2:503629/150‑1 (MQ=255)
                                                                        tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                 >  2:73604/1‑150 (MQ=255)
                                                                        tAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATCGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                                                 <  4:242328/150‑1 (MQ=255)
                                                                                  aaaaCCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCcgcgc                                                       <  2:184267/150‑1 (MQ=255)
                                                                                          cacaTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcg                                               >  2:225634/1‑150 (MQ=255)
                                                                                              ttttttAGTGCCGGGCGCGAAAATCTGCGAAAGGGAAAATACTCAAAATCAG‑‑GAGTATTGCAAATCGCAAGAGACGCGAGGAAGCTTTAAGCAACGGAAAAAAAAAAAAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg                                           <  4:144621/150‑1 (MQ=255)
                                                                                              tttttCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg                                           <  1:109227/150‑1 (MQ=255)
                                                                                                    tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTGTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTc                                     >  4:183704/1‑150 (MQ=255)
                                                                                                    tgtgACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTc                                     >  2:318417/1‑150 (MQ=255)
                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                 <  2:371986/150‑1 (MQ=255)
                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                 <  2:372007/150‑1 (MQ=255)
                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgcg                                 <  2:371994/150‑1 (MQ=255)
                                                                                                          tGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGGATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGcc                               <  4:57508/150‑1 (MQ=255)
                                                                                                           ggCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCg                              >  2:158163/1‑150 (MQ=255)
                                                                                                            gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGa                             <  4:312731/150‑1 (MQ=255)
                                                                                                             cgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAc                            >  3:3827/1‑150 (MQ=255)
                                                                                                             cgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAc                            <  1:67651/150‑1 (MQ=255)
                                                                                                               cTACAATCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGAccc                          <  2:394650/150‑1 (MQ=255)
                                                                                                                  cAATCGGCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccgg                       >  2:573386/1‑150 (MQ=255)
                                                                                                                   aaTCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggc                      >  1:340836/1‑150 (MQ=255)
                                                                                                                     tCTTCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCccggccg                    <  1:170857/150‑1 (MQ=255)
                                                                                                                       ttCCAAAGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGt                  >  2:552771/1‑150 (MQ=255)
                                                                                                                           aaaGTCACAATCCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCga              <  4:313014/150‑2 (MQ=255)
                                                                                                                            aaGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCggg             <  1:535818/150‑1 (MQ=255)
                                                                                                                             aGTCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGc            <  1:73604/150‑1 (MQ=255)
                                                                                                                               tCACAATTCTCAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGa          >  1:79801/1‑150 (MQ=255)
                                                                                                                                       ctcAAAATCAG‑‑GAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGct  <  3:219999/150‑1 (MQ=255)
                                                                                                                                                  ||                                                                                                                                           
TCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGCGCCGACCCCGGCCGGTTCGGGCGATATTCGCT  >  NC_000913_3_bme_pgi/4233599‑4233885

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