Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,965,782 T→G E389D (GAA→GAC)  tar ← methyl‑accepting chemotaxis protein II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,965,7820TG81.8% 24.6 / 0.6 11E389D (GAA→GAC) tarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (4/5);  total (5/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

TTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAC  >  NZ_CP009273/1965702‑1965856
                                                                                |                                                                          
ttGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCa                                                                   >  1:224997/1‑90 (MQ=255)
       ctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCgg                                                            >  1:14224/1‑90 (MQ=255)
         tttGCCGCCGGGGCGCTGCGACTGGCAAGATTACGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCt                                                          <  1:12439/90‑1 (MQ=255)
                              cTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                     <  2:477205/90‑1 (MQ=255)
                                     gATTCCGCACTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGatat                              <  2:297455/90‑1 (MQ=37)
                                             aCTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTgg                      <  2:493665/90‑1 (MQ=255)
                                              cTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGa                     >  2:274013/1‑90 (MQ=255)
                                              cTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGa                     >  1:244817/1‑90 (MQ=255)
                                              cTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGa                     >  2:239009/1‑90 (MQ=255)
                                                    ccGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCaa               <  2:161613/90‑1 (MQ=255)
                                                                 aaaaCCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc  <  1:274013/90‑1 (MQ=255)
                                                                                |                                                                          
TTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAC  >  NZ_CP009273/1965702‑1965856

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑