Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 209,716 T→C I300T (ATC→ACC)  tilS → tRNA lysidine(34) synthetase TilS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273209,7160TC85.7% 43.1 / 2.7 21I300T (ATC→ACC) tilStRNA lysidine(34) synthetase TilS
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (11/7);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCA  >  NZ_CP009273/209647‑209800
                                                                     |                                                                                    
ggATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAAc                                                                  <  2:147282/90‑1 (MQ=255)
ggATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTCACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAAc                                                                  <  1:378663/90‑1 (MQ=255)
     tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  2:726784/1‑90 (MQ=255)
     tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  1:239354/1‑90 (MQ=255)
     tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  1:154067/1‑90 (MQ=255)
     tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                             >  1:492194/1‑90 (MQ=255)
          ggAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGAt                                                        >  2:125564/1‑90 (MQ=255)
           gAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGAtt                                                       >  2:542873/1‑90 (MQ=255)
                ggcgctggcgcGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATc                                                  >  1:292266/1‑90 (MQ=255)
                            ggAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCa                                      <  2:541779/90‑1 (MQ=255)
                              aaGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCaaa                                    <  2:239354/90‑1 (MQ=255)
                              aaGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCaaa                                    <  2:492194/90‑1 (MQ=255)
                                   gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                               <  1:189373/90‑1 (MQ=255)
                                   gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                               <  1:392538/90‑1 (MQ=255)
                                             gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGc                     >  1:35266/1‑90 (MQ=255)
                                             gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGc                     >  2:403763/1‑90 (MQ=255)
                                                tACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGt                  >  2:442051/1‑90 (MQ=255)
                                                       gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                               <  1:20115/70‑1 (MQ=255)
                                                       gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                               >  2:20115/1‑70 (MQ=255)
                                                       gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc           >  1:452750/1‑90 (MQ=255)
                                                                tGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCa  <  2:583242/90‑1 (MQ=255)
                                                                     |                                                                                    
GGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACGTGGCTTCA  >  NZ_CP009273/209647‑209800

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