Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,591,677 T→C I7I (ATT→ATC tsr → methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,591,6770TC86.7% 28.2 / ‑3.4 15I7I (ATT→ATCtsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (4/9);  total (5/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGT  >  NC_000913/4591486‑4591902
                                                                                                                                                                                               |                                                                                                                                                                                                                                 
tAAATATGTGATTCCTATCACATATTTCTATTGTGAATAATTTATGTAAATAATGACTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAACGAGAAACCATGTTACAACGTATCAAAATTGTGACCAGCTTActg                                                                                                                                                                                                                    >  1:41130/1‑207 (MQ=255)
gaaaTATGTGATGCATATATAATATTTATAGTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATGCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAAATTGTGACCAGCTTActg                                                                                                                                                                                                                    <  2:41130/206‑1 (MQ=255)
                                                                                       ccatggatggcgggttacaccttttCATAAAGTTTTTGCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAg                                                                                 <  1:144058/228‑1 (MQ=255)
                                                                                               ggcgggttaccccttttCATAAAGTTTTTGCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTcg                                                                         <  2:173209/236‑1 (MQ=255)
                                                                                                             tttCATAAAGTTTTTGCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGAc                                                            <  1:483538/249‑1 (MQ=255)
                                                                                                                 aTAAAGATTTTGCTTTCAAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCATTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAACGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTAGCGTTGTTGCAGACGCGt                                                        >  2:19553/1‑250 (MQ=255)
                                                                                                                  aaaaGTTTTGGCTTTCCAGGCCGAAAAACTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTGCGTTAATGCCTTAAAGAATGTCAAAGAAAATTGCACTGTTTTACAAACCTTTCGCCAGCAGCAATCCTCGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTaa                                                      <  2:427648/250‑1 (MQ=255)
                                                                                                                      gTTTTTGCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACAcc                                                  <  1:393773/251‑1 (MQ=255)
                                                                                                                            gctttccAGGCCGAAAATCTTGCATCGGTCCACCGGCAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGCCGCGTAACCCCCTCAAc                                            >  2:453242/2‑251 (MQ=255)
                                                                                                                            gctttccAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTGTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCaa                                             <  2:148581/249‑1 (MQ=255)
                                                                                                                                     gCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGt                                   <  2:423884/251‑1 (MQ=255)
                                                                                                                                                tGCATCGGTCCACATGAAAGAGAAACCATGTTAAAACGTATCAAGATCGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTa                           <  1:250633/248‑1 (MQ=255)
                                                                                                                                                      ggTCCACAGGACAGAGAGACCAGGTTAAAAAGTATAAAGCTCGTGACCAGCTTACTGCTGATTTTGTCCGTTTTTGGCCTTTTACAACTGACATCAGGCGTTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATATCACTGTTTTACAAACCATTAGCCAGCAGCAATCCACGATGAATGGCCGCTGGATCGCGTTGTTGAAGACGCGTCACCCCCTCAAACGCGCAGGTATcccacacatggttgg                  >  2:303381/1‑238 (MQ=255)
                                                                                                                                                            caGGAAAGCGAAACCATGTTAAAACGTATCAAGATCGTGACAAGCTTAATGGTGGCTATGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCATTGTTTTACAACACATTCGCCAGCAGCAATCCACGCTGAATGGCAGTTAGGTCGCGTTGTTGCCGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGa             >  2:248287/1‑250 (MQ=255)
                                                                                                                                                                      aaaCCATTTTAAAACGTATCAAGATCGTGACCAGTTTACTGCTGGTTTTTGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGAAGACGCGTAACACTCTCAACCGCGCGGGTTTCCGCTACATGATGGATCAGAATAATATTGGt  <  2:360983/251‑1 (MQ=255)
                                                                                                                                                                                               |                                                                                                                                                                                                                                 
TAAATATGTGATTCATATCACATATTTATATTGTGAATAATTTATGTAAATAATGGCTTTTTAAATTCAGAGTGTGAATAAAATTCACTCGGCGTAATCTCCGCGGGATATTCATAAAGTTTTTCCTTTCCAGGCCGAAAATCTTGCATCGGTCCACAGGAAAGAGAAACCATGTTAAAACGTATCAAAATTGTGACCAGCTTACTGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGT  >  NC_000913/4591486‑4591902

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