Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,510,238 T→G intergenic (+445/+452) insO → / ← fecE IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,510,2380TG100.0% 63.5 / NA 20intergenic (+445/+452)insO/fecEIS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (13/7);  total (13/7)

GCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATT  >  NC_000913/4509997‑4510474
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
ncgcTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCAGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCttt                                                                                                                                                                                                                                     >  1:41984/2‑251 (MQ=255)
                                  ggCCCGTCTCCCTTCACTTTATGCCAGGTTCGGTAAGCTTGGTGAGACACCGACTCATGAAGATATCATTGATAAAAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTTAGGAAGCTTCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAAATGGATTAGCAATTCCGTTCCCGGGTTGCTTTTGGTTTTCGTTTGCGTCTGATCCTTGATTTAttc                                                                                                                                                                                                   >  1:135854/1‑249 (MQ=255)
                                    nccGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTc                                                                                                                                                                                                 >  1:113461/2‑251 (MQ=255)
                                                    ttATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCtn                                                                                                                                                                                   <  1:274510/249‑2 (MQ=255)
                                                      aTGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTg                                                                                                                                                                               <  1:384199/251‑1 (MQ=255)
                                                                              agaCACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGAAAACGCACACTTTTTAGAAAAATAGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGt                                                                                                                                                       >  1:198733/1‑251 (MQ=255)
                                                                                acacCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAtt                                                                                                                                                                                                        >  1:231738/1‑200 (MQ=255)
                                                                                                                 cGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCg                                                                                                                                                                                                                      >  1:24810/1‑153 (MQ=255)
                                                                                                                  gATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGaaa                                                                                                                                  >  1:410991/1‑236 (MQ=255)
                                                                                                                     ccATCAACTGGCATGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACagag                                                                                                                                       >  1:157943/1‑228 (MQ=255)
                                                                                                                            nnggNNNNNATANCGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTg                                                                                                                                                                                                                  <  2:267846/144‑1 (MQ=255)
                                                                                                                                                gCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTACTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCGGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGtt                                                                                     <  1:443710/251‑1 (MQ=255)
                                                                                                                                                  cGGTGATTAAGAAAGCTGCACCTACCTAAGTAGCAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCa                                                                                                                                                                    >  1:388630/1‑170 (MQ=255)
                                                                                                                                                   ggTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGttttt                                                                                  >  1:140821/1‑251 (MQ=255)
                                                                                                                                                   ggTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTGTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAGTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCGGATCATTGATTGATCCACAAAATCCCAACCTCATTGGTAATGAACCAGCTCCGTTAATGTCCGCTCTGGCACagag                                                                                                                                       <  1:250524/198‑1 (MQ=255)
                                                                                                                                                            gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAGCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTccc                                                                                                                    <  1:189572/208‑1 (MQ=255)
                                                                                                                                                              aaGCTGCACCGACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTATCAATTCCGTTCCAGGGTTGCTTTTGATTTa                                                                                                                                                                                                                              <  1:96861/100‑1 (MQ=255)
                                                                                                                                                                                 agtagCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAACCTGGATTAGCACTTACGCTCCAGGGATGCTTTTGATTTACGTAGGCGTCTGATCATTGATTTATCCTAAAAAGCCCAACCTCATTGGTAATGACCCAGCTCCGTGAATGTc                                                                                                                                                      >  1:257035/1‑153 (MQ=255)
                                                                                                                                                                                                                                  aGCAANTCCGTNCNNGGGCCGCTTTAGATTTACGNTGGCGTCTGATCATTNATTTNTCCTCnacngn                                                                                                                                                                                           >  2:175313/1‑63 (MQ=255)
                                                                                                                                                                                                                                   ncAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGACATGTCCGTTAAGCAGAGTTCAAAatt  >  1:265019/2‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
GCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTTCAAAATT  >  NC_000913/4509997‑4510474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: