Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,233,680 T→C intergenic (‑447/‑2227) 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,6800TC98.3% 190.8 / ‑5.2 60intergenic (‑447/‑2227)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (34/25);  minor base G (0/1);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCG  >  NC_000913_3_bme_pgi/4233536‑4233823
                                                                                                                                                |                                                                                                                                                 
gtaCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAGCAGTACGCTAAGGGCGCCaaaa                                                                                                                                               >  1:443229/1‑149 (MQ=255)
  aCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAAcc                                                                                                                                             <  1:297447/149‑1 (MQ=255)
   cATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCa                                                                                                                                            <  1:897660/149‑1 (MQ=255)
    aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                           <  2:613736/149‑1 (MQ=255)
    aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                           <  2:169444/149‑1 (MQ=255)
    aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCACCGTGAATTATGATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCAt                                                                                                                                           <  1:534605/149‑1 (MQ=255)
     tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATc                                                                                                                                          <  1:50402/149‑1 (MQ=255)
       tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCCCCAAAACCATcac                                                                                                                                        >  1:54193/1‑149 (MQ=255)
       tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATccc                                                                                                                                        >  1:17231/1‑147 (MQ=255)
       tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATcac                                                                                                                                        <  1:525026/149‑1 (MQ=255)
                 gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCtgt                                                                                                                              <  2:351856/149‑1 (MQ=255)
                 gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAACAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCtgt                                                                                                                              <  2:283652/149‑1 (MQ=255)
                         ggggCGGTTGTCAACGATGGGGTCATGCGGAGTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCACAACCATCACATTTGTCTTTGACgggtg                                                                                                                      >  2:318410/1‑147 (MQ=255)
                                    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)
                                                              aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATc                                                                                                                                          <  2:539683/92‑1 (MQ=255)
                                                              aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACCAAAACCATc                                                                                                                                          >  1:539683/1‑92 (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                         <  2:198987/144‑1 (MQ=255)
                                                                                                                         acTCAACATTACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGc                         >  1:198987/1‑144 (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                                               >  2:135550/1‑114 (MQ=255)
                                                                                                                                 ttACGCTAACGGCACCAAAACCATCACATTTTTCTGTGACTGGTGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAg                                               <  1:135550/114‑1 (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)
                                                                                                                                                |                                                                                                                                                 
GTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTA‑‑ATGGCAAGAGACGCGACGAAGCTTGAAGCAACGGTCGCAAAGCTGAAAAAGCATTGGGCGGAATCG  >  NC_000913_3_bme_pgi/4233536‑4233823

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: