Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,141,409 C→T 28.6% A79A (GCG→GCA hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,141,4090CT28.6% 87.9 / 40.7 56A79A (GCG→GCAhybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base C (23/17);  new base T (3/13);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

GCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCAC  >  W3110S.gb/3141342‑3141477
                                                                   |                                                                    
gCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCTCGTTAACCGGAACGTCGATATTCAGCGCa                                                                   >  1:3066549/1‑71 (MQ=255)
gCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCa                                                                   >  1:38989/1‑71 (MQ=255)
 cGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCa                                                                   <  1:3318166/70‑1 (MQ=255)
  gtggtgTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGctct                                                                 <  1:3944111/71‑1 (MQ=255)
  gtggtgTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACt                                                                 <  1:206333/71‑1 (MQ=255)
  gtggtgTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACt                                                                 <  1:1826507/71‑1 (MQ=255)
      tgtgCGCAGCCAGAATGATGTTACGGATTTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCt                                                             <  1:1564323/71‑1 (MQ=255)
      tgtgCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCt                                                             <  1:715083/71‑1 (MQ=255)
          cgcAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCg                                                         >  1:3868178/1‑71 (MQ=255)
              gCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACg                                                     <  1:965545/71‑1 (MQ=255)
              gCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACg                                                     <  1:1916544/71‑1 (MQ=255)
                   aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACgg                                                <  1:183271/71‑1 (MQ=255)
                   aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                <  1:1913173/71‑1 (MQ=255)
                   aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACg                                                 >  1:2526207/1‑70 (MQ=255)
                   aaTGATGTTACGGATGTATTGCGCGTTAACCGGAACGTAGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                <  1:131425/71‑1 (MQ=255)
                     tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:2656165/1‑69 (MQ=255)
                     tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:2598945/1‑69 (MQ=255)
                     tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:2842338/71‑1 (MQ=255)
                     tgatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:2659611/71‑1 (MQ=255)
                      gatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:910236/1‑68 (MQ=255)
                      gatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACgg                                                >  1:1684776/1‑68 (MQ=255)
                      gatgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              <  1:3649702/70‑1 (MQ=255)
                        tgTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCt                                                             <  1:1108980/53‑1 (MQ=255)
                            aCGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAgcg                                       >  1:2741870/1‑71 (MQ=255)
                            aCGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAgc                                        >  1:721700/1‑70 (MQ=255)
                                 tGTATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACt                                                                 <  1:665202/40‑1 (MQ=255)
                                   tATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGt                                <  1:1568077/71‑1 (MQ=255)
                                   tATTGCGCGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGt                                <  1:3381152/71‑1 (MQ=255)
                                     ttGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGaa                                              >  1:2710161/1‑55 (MQ=255)
                                       gcgcGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:1151767/1‑69 (MQ=255)
                                       gcgcGTTAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:3209269/1‑69 (MQ=255)
                                       gcgcGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTgg                              >  1:159574/1‑69 (MQ=255)
                                        cgcgTTAACCGGAACGTCGATATTCAGCGCACTTTCTGcag                                                         >  1:3105811/1‑39 (MQ=255)
                                        cgcgTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAg                                 >  1:3857124/1‑65 (MQ=255)
                                            ttAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACgg                                                <  1:2487899/46‑1 (MQ=255)
                                            ttAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACgg                                                <  1:1757990/46‑1 (MQ=255)
                                            ttAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACa                       <  1:2316352/71‑1 (MQ=255)
                                            ttAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACa                       <  1:1245332/71‑1 (MQ=255)
                                            ttAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACa                       <  1:544668/71‑1 (MQ=255)
                                             tAACCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACgg                                                <  1:505628/45‑1 (MQ=255)
                                              aaCCGGAACGTCGATATTCAGTGCACTTTCTGCCGCACGAACGGaa                                              >  1:2189148/1‑46 (MQ=255)
                                              aaCCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgta                            >  1:690787/1‑64 (MQ=255)
                                              aaCCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGgta                            >  1:70025/1‑64 (MQ=255)
                                                       gTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGa            <  1:240579/71‑1 (MQ=255)
                                                       ctCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGa            <  1:738940/70‑1 (MQ=255)
                                                          gATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATAc         >  1:494058/1‑71 (MQ=255)
                                                          gATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATAc         <  1:3840185/71‑1 (MQ=255)
                                                          gATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATAc         >  1:817801/1‑71 (MQ=255)
                                                           atatTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATAc         <  1:3419331/70‑1 (MQ=255)
                                                            tatTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGa            <  1:756591/66‑1 (MQ=255)
                                                             atTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGtt      >  1:4966/1‑71 (MQ=255)
                                                             atTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGtt      >  1:1979164/1‑71 (MQ=255)
                                                             atTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGtt      >  1:2565729/1‑71 (MQ=255)
                                                                cAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGt       >  1:3199180/1‑67 (MQ=255)
                                                                cAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGca   >  1:815239/1‑71 (MQ=255)
                                                                 aGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGcac  <  1:1476752/71‑1 (MQ=255)
                                                                   |                                                                    
GCGTGGTGTGCGCAGCCAGAATGATGTTACGGATGTATTGCGCGTTAACCGGAACGTCGATATTCAGCGCACTTTCTGCCGCACGAACGGAAGACAGCGCGTGAGTGGTAGTACATACGCCACAGATACGTTGCAC  >  W3110S.gb/3141342‑3141477

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: