Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,591,930 G→C E92Q (GAG→CAG)  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,9300GC95.2% 56.6 / ‑6.2 21E92Q (GAG→CAG) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (10/10);  minor base A (0/1);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATG  >  NC_000913/4591691‑4592158
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
ggctgGTTGTGGCCTTTTTTGGTCTTTTACAACCGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCTCTGTTTTACAATCCAATCGCCAGCAGCGATCCACTCTGAATGGCAGCTGGGGCGCGTTGTTGCAGTGGCGTAACACCCTCAACCGGGCGGGTTTCCGCTACATGAGAAATAATAATAATATTGGTAGCGGGTCAACCGTTGCTGAGCTGATGCAGAGtacaagt                                                                                                                                                                                                                           <  2:248928/250‑7 (MQ=255)
 gctgGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAAGAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGt                                                                                                                                                                                                                           <  2:190242/250‑1 (MQ=255)
          ggCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATAGTGGTAGCGGTTCAACCGTTGCTGAGCTGAGGCAGAGTGCCAGTATTTCGCTg                                                                                                                                                                                                                  <  2:26649/250‑1 (MQ=255)
                              aaCTGACAGCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAATAAAATTTCACTGTTTTACAAACCATGCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGTATCATAATAATATGTGTATCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTgg                                                                                                                                                                                             <  1:161866/251‑1 (MQ=255)
                                   aCATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGc                                                                                                                                                                                           >  2:9848/1‑248 (MQ=255)
                                         ggCGGTCTGTGCTTTAATGCCGTAAAGAATGACAAAGAAAATGTCACTGTTTTACAAACCAGTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGATGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCGGAGCTGATGCAGAGTGCCAGTATGTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGa                                                                                                                                                                                  <  2:72650/251‑1 (MQ=255)
                                                        aaTGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTg                                                                                                                                                                    >  1:121145/1‑250 (MQ=255)
                                                                                        gttgttccaacccATTCGCCAGCAGCAATCCCCGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATGGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATGACGAAGCGTTGCCTCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCagag                                                                                                                                    <  2:70799/239‑1 (MQ=255)
                                                                                                        cGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGTCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACg                                                                                                                    >  1:174832/1‑250 (MQ=255)
                                                                                                         gCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGa                                                                                                                   >  2:60450/1‑250 (MQ=255)
                                                                                                                    tCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGGCCCGCGTCAGAGAACCGCCGCAGCGGCAGAGATAAAACGCAATACCGATATTTATcaa                                                                                                         >  2:231897/1‑248 (MQ=255)
                                                                                                                            gAATGGCAGCGGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGTCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTg                                                                                                <  1:79590/250‑1 (MQ=255)
                                                                                                                                      tGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGCATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGCAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                       >  1:146777/1‑217 (MQ=255)
                                                                                                                                      tGGGTCGCGTTGTGGCAGACGCGTAACCCCCTCAACCGCGCGGGTATCCGCTACATGATGGATAAGAATAATATTGGTAGCGGTTCAACCGTGGCTGAGCTGATGCAGAGTGCCAGTATGTCGCTGAAACAGGCGGAAAAAAACGGGGCAGATTACGAGGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAAtt                                                                                                                       <  2:146777/217‑1 (MQ=255)
                                                                                                                                                 gttgCAGACGTGTAACACCCTAAACCGCGCGGGTATCCGCTACCTGATGGCTCAGAATAATATTGGCAGCGGTTCCACCGTTGCTGAGCTGCTGCAGCGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGAccc                                                                                                                                                                >  2:41312/1‑165 (MQ=255)
                                                                                                                                                 gttgCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCGGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGAccc                                                                                                                                                                <  1:41312/165‑1 (MQ=255)
                                                                                                                                                                                    tCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCCGGCAAAATCAACGAGTTCTTTGATCAGCCg                                       >  2:75145/1‑251 (MQ=255)
                                                                                                                                                                                      cGCTACATGATGGATCAGAATAATATTGGTAGCGGTGCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAc                                     <  1:9848/251‑1 (MQ=255)
                                                                                                                                                                                                      aGAATAATATTGGGAGCGGTCCAACCGTTGATGGGCAGAGTAAGAGTGCCAGTATTTCGCGGAAACAGGCGGAAAAAAACGGGGCAGATTACGAAGCGTGGCCGCGTGACCCGGGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGATCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGAACCAGGGATAGCAGGAc                     <  2:174832/251‑1 (MQ=255)
                                                                                                                                                                                                             tatTGGTAGCGGTTCAACCGTTGCTGAGCTGATGCAGCGTGCCAGTATTTCGCTGAAACAGGCGGACAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCg              >  1:236951/1‑251 (MQ=255)
                                                                                                                                                                                                                         ttCAACCGTTGCTGAGCTGATGCAGAGTGCCAGTATTTCGCTGAACCAGGCGGAAAAAAACTGGGCAGATTACGACGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCAAACTATTAGGTGCAGACAAAATCAACGAGTTTTTTGATCAGCCGACGCAGGGATATCAGGACGGTTTCGCGAAGCAGtctg  >  1:130615/1‑248 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
TGCTGGTTTTGGCCGTTTTTGGCCTTTTACAACTGACATCAGGCGGTCTGTTCTTTAATGCCTTAAAGAATGACAAAGAAAATTTCACTGTTTTACAAACCATTCGCCAGCAGCAATCCACGCTGAATGGCAGCTGGGTCGCGTTGTTGCAGACGCGTAACACCCTCAACCGCGCGGGTATCCGCTACATGATGGATCAGAATAATATTGGTAGCGGTTCAACCGTTGCTGAGCTGATGGAGAGTGCCAGTATTTCGCTGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATG  >  NC_000913/4591691‑4592158

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