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% 81.4 / NA 20intergenic (‑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 . (12/8);  total (12/8)

AGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATAT  >  NC_000913/1978188‑1978342
                                                                               |                                                                           
aGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcc                                                                            <  1:313609/81‑1 (MQ=255)
 ggTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGccc                                                                           <  1:313610/81‑1 (MQ=255)
     tGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGcccggt                                                                        <  2:313611/80‑4 (MQ=255)
       cATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTaaa                                                                    <  2:313612/81‑1 (MQ=255)
        attattCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTaaa                                                                    >  1:313635/1‑80 (MQ=255)
         ttattCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTaaaaa                                                                  <  1:313613/81‑1 (MQ=255)
                 acccaGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAg                                                          <  1:313614/81‑1 (MQ=255)
                       aaTAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCAt                                                       >  2:313632/1‑78 (MQ=255)
                        aTAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGa                                                   <  2:313615/81‑1 (MQ=255)
                         tAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGaa                                                  >  2:313637/1‑81 (MQ=255)
                          aaCCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGaa                                                  >  1:313644/1‑80 (MQ=255)
                            ccAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAg                                               <  2:313616/81‑1 (MQ=255)
                                        ttATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCt                                   >  1:313650/1‑81 (MQ=255)
                                           tGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGc                                 >  1:313633/1‑80 (MQ=255)
                                             cgctttCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAAt                              <  2:313617/78‑1 (MQ=255)
                                                  tCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAg                           >  2:313646/1‑79 (MQ=255)
                                                  tCACCGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAg                           >  1:313634/1‑79 (MQ=255)
                                                     ccGCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGaa                      >  1:313640/1‑81 (MQ=255)
                                                       gCACCCCGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAac                     >  2:313645/1‑80 (MQ=255)
                                                             cGTGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCttt              <  1:313618/81‑1 (MQ=255)
                                                              atGATGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCtttt             <  1:313619/80‑1 (MQ=255)
                                                                 aTGTCGCCGGGAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTgg          >  1:313643/1‑81 (MQ=255)
                                                                         gggAAGCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAatat  >  1:313642/1‑81 (MQ=255)
                                                                               |                                                                           
AGGTATGCATTATTCCCACCCAGAATAACCAACTTTATTTTTATGCGGTTTCACCGCACCCCGTGATGTCGCCGGGAAGCCCCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATAT  >  NC_000913/1978188‑1978342

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: