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,8850GA97.3% 115.0 / NA 37*90* (TAG→TAAfliQflagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (17/19);  minor base C (1/0);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTA  >  NC_000913/2022754‑2023006
                                                                                                                                   |                                                                                                                         
aTGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACt                                                                                                                    <  2:391029/139‑1 (MQ=255)
                   cGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCg                                                                                                 >  1:798711/1‑139 (MQ=255)
                           atcGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATgg                                                                                         <  2:798711/139‑1 (MQ=255)
                                   aTTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                              >  1:362703/1‑126 (MQ=255)
                                   aTTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                              <  2:362703/126‑1 (MQ=255)
                                      tATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGtt                                                                              <  2:512913/139‑1 (MQ=255)
                                       aTCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGACTATCCTGGTTa                                                                             >  2:866535/1‑139 (MQ=255)
                                            cATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTAt                                                                                                                  >  2:769854/1‑97 (MQ=255)
                                            cATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTAt                                                                                                                  <  1:769854/97‑1 (MQ=255)
                                                         ggACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTa                                                                                                                      >  2:694023/1‑80 (MQ=255)
                                                         ggACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTa                                                                                                                      <  1:694023/80‑1 (MQ=255)
                                                         ggACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCg                                                           >  1:181539/1‑139 (MQ=255)
                                                            ccGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTa                                                        <  2:461770/139‑1 (MQ=255)
                                                                           cTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGACTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGgcg                                         <  1:866535/139‑1 (MQ=255)
                                                                             gTTGCTGGATTACGTCCGCATCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                              >  2:633014/1‑84 (MQ=255)
                                                                             gTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAAc                                                                                              <  1:633014/84‑1 (MQ=255)
                                                                               tGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGt                                                                               >  2:116238/1‑97 (MQ=255)
                                                                               tGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGt                                                                               <  1:116238/97‑1 (MQ=255)
                                                                                  tGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGa                                                                                                       <  2:46218/70‑1 (MQ=255)
                                                                                  tGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGa                                                                                                       >  1:46218/1‑70 (MQ=255)
                                                                                      ttACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                    <  1:181077/85‑1 (MQ=255)
                                                                                      ttACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                    >  2:181077/1‑85 (MQ=255)
                                                                                              gCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGa                      >  1:1029394/1‑139 (MQ=255)
                                                                                                      ttCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTaac                                                                                    <  2:645455/69‑3 (MQ=255)
                                                                                                      ttCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATc                                                                                    >  1:645455/1‑69 (MQ=255)
                                                                                                      ttCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACcgcg                          <  2:119973/127‑1 (MQ=255)
                                                                                                      ttCACTAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACcgcg                          >  1:119973/1‑127 (MQ=255)
                                                                                                           tAACCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACg         <  2:1029394/139‑1 (MQ=255)
                                                                                                            aaCCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCt                                                                >  2:242488/1‑83 (MQ=255)
                                                                                                            aaCCTGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCt                                                                <  1:242488/83‑1 (MQ=255)
                                                                                                                tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACg         <  2:30776/134‑1 (MQ=255)
                                                                                                                tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACg         >  1:30776/1‑134 (MQ=255)
                                                                                                                tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCCATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGa                <  2:283095/127‑1 (MQ=255)
                                                                                                                tGCCGTATATCATCGGGTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCCATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGa                >  1:283095/1‑127 (MQ=255)
                                                                                                                  ccGTATATCATCGGGTACCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTa  >  2:614708/1‑139 (MQ=255)
                                                                                                                               ggTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTaa                                                                            >  2:358390/1‑52 (MQ=255)
                                                                                                                               ggTAACCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTaa                                                                            <  1:358390/52‑1 (MQ=255)
                                                                                                                                   |                                                                                                                         
ATGACGCTGTCGTTTATTCCGAAAATCATCGCCGTATTTATCGCCATTATTATTGCCGGACCGTGGATGCTCAATCTGTTGCTGGATTACGTCCGCACCTTGTTCACTAACCTGCCGTATATCATCGGGTAGCCGTACTATGTTGCAGGTGACAAGCGAACAATGGCTATCCTGGTTAAACCTGTACTTCTGGCCGTTACTGCGCGTGCTGGCGCTGATCTCCACCGCGCCGATTCTGAGCGAACGCAGCGTA  >  NC_000913/2022754‑2023006

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