Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,423,096 A→G V249A (GTT→GCT)  qorB ← NAD(P)H:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,423,0960AG86.7% 33.9 / ‑2.9 15V249A (GTT→GCT) qorBNAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (10/3);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGAT  >  NZ_CP009273/4423014‑4423180
                                                                                  |                                                                                    
ttCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGa                                                                               >  1:134014/1‑90 (MQ=255)
ttCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGa                                                                               >  1:184821/1‑90 (MQ=255)
                     gTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAg                                                          <  1:63262/90‑1 (MQ=255)
                       ggCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTc                                                        >  1:263817/1‑90 (MQ=255)
                       ggCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTc                                                        >  2:291977/1‑90 (MQ=255)
                        gCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTcc                                                       >  2:365302/1‑90 (MQ=255)
                               aatTTGCTAAGTTTTTTGCTGTCATCAAACAGCCCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGc                                                <  2:62459/90‑1 (MQ=255)
                                atTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCa                                               <  2:134014/90‑1 (MQ=255)
                                         gtgtTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGc                                      <  1:365302/90‑1 (MQ=255)
                                                   gTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGc                            >  1:164700/1‑90 (MQ=255)
                                                   gTCATCAAACAGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGc                            >  1:265837/1‑90 (MQ=255)
                                                             aGACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGaaa                  >  2:2592/1‑90 (MQ=255)
                                                              gACCGCCTTTCGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAAt                 >  1:435525/1‑90 (MQ=255)
                                                                        cGATGCGCCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCt       >  1:416167/1‑90 (MQ=255)
                                                                             cgcCAGCGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGAt  >  2:243332/1‑90 (MQ=255)
                                                                                  |                                                                                    
TTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCCAACGTCAGAATCCGCCAGCATATCCGCCAGTCCGTCGGGCAGTCCGACGCTTTTCAGTGCCGCGGCGAAATCGGCTTCGCTCAGAT  >  NZ_CP009273/4423014‑4423180

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