Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,572,102 T→C intergenic (‑57/+305) gadB ← / ← pqqL glutamate decarboxylase B/putative zinc peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,572,1020TC100.0% 60.6 / NA 18intergenic (‑57/+305)gadB/pqqLglutamate decarboxylase B/putative zinc peptidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (2/16);  total (2/16)

CGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATA  >  NC_000913/1572025‑1572174
                                                                             |                                                                        
cGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCaa                                                                        <  1:263028/80‑1 (MQ=255)
               aTCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAgg                                                           <  1:263030/78‑1 (MQ=255)
               aTCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAgg                                                           <  2:263029/78‑1 (MQ=255)
                       tAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGttt                                                 <  1:263031/80‑1 (MQ=255)
                         aaCTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTaa                                               <  1:263033/80‑1 (MQ=255)
                         aaCTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTaa                                               <  2:263032/80‑1 (MQ=255)
                                tAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACaa                                        >  1:263053/1‑80 (MQ=255)
                                    aTGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATcc                                   <  1:263034/81‑1 (MQ=255)
                                       aTTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTaa                                <  1:263035/81‑1 (MQ=255)
                                        ttGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAAt                               <  1:263036/81‑1 (MQ=255)
                                           gATCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGtt                            <  2:263037/81‑1 (MQ=255)
                                             tCGCATTAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGttat                          <  1:263038/81‑1 (MQ=255)
                                                  ttAAAAAAGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATc                     <  1:263039/81‑1 (MQ=255)
                                                     aaaaaGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGt                   <  1:263040/80‑1 (MQ=255)
                                                     aaaaaGTAGGATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGt                   <  2:263041/80‑1 (MQ=255)
                                                              gATTTATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAc         >  2:263055/1‑81 (MQ=255)
                                                                  tATCGATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTAcgcg      <  1:263042/80‑1 (MQ=255)
                                                                      gATAAAGCAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGaata  <  1:263043/80‑1 (MQ=255)
                                                                             |                                                                        
CGTTACTTGCTTCTTATCCATTTTAAACTCCTTAAAATGATTGGATCGCATTAAAAAAGTAGGATTTATCGATAAAGTAAGCAAGTGCAAAGGACTCGTGTTTAAATAACAAAATCCTAATGTTATTTATCGTGAGATATTACGCGAATA  >  NC_000913/1572025‑1572174

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