Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,901,765 G→A G152G (GGC→GGT)  bglH ← carbohydrate‑specific outer membrane porin, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,901,7650GA100.0% 134.9 / NA 40G152G (GGC→GGT) bglHcarbohydrate‑specific outer membrane porin, cryptic
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/19);  total (21/19)

AACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGA  >  NC_000913/3901631‑3901896
                                                                                                                                      |                                                                                                                                   
aaCCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCa                                                                                                                                 >  2:821069/1‑139 (MQ=255)
              gTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGcc                                                                                                                   <  2:536645/139‑1 (MQ=255)
                    ctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGccc                                                                                                                  >  2:855295/1‑134 (MQ=255)
                    ctgACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGccc                                                                                                                  <  1:855295/134‑1 (MQ=255)
                       aCCAACGTTACCATCCATCATCACGACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTa                                                                                                          >  1:619736/1‑139 (MQ=255)
                         cAACGTTACCATCCATTATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGATCCAATCGCCCATGATTTAgg                                                                                                        <  2:619736/139‑1 (MQ=255)
                                    tccatcATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATgag                                                                                             >  1:339860/1‑139 (MQ=255)
                                      catcatCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTATACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTTTTCGTTACCAAAGCGACCCAGAGAACCAATCGCTCATGATTTAGGTGAACCATGaaag                                                                                           >  1:683903/1‑139 (MQ=255)
                                      catcatCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCagag                                                                                                                              >  1:325572/1‑104 (MQ=255)
                                      catcatCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCagag                                                                                                                              <  2:325572/104‑1 (MQ=255)
                                              caACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTccc                                                                                   <  1:821069/139‑1 (MQ=255)
                                                    cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                             >  2:472658/1‑139 (MQ=255)
                                                    cATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAAcc                                                                             >  2:188754/1‑139 (MQ=255)
                                                      tCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAg                                                                           >  1:232860/1‑139 (MQ=255)
                                                       cAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGa                                                                          <  1:188754/139‑1 (MQ=255)
                                                        aaCCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         <  1:472658/139‑1 (MQ=255)
                                                                tGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATa                                                                 <  1:352773/139‑1 (MQ=255)
                                                                                     gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                        <  1:324683/95‑1 (MQ=255)
                                                                                     gttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCg                                                                                        >  2:324683/1‑95 (MQ=255)
                                                                                      ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATgag                                                                                             >  1:406218/1‑89 (MQ=255)
                                                                                      ttgttTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATgag                                                                                             <  2:406218/89‑1 (MQ=255)
                                                                                         ttttAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCa                                        >  2:516236/1‑139 (MQ=255)
                                                                                           ttAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAg                                                                                                         >  1:243217/1‑72 (MQ=255)
                                                                                           ttAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAg                                                                                                         <  2:243217/72‑1 (MQ=255)
                                                                                                   aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                       <  2:98056/50‑1 (MQ=255)
                                                                                                   aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAAt                                                                                                                       >  1:98056/1‑50 (MQ=255)
                                                                                                   aaaTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAAGAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACaa                              >  2:572206/1‑139 (MQ=255)
                                                                                                     aTCAAACCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAAt                            >  1:243292/1‑139 (MQ=255)
                                                                                                          aCCAGCCGGAGTATTCGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACt                       >  1:431337/1‑139 (MQ=255)
                                                                                                                     tATTCGTTACCAAAGCGACCCAGAGTACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCACCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCa            <  2:683903/139‑1 (MQ=255)
                                                                                                                         cGTTACCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGt        <  2:431337/139‑1 (MQ=255)
                                                                                                                             aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGtt                                                           >  1:869801/1‑84 (MQ=255)
                                                                                                                             aCCAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGtt                                                           <  2:869801/84‑1 (MQ=255)
                                                                                                                               cAAAGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGga  <  2:339860/139‑1 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATa                                                                 >  2:668562/1‑74 (MQ=255)
                                                                                                                                 aaGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATa                                                                 <  1:668562/74‑1 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                    <  2:745191/38‑1 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGaa                                                                                                    >  1:745191/1‑38 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         <  1:864579/65‑1 (MQ=255)
                                                                                                                                  aGCGACCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGaa                                                                         >  2:864579/1‑65 (MQ=255)
                                                                                                                                      |                                                                                                                                   
AACCAGCCGGTAGAGTACTGCTGACCAACGTTACCATCCATCATCACAACGGCATCAACCCGTTTGCCGTTTTCGTTGTAGACACGTTGTTTTAACTGCAAATCAAACCAGCCGGAGTATTCGTTACCAAAGCGGCCCAGAGAACCAATCGCCCATGATTTAGGTGAACCATGAGAGGCGGTCCCCCAACCAGAACGGTAATAGCCGTTGTAACTAAAGCCAATTTCATCTTTCACAAATTTACTGAAATCTTTCAGAGTCATGGA  >  NC_000913/3901631‑3901896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑