Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,295,693	C→T	14.3%	A18T (GCG→ACG)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,295,693	0	C	T	14.3%	168.5 / 16.3	63	A18T (GCG→ACG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (7/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTT  >  W3110S.gb/2295626‑2295758                                                                    |                                                                  tcttcAGTGAGCTGACGGAACTGTTGTTCCTGTTCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAg                                                                >  1:1575275/1‑71 (MQ=255) tcttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAg                                                                >  1:543013/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                               >  1:2655324/1‑71 (MQ=255)  cttcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                               >  1:3291662/1‑71 (MQ=255)    tcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcgc                                                             >  1:391990/1‑71 (MQ=255)    tcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcgc                                                             >  1:2300879/1‑71 (MQ=255)    tcAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgcgc                                                             >  1:1304119/1‑71 (MQ=255)     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                               <  1:3220236/68‑1 (MQ=255)     cAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAgc                                                               <  1:2065496/68‑1 (MQ=255)           gCTGACGGAACTTTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCg                                                      <  1:4033691/71‑1 (MQ=255)              gACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCg                                                      >  1:3745855/1‑68 (MQ=255)                 ggAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatc                                                <  1:2370778/71‑1 (MQ=255)                  gAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatca                                               <  1:2667077/71‑1 (MQ=255)                  gAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatca                                               <  1:3703362/71‑1 (MQ=255)                     cTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatc                                                <  1:4016821/67‑1 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCACCAGCGCTGAGCCGGAGatcat                                              >  1:4004625/1‑67 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCACCAGCGCTGAGCCGGAGatcat                                              >  1:3824595/1‑67 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCACCAGCGCTGAGCCGGAGatcat                                              >  1:3412515/1‑67 (MQ=255)                       gttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCACCAGCGCTGAGCCGGAGatcat                                              >  1:72531/1‑67 (MQ=255)                        ttgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATcatcag                                           <  1:2431121/69‑1 (MQ=255)                         tgttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagca                                         >  1:3800909/1‑70 (MQ=255)                           ttCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatc                                       <  1:828886/70‑1 (MQ=255)                             ccTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATcatcagcatcag                                     <  1:1821207/70‑1 (MQ=255)                               tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGat                                                 <  1:1129602/56‑1 (MQ=255)                               tgtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGat                                                 <  1:2041538/56‑1 (MQ=255)                                gtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCg                                                      <  1:3304176/50‑1 (MQ=255)                                gtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACg                                 >  1:2871285/1‑71 (MQ=255)                                gtgCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACg                                 >  1:1310105/1‑71 (MQ=255)                                  gCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGcc                               >  1:3102885/1‑71 (MQ=255)                                   cTTCATCTTTAAACTGCAACACATCGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGcc                               >  1:4028785/1‑70 (MQ=255)                                       aTCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATa                          <  1:1958256/71‑1 (MQ=255)                                        tCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgag                                                   <  1:2597939/45‑1 (MQ=255)                                        tCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgag                                                   <  1:2815239/45‑1 (MQ=255)                                          tttAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGga                                                    <  1:776481/42‑1 (MQ=255)                                          tttAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGatcatc                                             >  1:3081073/1‑49 (MQ=255)                                          tttAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGc                                >  1:2101407/1‑62 (MQ=255)                                          tttAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        >  1:1128320/1‑70 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGATCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:3727769/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                                  >  1:1481962/1‑42 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                                  >  1:1760887/1‑42 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGgaga                                                  >  1:1795751/1‑42 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:2408599/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:3047181/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:543731/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:3921222/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:3226327/68‑1 (MQ=255)                                            tAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:3241998/68‑1 (MQ=255)                                            tAAACTGCAACACATCAATGGTCGCCAGCGCTGAGCCGgag                                                   >  1:816238/1‑41 (MQ=255)                                              aaCTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaa                        <  1:425957/66‑1 (MQ=255)                                                          tCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCAt                >  1:2498573/1‑62 (MQ=255)                                                          tCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCAt                >  1:3746588/1‑62 (MQ=255)                                                          tCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcc       >  1:1918485/1‑71 (MQ=255)                                                           cGATGGTCGTCAGCGCTGAGCCGGAGATcatcagcat                                        >  1:481831/1‑37 (MQ=255)                                                           cGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                      >  1:804402/1‑55 (MQ=255)                                                           cGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                      >  1:2845925/1‑55 (MQ=255)                                                           cGATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATaaaaa                      >  1:2146508/1‑55 (MQ=255)                                                            gATGGTCGTCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctc     >  1:2041147/1‑71 (MQ=255)                                                            gATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAAc                     >  1:2056048/1‑55 (MQ=255)                                                              tGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctcct   <  1:3732700/71‑1 (MQ=255)                                                              tGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctcct   <  1:3976474/71‑1 (MQ=255)                                                               ggTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGcctcct   >  1:3829634/1‑70 (MQ=255)                                                                gTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCtt  <  1:738404/70‑1 (MQ=255)                                                                gTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCtt  <  1:2743070/70‑1 (MQ=255)                                                                    |                                                                  TCTTCAGTGAGCTGACGGAACTGTTGTTCCTGTGCTTCATCTTTAAACTGCAACACATCGATGGTCGCCAGCGCTGAGCCGGAGATCATCAGCATCAGCACGCCCAATAAAAACCTCATTGTGCGGCCTCCTT  >  W3110S.gb/2295626‑2295758

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