Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,708 C→T intergenic (‑475/‑2199) 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,7080CT98.3% 184.5 / ‑4.3 59intergenic (‑475/‑2199)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (37/21);  total (37/22)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

CAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGC  >  NC_000913_3_bme_pgi/4233572‑4233853
                                                                                                                                        |                                                                                                                                                   
cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTc                                                                                                                                         <  2:865981/149‑1 (MQ=255)
 aaCGATGGGGTCATGCGGATTTTTCATCCGCTCCTGGCGGTCAGAAGTTCCGCTGATTAATGCTTCACTGTGCTAAGGGTTTACACACAACATTACGCTAACGACACCAAAACCATCACATTTTTATGTGACTGGTGCTACAATGTTcc                                                                                                                                        >  1:845396/1‑149 (MQ=255)
 aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGCCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCCACAATCTTcc                                                                                                                                        >  2:858098/1‑149 (MQ=255)
    gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCaaa                                                                                                                                     <  1:681660/149‑1 (MQ=255)
     aTGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAACTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAg                                                                                                                                    <  1:996369/149‑1 (MQ=255)
      tGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGt                                                                                                                                   >  2:557109/1‑149 (MQ=255)
           tCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACaa                                                                                                                              >  2:959535/1‑149 (MQ=255)
                cGGATTTTTCATCCACTCCTGGCGGACAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTCCAATCTTCCAGATTAACAATtcta                                                                                                                         >  2:302967/1‑148 (MQ=255)
                   aTTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                                                                                      <  1:488669/149‑1 (MQ=255)
                    tttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCaaaa                                                                                                                     >  2:46183/1‑149 (MQ=255)
                       ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                                                                                  >  2:428084/1‑149 (MQ=255)
                             cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACCTTTATTGAATACCTGTGTGACAATCTTCCAAGAACGCAATTCTCAAAATCAGAagag                                                                                                            >  1:976229/1‑149 (MQ=255)
                             cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCAACAATCTTCCAAAGTCACAAATCTCAAAAACAGAagag                                                                                                            >  1:309844/1‑149 (MQ=255)
                                   tGGCGGTCGGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATAAGAAGAGTATTGc                                                                                                      >  2:27825/1‑149 (MQ=255)
                                   tGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGc                                                                                                      >  2:658065/1‑149 (MQ=255)
                                        gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATg                                                                                               <  2:165802/149‑1 (MQ=255)
                                        atCAGTAGGTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATg                                                                                               <  1:221222/148‑1 (MQ=255)
                                             tagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAag                                                                                          >  1:630181/1‑149 (MQ=255)
                                             tagtTCAGCTAATAAATGCTGCACTGCGCTAAGGTTTTACACTCAACATTACGCTAACGGCACCACAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCATAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAag                                                                                          <  2:443229/149‑1 (MQ=255)
                                               gtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAagag                                                                                        >  1:702211/1‑149 (MQ=255)
                                                ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAgaga                                                                                       >  1:848708/1‑149 (MQ=255)
                                                    gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGAcgcg                                                                                   >  2:867417/1‑149 (MQ=255)
                                                        ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGCCGccacaa                                                                               >  2:996933/1‑147 (MQ=255)
                                                           aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGc                                                                            <  1:959535/149‑1 (MQ=255)
                                                                ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAGATCAGAAGAGTATTGTCA‑‑ATGGCAAGAGATGCGACGAAGCTTGaa                                                                       <  2:181782/149‑1 (MQ=255)
                                                                ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGaa                                                                       >  2:630179/1‑149 (MQ=255)
                                                                     tGCGCTAAGGGTTTACACTGAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAAc                                                                  >  2:319882/1‑149 (MQ=255)
                                                                             gggTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTTCTA‑‑ATGGCAAGATACGCGACGAAGCTTGAAGCAACGGTCGCaa                                                          >  2:491993/1‑149 (MQ=255)
                                                                              ggTTGACACTCAACATTACGTTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCaaa                                                         <  2:378441/149‑1 (MQ=255)
                                                                                tttACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAACGACACCAATCTCAACAT‑‑GAATGGCAGTACCAGCAGGGGATGAGAGGCGACGAAGCTTGACGCAACGGTCGCagatc                                                       >  1:444572/1‑145 (MQ=255)
                                                                                    cacTCAACATTACGCTATCGGCTCGAAGACCATCACATTTTTCTGTGACTGGTGGTTCAATCTTCCAAAGTCGCAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCTAAGCTGaa                                                   <  2:845396/149‑1 (MQ=255)
                                                                                     acTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                       >  1:198987/1‑144 (MQ=255)
                                                                                     acTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                                                       <  2:198987/144‑1 (MQ=255)
                                                                                      cttaaCATTACGCTAACGGCACCAAAACCATCACATGTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAACCCCAAAAAGCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaa                                                 <  2:418431/146‑1 (MQ=255)
                                                                                       tCAACATTTCGCTAACGGCACCAAACCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                <  1:179536/149‑1 (MQ=255)
                                                                                       tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGaaaaa                                                >  1:804889/1‑149 (MQ=255)
                                                                                       tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGATTGAACCAACGGTCGCAAAGCTGaaaaa                                                >  2:578824/1‑149 (MQ=255)
                                                                                       tCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTATCAAAATCAGAGGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGACGAAACGGTCGCAAAGATggcaga                                                >  2:23733/1‑144 (MQ=255)
                                                                                         aaCATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGc                                              >  2:448676/1‑149 (MQ=255)
                                                                                             ttACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAg                                                                             <  1:135550/114‑1 (MQ=255)
                                                                                             ttACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAg                                                                             >  2:135550/1‑114 (MQ=255)
                                                                                                 gCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCGAAGCTGAAAAAGCATTGGGCg                                      >  1:764003/1‑149 (MQ=255)
                                                                                                   tAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTACAAAGTCACAAGTCTAAAAATTAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGAAACGCTGAAAAAGCAttagcggga                                    >  2:166051/1‑142 (MQ=255)
                                                                                                      cGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATc                                 <  1:587199/149‑1 (MQ=255)
                                                                                                       ggCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCg                                >  1:617681/1‑149 (MQ=255)
                                                                                                          accaaaaCCATCACATTTTTTTTTTTTTTTTGCTACAATCTTCCAAAGTCACAATTCTAAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGATACGCGACGAAGCTTGAAGCAACGGTCTTATAGCTGTAAATTCATTGGGCGGAATCGGcc                             >  2:418960/4‑149 (MQ=255)
                                                                                                             aaaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACTAAGCACAAAGCAACGGTCGCAAGGCTGAAAAAGCATTGGGCGGAATCGGccccc                          >  1:693107/1‑148 (MQ=255)
                                                                                                              aaaCCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGACCCcgc                         >  2:211890/1‑149 (MQ=255)
                                                                                                                  cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                     <  1:106156/149‑1 (MQ=255)
                                                                                                                  cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                     >  1:180973/1‑149 (MQ=255)
                                                                                                                  cATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGCATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGa                     <  2:804889/149‑1 (MQ=255)
                                                                                                                   aTCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTGTTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGat                    >  1:868997/1‑149 (MQ=255)
                                                                                                                   aTCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGat                    <  2:17231/149‑1 (MQ=255)
                                                                                                                    tCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAATATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCCGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGata                   <  2:947209/149‑1 (MQ=255)
                                                                                                                           tttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATAAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATgcgcgc            <  2:54193/149‑1 (MQ=255)
                                                                                                                            ttCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGcgcgcg           <  1:46183/149‑1 (MQ=255)
                                                                                                                              cTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGc         >  2:33135/1‑149 (MQ=255)
                                                                                                                                     tGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc  >  2:574728/1‑149 (MQ=255)
                                                                                                                                     tGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCACAAAACTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAgc  >  1:264107/1‑149 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
CAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCGGCCCCGCGCGATATGCGCGCGGCATTCAGC  >  NC_000913_3_bme_pgi/4233572‑4233853

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: