Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,433,164 C→T *287* (TAG→TAA qorB ← NAD(P)H:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,433,1640CT95.5% 58.8 / ‑3.7 22*287* (TAG→TAAqorBNAD(P)H:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (12/9);  total (13/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 = 8.09e-01

GGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAACTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCC  >  NC_000913/4433033‑4433277
                                                                                                                                   |                                                                                                                 
ggCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttat                                                                                                            <  1:500064/139‑1 (MQ=255)
  cGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGttatta                                                                                                          >  1:707651/1‑139 (MQ=255)
           cgccgcTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAACTAGTTATTAACATTaaaa                                                                                                 >  1:805478/1‑139 (MQ=255)
              cgcTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTaa                                                                                                         <  2:522127/128‑1 (MQ=255)
              cgcTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTaa                                                                                                         >  1:522127/1‑128 (MQ=255)
                                  cATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTaa                                                                          >  2:332505/1‑139 (MQ=255)
                                                      gTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACg                                                                        <  2:638990/121‑1 (MQ=255)
                                                      gTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACg                                                                        >  1:638990/1‑121 (MQ=255)
                                                                 aGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCt                                                                                   >  1:423016/1‑99 (MQ=255)
                                                                 aGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCt                                                                                   <  2:423016/99‑1 (MQ=255)
                                                                             tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTaa                                                                                                         >  1:50000/1‑65 (MQ=255)
                                                                             tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTaa                                                                                                         <  2:50000/65‑1 (MQ=255)
                                                                             tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGt                                                                    <  1:597625/102‑1 (MQ=255)
                                                                             tCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGt                                                                    >  2:597625/1‑102 (MQ=255)
                                                                                  gCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGtcatc                          <  1:697510/139‑1 (MQ=255)
                                                                                                    agagGATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGa        >  2:609444/1‑139 (MQ=255)
                                                                                                       ggATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAg                                                                                               <  1:596987/49‑1 (MQ=255)
                                                                                                       ggATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAg                                                                                               >  2:596987/1‑49 (MQ=255)
                                                                                                        gATGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATgcg    >  1:207982/1‑139 (MQ=255)
                                                                                                          tGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACg                                                                        <  2:71435/69‑1 (MQ=255)
                                                                                                          tGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACg                                                                        >  1:71435/1‑69 (MQ=255)
                                                                                                          tGCGGGAGGATGCCACTTTAATTAATTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGcc  >  1:331235/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                 
GGCGTACAGCACGCCGCTAATCATTACTGGCATCCATTGATCTGTCAGAAGAAAGTTCAGGAAGCAGACATGATTGCTCCCGGCATCCCGTCATTATCAGAGAGGATGCGGGAGGATGCCACTTTAATTAACTAGTTATTAACATTAAAAAGATGGCTTACGCTTTCGGCTAACGTTGTCGTTGGGTGGCCAATCAATTTGCTAAGTGTTTTGCTGTCATCAAACAGACCGCCTTTCGATGCGCC  >  NC_000913/4433033‑4433277

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