Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,750,169	C→T	12.0%	T228M (ACG→ATG)	yfjB →	NAD kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,750,169	0	C	T	12.0%	216.6 / 16.0	75	T228M (ACG→ATG)	yfjB	NAD kinase
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (6/3);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCACTGCC  >  W3110S.gb/2750104‑2750238                                                                  |                                                                      tACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATc                                                                   >  1:1937118/1‑70 (MQ=255)      tGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCt                                                              >  1:662686/1‑70 (MQ=255)             cATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGttttt                                                      >  1:1524916/1‑71 (MQ=255)             cATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGATCCGTCTGCGtttt                                                       >  1:86778/1‑70 (MQ=255)                 ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGtttt                                                       >  1:696965/1‑66 (MQ=255)                 ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCa                                                  >  1:694885/1‑71 (MQ=255)                 ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCa                                                  >  1:582312/1‑71 (MQ=255)                 ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCa                                                  <  1:340570/71‑1 (MQ=255)                 ttCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCa                                                  >  1:525671/1‑71 (MQ=255)                 ttCACGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCa                                                  >  1:1223029/1‑71 (MQ=255)                   cccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATc                                                >  1:2407305/1‑71 (MQ=255)                   cccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATc                                                >  1:2690074/1‑71 (MQ=255)                    ccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATcc                                                                  <  1:1863691/52‑1 (MQ=255)                    ccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATcc                                                                  <  1:2699566/52‑1 (MQ=255)                    ccGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATcc                                                                  <  1:2775465/52‑1 (MQ=255)                     cGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGc                                              <  1:713330/71‑1 (MQ=255)                      gCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGcc                                             <  1:1698359/71‑1 (MQ=255)                        aTACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                                           <  1:1560358/71‑1 (MQ=255)                        aTACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                                           <  1:704736/71‑1 (MQ=255)                        aTACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                                           <  1:193846/71‑1 (MQ=255)                        aTACGTTGT‑AGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                                           <  1:2638225/70‑1 (MQ=255)                         tACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATcc                                                                  <  1:1599222/47‑1 (MQ=255)                           cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCATGAt                                                                    >  1:238343/1‑43 (MQ=255)                           cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGc                                                            >  1:1358218/1‑51 (MQ=255)                           cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc                                        >  1:2430375/1‑71 (MQ=255)                           cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc                                        >  1:2547942/1‑71 (MQ=255)                           cGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAAc                                        >  1:2558702/1‑71 (MQ=255)                               gTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACg                                       <  1:775696/68‑1 (MQ=255)                                 cAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGt                                           <  1:2679687/62‑1 (MQ=255)                                  aGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTg                                  >  1:2190214/1‑70 (MQ=255)                                  aGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTg                                  >  1:1912529/1‑70 (MQ=255)                                    cACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGaa                               <  1:2235610/71‑1 (MQ=255)                                      cGACCACTGGTCATAAACAGCAGCAGCATTATCCGTCTGCGTTTTTCGCATCGCCGTAACGAc                                     <  1:465833/63‑1 (MQ=255)                                         ccACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGttttt                                                      >  1:2428207/1‑43 (MQ=255)                                          cACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                         >  1:2570558/1‑71 (MQ=255)                                          cACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                         <  1:828562/71‑1 (MQ=255)                                          cACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGt                         <  1:2622022/71‑1 (MQ=255)                                           aCTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGtt                        <  1:2367800/71‑1 (MQ=255)                                               gTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa                    >  1:2658512/1‑71 (MQ=255)                                                tCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAAt                             <  1:1044962/61‑1 (MQ=255)                                                tCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGc                      <  1:580020/68‑1 (MQ=255)                                                 cATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa                    >  1:2109776/1‑69 (MQ=255)                                                 cATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGa                    >  1:271004/1‑69 (MQ=255)                                                 cATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTAGCATCGCCGTAACGACCTGGAAATCAGTTGCGa                    >  1:1551506/1‑69 (MQ=255)                                                   tAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGa                                <  1:1490988/55‑1 (MQ=255)                                                   tAAACAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGa                                <  1:2062752/55‑1 (MQ=255)                                                   tAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATc                            <  1:2644903/59‑1 (MQ=255)                                                   tAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATc                            <  1:537529/59‑1 (MQ=255)                                                     aaCAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:1010845/1‑69 (MQ=255)                                                     aaCAGCAGCAGCATGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:332301/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGAc                   >  1:483977/1‑66 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:428890/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:874582/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:2600873/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:1708228/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:2692757/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGc                >  1:2077020/1‑69 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa              <  1:1147556/71‑1 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa              <  1:17990/71‑1 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa              <  1:1442194/71‑1 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa              <  1:693277/71‑1 (MQ=255)                                                     aaCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGAGACAGc                >  1:1352260/1‑69 (MQ=255)                                                      aCAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGAc                                     <  1:1352166/47‑1 (MQ=255)                                                       cagcagcagCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCTACAGCCAGa            <  1:1223776/71‑1 (MQ=255)                                                         gcagcagcaTGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAg                          >  1:336616/1‑55 (MQ=255)                                                           agcagcaCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGc        <  1:1582093/71‑1 (MQ=255)                                                           agcagcaCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGc        <  1:1571723/71‑1 (MQ=255)                                                            gcagcaCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa              <  1:2332781/64‑1 (MQ=255)                                                            gcagcaCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCa              <  1:582219/64‑1 (MQ=255)                                                              agcaCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCACt     <  1:217171/71‑1 (MQ=255)                                                              agcaCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCACt     >  1:415031/1‑71 (MQ=255)                                                               gcaCGATCCGTCTGCGTTTTTAGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCa       >  1:1380331/1‑68 (MQ=255)                                                                caCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGaaa                              >  1:375717/1‑44 (MQ=255)                                                                 aCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCACTGcc  <  1:2773405/71‑1 (MQ=255)                                                                 aCGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCACTGcc  <  1:262556/71‑1 (MQ=255)                                                                  cGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAg         <  1:1896603/63‑1 (MQ=255)                                                                  |                                                                      TACCCTGGTGCCCATGTTCCCGCATACGTTGTCAGCACGACCACTGGTCATAAACAGCAGCAGCACGATCCGTCTGCGTTTTTCGCATCGCCGTAACGACCTGGAAATCAGTTGCGACAGCCAGATAGCACTGCC  >  W3110S.gb/2750104‑2750238

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