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,7160TC96.2% 77.5 / ‑3.6 26I300T (ATC→ACC) tilStRNA lysidine(34) synthetase TilS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (16/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CCGCGACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACG  >  NZ_CP009273/209630‑209792
                                                                                      |                                                                            
ccGCGACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCg                                                                           <  2:487550/90‑1 (MQ=255)
       gcgTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTAt                                                                    <  2:631279/90‑1 (MQ=255)
                  gATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACt                                                         >  1:668384/1‑90 (MQ=255)
                      tGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtgg                                                     >  1:139955/1‑90 (MQ=255)
                       ggCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTgtggt                                                    <  2:668384/90‑1 (MQ=255)
                                                 gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGc                          >  1:405767/1‑90 (MQ=255)
                                                 gATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGc                          >  1:662830/1‑90 (MQ=255)
                                                    gCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                       <  1:332033/90‑1 (MQ=255)
                                                       tCACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAAc                    >  2:57867/1‑90 (MQ=255)
                                                        cACCCTGTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACa                   >  2:93560/1‑90 (MQ=255)
                                                              gTTTACGTTTGGGTGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                       >  2:163961/1‑80 (MQ=255)
                                                              gTTTACGTTTGGGTGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaa                       <  1:163961/80‑1 (MQ=255)
                                                              gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCa                                                                  >  2:11526/1‑37 (MQ=38)
                                                              gTTTACGTTTGGGCGCGTTTGAAACCCGACGCTATCa                                                                  <  1:11526/37‑1 (MQ=38)
                                                              gTTTACGCTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGc             >  1:598527/1‑90 (MQ=255)
                                                                    gTTTGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGc       <  2:373650/90‑1 (MQ=255)
                                                                      ttGGGCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAg     >  1:207986/1‑90 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaaa                      >  2:277523/1‑71 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGaaa                      <  1:277523/71‑1 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGc       <  2:174475/86‑1 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGc       >  1:174475/1‑86 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc   >  1:457228/1‑90 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc   >  2:312200/1‑90 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc   >  2:349059/1‑90 (MQ=255)
                                                                        gggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGAc   >  2:681012/1‑90 (MQ=255)
                                                                         ggCGCGTTTGAAACCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACg  <  1:559871/90‑1 (MQ=255)
                                                                                      |                                                                            
CCGCGACGCGTTGGTGAGGATCTGGCAGGAAGTGGCGCTGGCGCGGGAAGATGCCTCACCCTGTTTACGTTTGGGCGCGTTTGAAATCCGACGCTATCAGTCGCAACTGTGGTGGATTAAATCCGTCACCGGGCAAAGCGAAAACATTGTGCCGTGGCAGACG  >  NZ_CP009273/209630‑209792

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