Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,400,870 T→G 19.0% D84A (GAC→GCC)  nuoI ← NADH:ubiquinone oxidoreductase, chain I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,400,8700TG19.0% 26.8 / 6.0 21D84A (GAC→GCC) nuoINADH:ubiquinone oxidoreductase, chain I
Reads supporting (aligned to +/- strand):  ref base T (13/4);  new base G (2/2);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTACCGCGCAGAGGTTACAGGCTA  >  W3110S.gb/2400803‑2400939
                                                                   |                                                                     
ttcGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGTCtt                                                                    >  1:711213/1‑71 (MQ=255)
ttcGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGGCtt                                                                    >  1:28126/1‑71 (MQ=255)
  cGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGTCttt                                                                   <  1:1947829/70‑1 (MQ=255)
            acaGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCtt                                                        <  1:2185155/71‑1 (MQ=255)
                 aaaTGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTg                                                    >  1:4434908/1‑70 (MQ=255)
                                  ttGATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGAc                                  >  1:977134/1‑71 (MQ=255)
                                    gATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCgaga                                            >  1:1266758/1‑59 (MQ=255)
                                    gATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGAc                                  >  1:3933558/1‑69 (MQ=255)
                                    gATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCg                                >  1:3722818/1‑71 (MQ=255)
                                     aTGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATa                                          >  1:4623492/1‑60 (MQ=255)
                                     aTGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATa                                          >  1:2139692/1‑60 (MQ=255)
                                        cGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCggg                              >  1:2013031/1‑69 (MQ=255)
                                          gAAAAATTCCGGGTACCAGCGACCGGCTTTGGTTTCt                                                            <  1:3941345/37‑1 (MQ=255)
                                               aTTCCGGGTACCAGCGACCGGCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTa                       >  1:3999218/1‑69 (MQ=255)
                                                ttCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATa                                          >  1:3657782/1‑49 (MQ=255)
                                                ttCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATa                                          >  1:1113431/1‑49 (MQ=255)
                                                 tCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGAc                                  >  1:1789673/1‑56 (MQ=255)
                                                   cGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTACcgcg                  <  1:1029490/70‑1 (MQ=255)
                                                    gggTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTACCGCGCa                <  1:3809291/71‑1 (MQ=255)
                                                     ggTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGAc                                  >  1:4134529/1‑52 (MQ=255)
                                                          cAGCGACCGGCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTACCGCGCAGAGGtt          <  1:1055749/71‑1 (MQ=255)
                                                                   tCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTACCGCGCAGAGGTTACAGGCTa  >  1:1972176/1‑70 (MQ=255)
                                                                   |                                                                     
TTCGCACAGACCACAGAAAATGCAGCGTGAGAAGTTGATGCGGAAAAATTCCGGGTACCAGCGACCGTCTTTGGTTTCTGCTTTTTGCAGCGAGATACAGCCGACCGGGCAGGCTACCGCGCAGAGGTTACAGGCTA  >  W3110S.gb/2400803‑2400939

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