Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,064,803	T→C	D105G (GAT→GGT)	yihS ←	sulfoquinovose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,064,803	0	T	C	93.6%	144.5 / ‑0.1	47	D105G (GAT→GGT)	yihS	sulfoquinovose isomerase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (22/22);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

CTTCGCTCCAGAAATATTTCTCGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGG  >  NZ_CP009273/4064669‑4064935                                                                                                                                       |                                                                                                                                     cttcGCTCCAGAAATATTTCTCGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACccac                                                                                                                                  >  2:752880/1‑139 (MQ=255)     gCTCCAGAAATATTTCTCGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACg                                                                                                                              <  1:752880/139‑1 (MQ=255)              atatTTCTCGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCa                                                                                                                     >  1:855063/1‑139 (MQ=255)                    ctcGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCAc                                                                                                                <  2:1088697/138‑1 (MQ=255)                    ctcGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCAc                                                                                                                <  1:599177/138‑1 (MQ=255)                       gATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACg                                                                                                               <  1:912716/136‑1 (MQ=255)                       gATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACg                                                                                                               <  2:957435/136‑1 (MQ=255)                                     tgtAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCCCCACGCCCTCGTCATTCACGCAGGCATCCCAgccgcc                                                                                              <  2:22846/139‑1 (MQ=255)                                            gaccagcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGCCGGCCCCCAGCAGAGCAAAGAAGTGTTGATAGCCTTGTTTGGAGCCCCCCACCACGCCCTCGCAATTAACGAAGGCATACCAGCCGCCATAtttt                                                                                       <  1:1237/136‑1 (MQ=255)                                             agcagcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGACCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATAttttt                                                                                      >  2:1237/1‑139 (MQ=255)                                               cagcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTa                                                                                    >  2:604054/1‑139 (MQ=255)                                                 gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACg                                                                                                               >  1:699526/1‑110 (MQ=255)                                                 gcTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACg                                                                                                               <  2:699526/110‑1 (MQ=255)                                                     gcgcgcTTCCGGGTTACCCGTTGTGACGGCGCTGGCGGCCCCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATCCCAGCCGCCATATTTTTTATCACg                                                                               <  1:27749/138‑1 (MQ=255)                                                      cgcgcTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCa                                                                             <  1:717275/139‑1 (MQ=255)                                                             ccGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCCCCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTgcg                                                                        <  1:477255/137‑1 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCCCCACGCCCTCGTCATTCACGCAGGCATACCAgccgcc                                                                                              <  2:228305/107‑1 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCAc                                                                                                                >  1:999846/1‑89 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCAc                                                                                                                <  2:999846/89‑1 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAgccgcc                                                                                              >  1:228305/1‑107 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATgg                                                              >  2:655290/1‑139 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATgg                                                              >  1:1250259/1‑139 (MQ=255)                                                                     cccGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATgg                                                              >  1:374345/1‑139 (MQ=255)                                                                         ttGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCAtt                                                                                                                   >  1:300683/1‑82 (MQ=255)                                                                         ttGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCAtt                                                                                                                   <  2:300683/82‑1 (MQ=255)                                                                              aCGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATGGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATAc                                                     <  2:569788/139‑1 (MQ=255)                                                                                ggCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCg                                                   >  1:953224/1‑139 (MQ=255)                                                                                ggCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCg                                                   >  2:292426/1‑139 (MQ=255)                                                                                ggCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCg                                                   >  1:27250/1‑139 (MQ=255)                                                                                             cccAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATcaaccaac                                       >  2:983991/1‑138 (MQ=255)                                                                                                accggagCAAAGAGGTGTTGATAGCCCTGTTTGGAGGCACCCACCCCGCCCTCGTCATTCCCGCAGGCTTACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGt                                   <  2:410440/135‑1 (MQ=255)                                                                                                 gcagAGCAAAGAAGTGTTGATAGCCCTGTTTGGCGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGt                                   >  2:578134/1‑138 (MQ=255)                                                                                                     agCAAAGAAGTTTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAgcg                              <  2:374345/139‑1 (MQ=255)                                                                                                     agCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGCTACCGTGATCAACCAACGAGTAAgcg                              >  2:174588/1‑139 (MQ=255)                                                                                                       cAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATa                                                      <  1:670731/113‑1 (MQ=255)                                                                                                       cAAAGAAGTGTTGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATa                                                      >  2:670731/1‑113 (MQ=255)                                                                                                         aagaagTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAgcgc                             >  1:897335/1‑136 (MQ=255)                                                                                                         aagaagTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAgcgc                             <  2:897335/136‑1 (MQ=255)                                                                                                           gatgtgtTTATAGCCCTGTTTGGAGGCTCCCCCCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGt                        <  2:957325/136‑1 (MQ=255)                                                                                                                 ttGATAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGAcc                   >  2:1269206/1‑138 (MQ=255)                                                                                                                    aTAGCCCTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGTCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATc               >  2:1075777/1‑139 (MQ=255)                                                                                                                         ccTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATcaaccaac                                       >  2:964423/1‑110 (MQ=255)                                                                                                                         ccTGTTTGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATcaaccaac                                       <  1:964423/110‑1 (MQ=255)                                                                                                                             tttGGAGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCgcagc       <  1:1125239/138‑1 (MQ=255)                                                                                                                               tGGAGGCCCCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCaa     <  1:874733/138‑1 (MQ=255)                                                                                                                                  aGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACgg  >  2:1006999/1‑138 (MQ=255)                                                                                                                                  aGGCACCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACgg  >  1:909642/1‑138 (MQ=255)                                                                                                                                       |                                                                                                                                     CTTCGCTCCAGAAATATTTCTCGATAATTTCAATGGTGTAATCGAGCAGCTTGCGCGCTTCCGGGTGACCCGTTGTGACGGCGCTGGCGGCACCCAGCAGAGCAAAGAAGTGTTGATAGCCCTGTTTGGAGGCATCCACCACGCCCTCGTCATTCACGCAGGCATACCAGCCGCCATATTTTTTATCACGCAGTGCGCCGTTCATGGCTTTGATACCGTGATCAACCAACGAGTAAGCGCCAGGTCGACCCATCGCCGCAGCAACGG  >  NZ_CP009273/4064669‑4064935

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑