Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 639,164 G→A 61.9% P63P (CCG→CCA ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb639,1640GA61.9% 10.2 / 20.1 21P63P (CCG→CCAahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base A (4/9);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACG  >  W3110S.gb/639103‑639231
                                                             |                                                                   
tCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTg                                                            <  1:1715499/71‑1 (MQ=255)
       aCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCa                                                     <  1:826893/71‑1 (MQ=255)
       aCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCa                                                     <  1:690602/71‑1 (MQ=255)
              gACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAgg                                              <  1:1243046/71‑1 (MQ=255)
                cAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCtt                                                                 <  1:905712/50‑1 (MQ=255)
                 aaaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTc                                           >  1:981858/1‑71 (MQ=255)
                  aaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGt                                                                    >  1:1036831/1‑45 (MQ=255)
                     tCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAc                                       <  1:1440046/71‑1 (MQ=255)
                     tCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAc                                       <  1:1748183/71‑1 (MQ=255)
                        ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                                    <  1:841054/71‑1 (MQ=255)
                           ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAgggg                                 >  1:404619/1‑71 (MQ=255)
                             aaaGAAGATAACAGCTTGCCGGTGCGTAAGCCATc                                                                   >  1:827010/1‑35 (MQ=255)
                              aagaagATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCa                              >  1:1234662/1‑71 (MQ=255)
                                    aTAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACg                            <  1:1224213/67‑1 (MQ=255)
                                            ttGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACg                            <  1:633613/59‑1 (MQ=255)
                                            ttGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACg                            <  1:1719067/59‑1 (MQ=255)
                                               ccGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTccc              <  1:704920/70‑1 (MQ=255)
                                                    gCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                       >  1:137578/1‑56 (MQ=255)
                                                    gCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg        >  1:797524/1‑71 (MQ=255)
                                                        aaGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAgg                   >  1:1613413/1‑56 (MQ=255)
                                                          gCCGTCCTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACg  <  1:448320/71‑1 (MQ=255)
                                                             |                                                                   
TCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGGCCACG  >  W3110S.gb/639103‑639231

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