Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,140,619 C→T 12.1% A343T (GCA→ACA)  hybC ← hydrogenase 2, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,140,6190CT12.1% 90.9 / 4.3 33A343T (GCA→ACA) hybChydrogenase 2, large subunit
Reads supporting (aligned to +/- strand):  ref base C (11/18);  new base T (2/2);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCCAGGAATATTTCCCGTCGTCAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTA  >  W3110S.gb/3140551‑3140682
                                                                    |                                                               
cACCCAGGAATATTTCCCGTCGTCAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGcc                                                                <  1:1138623/70‑1 (MQ=255)
    cAGGAATATTTCCCGTCGTCAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGTCTgcgg                                                           >  1:2621911/1‑71 (MQ=255)
              tCCCGTCGTCAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTc                                                  <  1:1902047/70‑1 (MQ=255)
                   tcgtcAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCtt                                                <  1:162300/67‑1 (MQ=255)
                     gtcAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGTCTGCGGCGCTTCGTCTTta                                              >  1:2697393/1‑67 (MQ=255)
                     gtcAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATAc                                           >  1:514444/1‑70 (MQ=255)
                        aGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAgg                                       <  1:1837343/71‑1 (MQ=255)
                               ccaTCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCtt                                >  1:2570614/1‑71 (MQ=255)
                                 atcatAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCtt                                                <  1:686295/53‑1 (MQ=255)
                                 atcatAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCtt                                                <  1:2548656/53‑1 (MQ=255)
                                 atcatAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTtata                                            >  1:739229/1‑57 (MQ=255)
                                  tcatAAGCCGGAATGGGGGTGCCTTCCCACGGTGCCTGCGgcgc                                                        <  1:609974/44‑1 (MQ=255)
                                   catAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATAcc                                          >  1:500377/1‑57 (MQ=255)
                                   catAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCt                                 <  1:1305228/66‑1 (MQ=255)
                                   catAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCt                                 <  1:2089805/66‑1 (MQ=255)
                                    atAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCgcgc                           >  1:23479/1‑71 (MQ=255)
                                           ggAATGGTGGTGCCTTCCCACGGTGTCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTg                    <  1:1198096/71‑1 (MQ=255)
                                           ggAATGGTGGTGCCTTCCCACGGTGTCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTg                    <  1:209655/71‑1 (MQ=255)
                                           ggAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTg                    <  1:82984/71‑1 (MQ=255)
                                           ggAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTg                    <  1:2510255/71‑1 (MQ=255)
                                            gAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTg                    <  1:1655380/70‑1 (MQ=255)
                                              aTGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTta                                              <  1:2420427/42‑1 (MQ=255)
                                               tggtggTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTc                               <  1:1449146/56‑1 (MQ=255)
                                               tggtggTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATg                >  1:50158/1‑71 (MQ=255)
                                                 gtggtgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGcc              >  1:142857/1‑71 (MQ=255)
                                                    gtgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTg                    <  1:2204964/62‑1 (MQ=255)
                                                    gtgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGaa                  >  1:1522282/1‑64 (MQ=255)
                                                     tgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTg          <  1:307323/71‑1 (MQ=255)
                                                     tgCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTg          <  1:1087266/71‑1 (MQ=255)
                                                         ttCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCttt                        <  1:1599701/53‑1 (MQ=255)
                                                         ttCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTg          <  1:1214801/67‑1 (MQ=255)
                                                           cccACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGAt        >  1:1655925/1‑67 (MQ=255)
                                                             cACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTa  >  1:1820950/1‑71 (MQ=255)
                                                                ggTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCt                                 >  1:344613/1‑37 (MQ=255)
                                                                    |                                                               
CACCCAGGAATATTTCCCGTCGTCAGACCAACCATCATAAGCCGGAATGGTGGTGCCTTCCCACGGTGCCTGCGGCGCTTCGTCTTTATACCAGGAGTGCTTCGCGCTTTCCTGAATGCCTTTGATCAGGTA  >  W3110S.gb/3140551‑3140682

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: