Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,978,270 Δ1 bp intergenic (‑73/+248) flhD ← / ← insB‑5 DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,978,2670C.100.0% 58.5 / NA 15intergenic (‑70/+251)flhD/insB‑5DNA‑binding transcriptional dual regulator FlhD/IS1 protein InsB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (4/11);  total (4/11)

GTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTT  >  NC_000913/1978190‑1978337
                                                                             |                                                                      
gTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcccg                                                                     <  2:352542/81‑2 (MQ=255)
 tATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcccgg                                                                    <  1:352543/81‑3 (MQ=255)
   tGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcccggt                                                                   <  2:352544/80‑4 (MQ=255)
    gCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcccggta                                                                  >  2:352564/1‑76 (MQ=255)
       ttattCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTaaaa                                                              >  1:352565/1‑80 (MQ=255)
          ttCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAaat                                                           <  1:352545/80‑1 (MQ=255)
                cccaGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGc                                                    <  1:352546/81‑1 (MQ=255)
                cccaGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGc                                                    <  1:352547/81‑1 (MQ=255)
                 ccaGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCa                                                   <  1:352548/81‑1 (MQ=255)
                          ccAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAg                                          <  2:352549/81‑1 (MQ=255)
                            aaCTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGTCCGGTAAAAAATAATTAGCATTAGAATAGtt                                        >  2:352568/1‑81 (MQ=255)
                              cTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGc                                      <  1:352550/81‑1 (MQ=255)
                                 tATTTTTATGCGGTTTCACCGAACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCg                                     <  1:352551/79‑1 (MQ=255)
                                       tATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCt                              <  2:352552/80‑1 (MQ=255)
                                            ggTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATa                        <  1:352553/81‑1 (MQ=255)
                                                  aCCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAg                   <  1:352554/80‑1 (MQ=255)
                                                    cGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAac                >  1:352567/1‑81 (MQ=255)
                                                            gTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCtttt        >  1:352572/1‑81 (MQ=255)
                                                                  tCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGttt  <  1:352555/81‑1 (MQ=255)
                                                                             |                                                                      
GTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTT  >  NC_000913/1978190‑1978337

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