Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,392,826 T→C Q18Q (CAA→CAG uspE ← universal stress protein UspE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,392,8260TC94.4% 49.8 / ‑3.4 18Q18Q (CAA→CAGuspEuniversal stress protein UspE
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (4/13);  total (4/14)
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.58e-01

TAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCC  >  NZ_CP009273/1392755‑1392890
                                                                       |                                                                
tAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGa                                                <  2:190548/90‑1 (MQ=255)
tAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGa                                                <  2:10650/90‑1 (MQ=255)
     cGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGAt                                           <  1:454306/90‑1 (MQ=255)
     cGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGAt                                           >  1:584967/1‑90 (MQ=255)
      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                          >  1:242487/1‑90 (MQ=255)
      ggCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                          >  2:156963/1‑90 (MQ=255)
      ggCAAAAAGGCTCTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATa                                          <  1:65147/90‑1 (MQ=255)
           aaaGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATaacaac                                     <  2:42938/90‑1 (MQ=255)
                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                               <  1:611623/90‑1 (MQ=255)
                 tttAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCAt                               <  1:253475/90‑1 (MQ=255)
                   tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                             <  1:706024/90‑1 (MQ=255)
                   tAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGt                             <  2:280664/90‑1 (MQ=255)
                     aTTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTc                           <  2:316952/90‑1 (MQ=255)
                             accaATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACAt                   <  2:483727/90‑1 (MQ=255)
                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                 <  1:297567/62‑1 (MQ=255)
                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGc                                 >  2:297567/1‑62 (MQ=255)
                                           gTAAATAAACAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGtc     <  2:584967/90‑1 (MQ=255)
                                              aataaaCAGCTCGCCGCAATGCTGGCTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTcc  <  1:100323/90‑1 (MQ=255)
                                                                       |                                                                
TAGATCGGCAAAAAGGCTTTAATTTTGCCACCAATCCGTTGATGTAAATAAACAGCTCGCCGCAATGCTGGTTGGTCGTCCTGGTTAGGATCGATAACAACGAGCATGTTCTGATACATAGCCATACAGGGTCTCC  >  NZ_CP009273/1392755‑1392890

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