Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,714,423 C→T 14.3% H61Y (CAC→TAC)  atpI → ATP synthase, membrane‑bound accesory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,714,4230CT14.3% 71.7 / 4.6 28H61Y (CAC→TAC) atpIATP synthase, membrane‑bound accesory subunit
Reads supporting (aligned to +/- strand):  ref base C (15/9);  new base T (2/2);  total (17/11)
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.98e-01

TCTCTGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAGGCCGGGTGGCCTGGACATTCGCATTTGGCGAAGCTTTCA  >  W3110S.gb/3714355‑3714486
                                                                    |                                                               
tctcTGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTTAc                                                               <  1:1713845/71‑1 (MQ=255)
tctcTGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTTAc                                                               <  1:1718635/71‑1 (MQ=255)
tctcTGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCAc                                                               <  1:753125/71‑1 (MQ=255)
tctcTGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCAc                                                               <  1:1508541/71‑1 (MQ=255)
  tctGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTTAcg                                                              >  1:2906382/1‑69 (MQ=255)
  tctGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCAc                                                               >  1:1788800/1‑69 (MQ=255)
    tGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAgg                                                           >  1:473462/1‑71 (MQ=255)
       aaTAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCa                                                             >  1:2144168/1‑66 (MQ=255)
        aTAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCa                                                       <  1:3445807/71‑1 (MQ=255)
        aTAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCa                                                       <  1:1884040/71‑1 (MQ=255)
           aGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATa                                                     >  1:1381978/1‑70 (MQ=255)
                   ccTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCa                                                                >  1:1247957/1‑51 (MQ=255)
                   ccTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCa                                                                >  1:2804118/1‑51 (MQ=255)
                      ggCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAg                                         >  1:1896265/1‑71 (MQ=255)
                         aGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAg                                                            >  1:4003515/1‑49 (MQ=255)
                         aGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAg                                                            >  1:2937070/1‑49 (MQ=255)
                          gTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTc                                                                 <  1:1221526/43‑1 (MQ=255)
                          gTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTc                                                                 <  1:1813336/43‑1 (MQ=255)
                           tCTTTCTGCCTAACGTTTTCTTTATGATATTTGCCTGGCGTCAcc                                                              <  1:1704739/45‑1 (MQ=255)
                             tttCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCAc                                                               >  1:2675015/1‑42 (MQ=255)
                             tttCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAGGCCGGGt                                  >  1:2033863/1‑71 (MQ=255)
                             tttCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAGGCCGGGt                                  >  1:277683/1‑71 (MQ=255)
                                     tAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGgcg                                                         <  1:71484/40‑1 (MQ=255)
                                     tAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGa                                            <  1:2060453/53‑1 (MQ=255)
                                       aagTTTTGTTTATGATATTTGCCTGGCGTTACCAGGCGCATACACCAGCGAAAGGCCGGGTGGCCTGGACa                        >  1:2863323/3‑71 (MQ=255)
                                       aCGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGaaa                                          >  1:2156774/1‑53 (MQ=255)
                                       aCGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAGGCC‑‑GTGGCCTGGACAtt                      >  1:2475082/1‑71 (MQ=255)
                                         gTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAGGCCGGGTGGCCTGGACAtt                      <  1:2300216/71‑1 (MQ=255)
                                                        tttGCCTGGCGTCACCAGGCGCATACACCAGCGAAAg                                         >  1:2664927/1‑37 (MQ=255)
                                                             cTGGCGTCACCAGGCGCATACACCAGCGAAAGGCCGGGTGGCCTGGACATTCGCATTTGGCGAAGCTTTCa  <  1:206546/71‑1 (MQ=255)
                                                                    |                                                               
TCTCTGCAATAAGCGGGGGCCTGGCAGTCTTTCTGCCTAACGTTTTGTTTATGATATTTGCCTGGCGTCACCAGGCGCATACACCAGCGAAAGGCCGGGTGGCCTGGACATTCGCATTTGGCGAAGCTTTCA  >  W3110S.gb/3714355‑3714486

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