Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,967,130 G→A A103V (GCC→GTC)  cheY ← chemotaxis protein CheY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,967,1300GA84.6% 20.9 / ‑7.0 13A103V (GCC→GTC) cheYchemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/6);  minor base C (1/1);  total (6/7)
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

GGCTTCCGCCGCTTCGGCAATGGCCTGATCCAGCCCCAGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGG‑CTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATAT‑TGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTAC  >  NC_000913/1966897‑1967373
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
ggCTTCCGCCTCTTCGGCAATGGCCTGATCCAGCCCCAGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCCCCAGCGTCGCGGCGGTAAATGG‑CGTCACCACATAGCCACTGACCCCCGCTTGcgccgc                                                                                                                                                                                                                                      <  1:144038/250‑1 (MQ=255)
                                                                              gATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGG‑CTTCACCACATAGCCACTGACCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAAtn                                                                                                                                                                                 <  1:107397/225‑2 (MQ=255)
                                                                                               cGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGG‑CTTCACCACATAGCCACTGACCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTACCACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTcccg                                                                                                                                      >  1:22/1‑249 (MQ=255)
                                                                                                 ccAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGG‑CTTCACCACATAGCCACTGACCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAggn                                                                                                                                    <  1:49816/251‑2 (MQ=255)
                                                                                                                               gatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGAGGTAAATGG‑CTTCACCACATAGCCACTGAccccc                                                                                                                                                                                                                                                 >  1:64105/1‑112 (MQ=255)
                                                                                                                                                       tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGG‑CTTCACCACATAGCCACTGACCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAg                                                                                                                                               <  1:36470/186‑1 (MQ=255)
                                                                                                                                                       tCACATGCCCAGTTTCTCAAAGATTTNGTTGAGNNTTNCCTCCAGCGTCGCGGCGGTNAANNN‑CNNNANAACATNGCNNCTGACCCCCNCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCCCCATTAACNCTGGCAATGNNGACATCGCGCCATCCGCCCGANTTGNNTNNanggancg                                                                                                                                         >  2:36470/1‑185 (MQ=255)
                                                                                                                                                                          aaGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGCAAATGG‑CTTCACCCCCTACACACAGCCACCCGCTTGCGCAGCAGCAAAGATGTTATCTTTCTTCGCTTCTGCCGTCAACATTAACAATGGCCATGCCGACATCGCGACATTCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATAT‑TGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGACTGCAACTTATTgaggccttt                                                          >  1:174242/1‑245 (MQ=255)
                                                                                                                                                                            nanatngttgANTTTTNCCTCCAGCGTCGCGGNGGTAAATGG‑CTTCACCACANNGNCACGGNccccc                                                                                                                                                                                                                                                 <  2:64105/63‑1 (MQ=255)
                                                                                                                                                                                             cctccAGCGTCGCGGCGGTAAATGG‑CTTCCCCACATAGCCACTGACCCCCGCTTGCGCCGCAGCAATTATGTTCTCTTTCTTCGCTTCTGCATTCCCCATTAACACTGGCAATGCCGACATCGCCCCACCCGCACGAATTGTTTTCAGCAATTCGAGGCCATCCATAG‑TGGGCAGTTTCCAGTCGGAGATAAAAAAGCCATACCCGCCTGCCTGCAACGTATTGAGAGCGTCGACGCCATCTTCCGCTTc                                        <  1:68045/250‑1 (MQ=255)
                                                                                                                                                                                                   gCGTCGCGGCGG‑AAATGGGCTGCCCCACATACCCCCCGCCCCCCGCTGGCGCCGCAGCAATGGTGTTCTCTTTCACCCCTTCTGCAGACACCATAAACAGTGGCATGGCCACCATCGCGCCATCCGCCCAAATTTGTTACAGCAATCCCAGCCCAGCCA‑ATGTGGGCGTTTTCCTGCCGGAGAGAACAAATCCTTAACCGCCTGCCTGCAACTTATTGAGAGc                                                             <  1:169152/223‑1 (MQ=255)
                                                                                                                                                                                                           gcggTNNATNN‑NTTNNCCACATAGCCACTGACCCCCGCTNGCGCCGCAGCAATGATGTTCTCTTTCTTCNCNTNTGCANNCACNATTAACACTGNCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGNCNATCCATAT‑TGGGCATGTTCCAGTCGgag                                                                                                      >  2:127889/1‑174 (MQ=255)
                                                                                                                                                                                                           gcggTAAATGG‑CTTCACCACATAGCCACTGACCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTGCTTCGCCTCTGCAGTCACCATTAACACGGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATAT‑TGGGCATGTTCCAGTCGgag                                                                                                      <  1:127889/174‑1 (MQ=255)
                                                                                                                                                                                                                                   nccACTGACCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATAT‑TGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTAc  >  1:35446/2‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
GGCTTCCGCCGCTTCGGCAATGGCCTGATCCAGCCCCAGTTCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGG‑CTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATAT‑TGGGCATGTTCCAGTCGGAGATAACAAATCCATAACCGCCTGCCTGCAACTTATTGAGAGCGTCGACGCCATCTTCCGCTTCCTCAACATTATTGAATCCCAGCTCTTTCAGCAGGTTAC  >  NC_000913/1966897‑1967373

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: ‑