Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,065,426 T→G T129P (ACC→CCC)  rutF ← pyrimidine utilization flavin reductase protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,065,4260TG93.8% 44.4 / ‑3.1 17T129P (ACC→CCC) rutFpyrimidine utilization flavin reductase protein F
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (8/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTG  >  NZ_CP009273/1065341‑1065510
                                                                                     |                                                                                    
tGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCg                                                                                  >  1:274613/1‑90 (MQ=255)
  cgcaCTGCGATCGACCCACACCAGCCCGTAGGTTGTGGGGGGACGATGAACCGCTTCCATGGCGCGAAACACAATGTCGTGGGGGCCGAc                                                                                <  1:319977/87‑1 (MQ=255)
                              tAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGa                                                                                 <  1:366009/61‑1 (MQ=255)
                              tAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGa                                                                                 >  2:366009/1‑61 (MQ=255)
                                   tgtggtgtgACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGa                                               >  2:171704/1‑90 (MQ=255)
                                         gtgACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCag                                         >  2:535545/1‑90 (MQ=255)
                                                 gAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCtt                                 >  2:116421/1‑90 (MQ=255)
                                                           gATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCgg                       <  1:375353/90‑1 (MQ=255)
                                                                cgcAGAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATc                  <  2:313441/90‑1 (MQ=255)
                                                                cgcAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATc                  <  2:319498/90‑1 (MQ=255)
                                                                cgcAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATc                  <  1:583640/90‑1 (MQ=255)
                                                                         gAATGTCGTGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCAcccc        >  1:79237/1‑90 (MQ=255)
                                                                              tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGaa                                              >  1:459849/1‑48 (MQ=255)
                                                                              tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGaa                                              <  2:459849/48‑1 (MQ=255)
                                                                              tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCAcccc        <  2:190345/86‑1 (MQ=255)
                                                                              tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCAcccc        >  1:190345/1‑86 (MQ=255)
                                                                                gTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTg  >  2:428333/1‑90 (MQ=255)
                                                                                     |                                                                                    
TGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTG  >  NZ_CP009273/1065341‑1065510

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