Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,022,885 G→A *90* (TAG→TAA fliQ → flagellar biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,022,8850GA100.0% 80.9 / NA 26*90* (TAG→TAAfliQflagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/13);  total (13/13)

CGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGG  >  NC_000913/2022773‑2023016
                                                                                                                |                                                                                                                                   
cGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCg                                                                                                           >  2:308415/1‑139 (MQ=255)
  aaaaTCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGaa                                                                                                         <  1:391852/139‑1 (MQ=255)
                                        aCCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                                        >  2:206398/1‑102 (MQ=255)
                                        aCCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                                        <  1:206398/102‑1 (MQ=255)
                                                  cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATgg                                                                                                   >  1:1137352/1‑97 (MQ=255)
                                                  cTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATgg                                                                                                   <  2:1137352/97‑1 (MQ=255)
                                                    cAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTa                                                                  <  2:620375/128‑1 (MQ=255)
                                                    cAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTa                                                                  >  1:620375/1‑128 (MQ=255)
                                                                 gATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                              <  1:620920/87‑1 (MQ=255)
                                                                 gATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                              >  2:620920/1‑87 (MQ=255)
                                                                       gTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACcgcg                                    <  1:1158724/139‑1 (MQ=255)
                                                                           gCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                                   >  1:830316/1‑120 (MQ=255)
                                                                           gCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                                   <  2:830316/120‑1 (MQ=255)
                                                                                   ttCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGc                        >  1:1108102/1‑139 (MQ=255)
                                                                                         aaCCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTa                                                                                       >  1:220288/1‑70 (MQ=255)
                                                                                         aaCCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTa                                                                                       <  2:220288/70‑1 (MQ=255)
                                                                                         aaCCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGc                  >  1:513603/1‑139 (MQ=255)
                                                                                           ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTa                                                                                       >  2:1259453/1‑68 (MQ=255)
                                                                                           ccTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTa                                                                                       <  1:1259453/68‑1 (MQ=255)
                                                                                               ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                                        >  1:441407/1‑47 (MQ=255)
                                                                                               ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                                        <  2:441407/47‑1 (MQ=255)
                                                                                               ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGc                                                        <  2:985793/95‑1 (MQ=255)
                                                                                               ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGc                                                        >  1:985793/1‑95 (MQ=255)
                                                                                               ccGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTa            <  1:308415/139‑1 (MQ=255)
                                                                                                         tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg  >  2:595274/1‑139 (MQ=255)
                                                                                                         tcGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACggg  <  2:513603/139‑1 (MQ=255)
                                                                                                                |                                                                                                                                   
CGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTACCGAAACGGG  >  NC_000913/2022773‑2023016

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