Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 954,052 C→T 6.2% L279L (CTG→CTA)  focA ← formate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb954,0520CT 6.2% 264.6 / 4.7 80L279L (CTG→CTA) focAformate transporter
Reads supporting (aligned to +/- strand):  ref base C (42/33);  new base T (2/3);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACC  >  W3110S.gb/953983‑954118
                                                                     |                                                                  
cctacctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGc                                                                    <  1:1968087/70‑1 (MQ=255)
 ctacctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCa                                                                   <  1:3637542/70‑1 (MQ=255)
  tacctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCa                                                                   >  1:1086863/1‑69 (MQ=255)
  tacctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAgg                                                                 >  1:3534759/1‑71 (MQ=255)
  tacctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAgg                                                                 >  1:1173708/1‑71 (MQ=255)
   acctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAgg                                                                 >  1:5032734/1‑70 (MQ=255)
    cctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTa                                                               <  1:56712/71‑1 (MQ=255)
    cctTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACCCAGGTa                                                               <  1:8878/71‑1 (MQ=255)
     cttcttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTaa                                                              <  1:3201680/71‑1 (MQ=255)
     cttcttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTaa                                                              <  1:2098606/71‑1 (MQ=255)
     cttcttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTa                                                               <  1:352566/70‑1 (MQ=255)
       tcttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAAt                                                            >  1:314025/1‑71 (MQ=255)
        cttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAAt                                                            >  1:1876698/1‑70 (MQ=255)
        cttAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAAt                                                            >  1:912588/1‑70 (MQ=255)
           aaGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATg                                                           <  1:4825359/68‑1 (MQ=255)
             gTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCa                                                      <  1:2112765/71‑1 (MQ=255)
              tGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAg                                                     <  1:167707/71‑1 (MQ=255)
                  ttttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGTAgg                                                                 <  1:4368815/55‑1 (MQ=255)
                  ttttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGTAgg                                                                 <  1:3783870/55‑1 (MQ=255)
                  ttttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGt                                                    <  1:4185238/68‑1 (MQ=255)
                    ttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGTAGGTa                                                               >  1:4618879/1‑55 (MQ=255)
                    ttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGTAGGTa                                                               >  1:1279234/1‑55 (MQ=255)
                    ttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGAcc                                                        >  1:2510777/1‑62 (MQ=255)
                    ttttATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGt                                                >  1:2145708/1‑70 (MQ=255)
                      ttatttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTca                                              <  1:266905/70‑1 (MQ=255)
                         tttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTcaaccc                                          <  1:2289161/71‑1 (MQ=255)
                         tttaCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTcaacc                                           <  1:203041/70‑1 (MQ=255)
                           taCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTcaaccca                                         <  1:4478532/70‑1 (MQ=255)
                            aCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTcaacccaac                                       >  1:285657/1‑71 (MQ=255)
                            aCTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTcaacccaac                                       >  1:904557/1‑71 (MQ=255)
                                 gTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATg                                                           <  1:636892/46‑1 (MQ=255)
                                  tACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACcaacaa                                 <  1:1217446/71‑1 (MQ=255)
                                  tACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACcaaca                                  >  1:177642/1‑70 (MQ=255)
                                  tACTTCGACAACCATTAATGGTGGTCGCTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACcaaca                                  >  1:1601450/1‑70 (MQ=255)
                                     ttCGACAACCATTAATGGTGGTCGTTTTCACGTAGGTAAATGACCCAg                                                     <  1:3454546/48‑1 (MQ=255)
                                      tCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:378328/1‑71 (MQ=255)
                                       cGACAACCATTAATGGTGGTCGTTTTCACGCAGGTaaa                                                             >  1:4981423/1‑38 (MQ=255)
                                       cGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:4965599/1‑70 (MQ=255)
                                       cGACAACCATTAATGGTGGTAGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccac                            >  1:4696326/1‑71 (MQ=255)
                                        gACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccac                            >  1:5023058/1‑70 (MQ=255)
                                        gACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccac                            >  1:2388846/1‑70 (MQ=255)
                                        gACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:900164/1‑69 (MQ=255)
                                             ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATg                                                 <  1:2865484/44‑1 (MQ=255)
                                             ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATg                                                 <  1:4416035/44‑1 (MQ=255)
                                             ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATg                                                 <  1:2201423/44‑1 (MQ=255)
                                             ccATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATg                                                 <  1:4367342/44‑1 (MQ=255)
                                              cATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATg                                                 <  1:4747701/43‑1 (MQ=255)
                                              cATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGa                      >  1:3122549/1‑70 (MQ=255)
                                              cATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGAt                     <  1:1956035/71‑1 (MQ=255)
                                               aTTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGt                                                >  1:3107091/1‑43 (MQ=255)
                                               aTTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccacc                           <  1:3466430/64‑1 (MQ=255)
                                               aTTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccacc                           <  1:3989902/64‑1 (MQ=255)
                                               aTTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGAt                     >  1:3850494/1‑70 (MQ=255)
                                                 tAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTcaa                                             >  1:2740002/1‑44 (MQ=255)
                                                  aaTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:225140/1‑59 (MQ=255)
                                                  aaTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:1739847/1‑59 (MQ=255)
                                                  aaTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATg                 >  1:1159387/1‑71 (MQ=255)
                                                  aaTGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATg                 >  1:1091888/1‑71 (MQ=255)
                                                         gtCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCg           <  1:3958517/70‑1 (MQ=255)
                                                          tCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGa          <  1:881167/70‑1 (MQ=255)
                                                            gTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAaccaccac                         <  1:564719/53‑1 (MQ=255)
                                                            gTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTg              <  1:5049832/64‑1 (MQ=255)
                                                            gTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTg              <  1:3636040/64‑1 (MQ=255)
                                                             ttttCACGCAGGTAAATGACGCAGTATGTCAACCCAACCAACAAACCACCACCGATAAt                  <  1:3627405/59‑1 (MQ=255)
                                                             ttttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGcc            <  1:2205598/65‑1 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACccaaccaa                                    >  1:4976857/1‑39 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACCcaaccaac                                   >  1:4844142/1‑40 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACCcaaccaac                                   >  1:487713/1‑40 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACaaa                                >  1:2911010/1‑43 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTaa    >  1:4323468/1‑71 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTaa    >  1:3505487/1‑71 (MQ=255)
                                                               ttCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTa     <  1:372548/70‑1 (MQ=255)
                                                                tCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:3972388/1‑45 (MQ=255)
                                                                tCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:3165580/1‑45 (MQ=255)
                                                                tCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:275203/1‑45 (MQ=255)
                                                                tCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAacca                             >  1:964051/1‑45 (MQ=255)
                                                                 cACGCAGGTAAATGACCCAGTATGTCATCCCAACCAACAAACCACCACCGATAATGTAGCCGATCGTAAc   >  1:4155760/1‑70 (MQ=255)
                                                                 cACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGtt               >  1:856403/1‑58 (MQ=255)
                                                                 cACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAAcc  >  1:38659/1‑71 (MQ=255)
                                                                 cACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAAcc  >  1:459082/1‑71 (MQ=255)
                                                                  aCGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCg           <  1:3541263/61‑1 (MQ=255)
                                                                     |                                                                  
CCTACCTTCTTAAGTGGATTTTTTATTTACTGCGTACTTCGACAACCATTAATGGTGGTCGTTTTCACGCAGGTAAATGACCCAGTATGTCAACCCAACCAACAAACCACCACCGATAATGTTGCCGATCGTAACC  >  W3110S.gb/953983‑954118

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: ‑