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,7820TG100.0% 53.6 / NA 17E389D (GAA→GACtarmethyl‑accepting chemotaxis protein II
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (8/9);  total (8/9)

AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGAT  >  NZ_CP009273/1965696‑1965862
                                                                                      |                                                                                
aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                               >  1:135681/1‑90 (MQ=255)
aGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCa                                                                               >  2:39077/1‑90 (MQ=255)
  ggCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCAcc                                                                             >  1:319612/1‑90 (MQ=255)
   gCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCAccc                                                                            >  1:50238/1‑90 (MQ=255)
          tctcTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCCCTTCCCCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcg                                                                     <  2:5050/90‑1 (MQ=25)
           ctctTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCAcgcgc                                                                    <  2:50238/90‑1 (MQ=255)
                  gccgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCa                                                             >  2:263474/1‑90 (MQ=255)
                    cgccTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAc                                                           <  1:92505/90‑1 (MQ=255)
                                    cTGGCAAGATTCCGCCCTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                           <  1:184824/90‑1 (MQ=255)
                                    cTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGcgcg                                           <  1:69795/90‑1 (MQ=255)
                                     tGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGa                                          <  1:402263/90‑1 (MQ=255)
                                                     ttCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAc                                                           <  1:32234/57‑1 (MQ=255)
                                                     ttCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAAc                                                           >  2:32234/1‑57 (MQ=255)
                                                      tCACCCGCCACCACGGCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAg                         >  2:58387/1‑90 (MQ=255)
                                                                     gCAAAACCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATa          <  1:165073/90‑1 (MQ=255)
                                                                       aaaaCCACGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAAc        <  1:386942/90‑1 (MQ=255)
                                                                             aCGGCCCTGGTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGAt  >  1:318699/1‑90 (MQ=255)
                                                                                      |                                                                                
AGGGCTTTGATCTCTTTTGCCGCCTGGGCGCTGCGACTGGCAAGATTACGCACTTCACCCGCCACCACGGCAAAACCACGGCCCTGTTCACCCGCACGCGCGGCTTCAACCGCGGCATTCAGCGCGAGGATATTAGTCTGGAAGGCAATACCGTCGATAACGCTGAT  >  NZ_CP009273/1965696‑1965862

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