Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,625,150 T→C R267R (AGA→AGG yhiJ ← YhiJ/YhiL family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,625,1500TC92.9% 39.7 / ‑3.3 14R267R (AGA→AGGyhiJYhiJ/YhiL family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (5/8);  total (5/9)
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

TGGTCGCCCAGAAACGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCA  >  NZ_CP009273/3625067‑3625234
                                                                                   |                                                                                    
tGGTCGCCCAGAAACGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTc                                                                                <  1:229068/90‑1 (MQ=255)
tGGTCGCCCAGAAACGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTc                                                                                <  2:365546/90‑1 (MQ=255)
 ggTCGCCCAGAAACGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCa                                                                               <  1:17234/90‑1 (MQ=255)
       ccAGAAACGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATAc                                                                         <  2:1726/90‑1 (MQ=255)
           aaaCGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATACTGat                                                                     <  1:59371/90‑1 (MQ=255)
               gTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATACTGattatt                                                                 <  2:153189/90‑1 (MQ=255)
                      gCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATACTGATTATTGTATTaa                                                          >  1:167032/1‑90 (MQ=255)
                      gCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATACTGATTATTGTATTaa                                                          >  1:303298/1‑90 (MQ=255)
                         tCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATACTGATTATTGTATTAACAt                                                       >  2:450333/1‑90 (MQ=255)
                                                  aTAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGtt                              <  1:31792/90‑1 (MQ=255)
                                                          gAAGACTTTCTAAAATAATGTCTCTCCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTa                      >  1:400295/1‑90 (MQ=255)
                                                                   cTAAAATAATGTCTCTCCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGt             <  2:181651/90‑1 (MQ=255)
                                                                           atGTCTCTCCTGTTCAGGATACTGATTATTGTATTa                                                           >  1:286692/1‑36 (MQ=38)
                                                                           atGTCTCTCCTGTTCAGGATACTGATTATTGTATTa                                                           <  2:286692/36‑1 (MQ=38)
                                                                              tctctcCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCa  <  2:279897/90‑1 (MQ=255)
                                                                              tctctcCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCa  <  2:303298/90‑1 (MQ=255)
                                                                                   |                                                                                    
TGGTCGCCCAGAAACGTACATGGCGTCACAGCGTCTTTATTAATAACTAAATAAGGCTGAAGACTTTCTAAAATAATGTCTCTTCTGTTCAGGATACTGATTATTGTATTAACATTAGTGTTAAGTTGCTCATTAAAGTTCTCCCTTAGCAATACGTTGAGAACTTCA  >  NZ_CP009273/3625067‑3625234

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