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,4260TG91.3% 69.5 / ‑0.5 23T129P (ACC→CCC) rutFpyrimidine utilization flavin reductase protein F
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (12/9);  total (13/10)
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.98e-01

GATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTG  >  NZ_CP009273/1065342‑1065510
                                                                                    |                                                                                    
gATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGa                                                                                 <  2:427313/90‑1 (MQ=255)
 aTAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGAc                                                                                <  2:375210/90‑1 (MQ=255)
 aTAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGAc                                                                                <  2:109313/90‑1 (MQ=255)
    aCTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCt                                                                             >  1:1001395/1‑90 (MQ=255)
     cTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTc                                                                            <  2:878326/90‑1 (MQ=255)
         gATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTcacca                                                                        >  1:603429/1‑90 (MQ=255)
         gATCGAACCACACCAGCCCGTAGGGCGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcacca                                                                        >  1:791016/1‑90 (MQ=255)
              aaCCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTgg                                                                   <  2:920718/90‑1 (MQ=255)
                            gTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTc                                                     >  2:432071/1‑90 (MQ=255)
                                   gtggtgtgACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGaa                                              >  1:106691/1‑90 (MQ=255)
                                   gtggtgtgACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGaa                                              >  1:402683/1‑90 (MQ=255)
                                          gACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCagag                                       >  1:961/1‑90 (MQ=255)
                                           aCGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGc                                      >  1:169576/1‑90 (MQ=255)
                                                gAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCtt                                 >  2:93485/1‑90 (MQ=255)
                                                     gCTTCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGt                            >  1:147283/1‑90 (MQ=255)
                                                        tCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGc                         <  2:391542/90‑1 (MQ=255)
                                                        tCGATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGc                         <  2:603429/90‑1 (MQ=255)
                                                          gATGGCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCgg                       <  1:929374/90‑1 (MQ=255)
                                                             ggCGCAAAACAGAATGTCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCa                    <  1:210636/90‑1 (MQ=255)
                                                                             tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGaa                                              >  1:91809/1‑48 (MQ=255)
                                                                             tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGaa                                              <  2:91809/48‑1 (MQ=255)
                                                                             tCGTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTc    >  2:452898/1‑90 (MQ=255)
                                                                               gTGGGGGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTg  >  1:726090/1‑90 (MQ=255)
                                                                                    |                                                                                    
GATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTG  >  NZ_CP009273/1065342‑1065510

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