Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_KO_pgi	2,513,351	C→A	L104M (CTG→ATG)	nupC →	nucleoside (except guanosine) transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_KO_pgi	2,513,351	0	C	A	94.9%	118.8 / ‑3.8	40	L104M (CTG→ATG)	nupC	nucleoside (except guanosine) transporter
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (19/18);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCTCCGAAATGTTCGAAAAACTGCTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCACTGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAAAACTTTA  >  NC_000913_3_KO_pgi/2513216‑2513495                                                                                                                                        |                                                                                                                                                 ggCTTCTCCGAAATGTTCGAAAAACTGCTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATtct                                                                                                                                     >  1:489389/1‑149 (MQ=255)              tGTTCGAAAAACTGCTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATAGGTATTCTCCAGCACATTcgc                                                                                                                        >  2:95469/1‑149 (MQ=255)                ttCGAAAAACTGCTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCgt                                                                                                                      <  2:489389/149‑1 (MQ=255)                           gCTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGa                                                                                                           <  1:536181/149‑1 (MQ=255)                            cTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGat                                                                                                          >  2:99620/1‑149 (MQ=255)                            cTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGat                                                                                                          >  1:510287/1‑149 (MQ=255)                                       cAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCTTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCaa                                                                                               <  2:672041/149‑1 (MQ=255)                                        aaCGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGAACAAGGCCTAGCCTTATCCTCCCTGACGATAATCTGCGCAATCACCCGAATCGCTGCATTCATCCGTTTTTTCCAGCACATTCGCTTGTTCACGTTGATCATCCGCGCAAt                                                                                              >  2:355521/1‑149 (MQ=255)                                              ggaaCTAATTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGttt                                                                                        <  2:214802/146‑1 (MQ=255)                                               ggACTTACTTCGGCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTc                                                                                       <  1:403011/149‑1 (MQ=255)                                                           tCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCaaa                                                                           >  2:137093/1‑149 (MQ=255)                                                           tCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCaaa                                                                           <  2:510287/149‑1 (MQ=255)                                                           tCTTTGCGAGCATGAAGGAGCAAGGGCTGGCGTTCTTCTTCCTGAAAGGGCTGTGCCCAATCGTCTTTATCTCTGCAAGGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATGCGCGCAATTGGTTTCCTGCTCTCCaaa                                                                           <  2:84431/149‑1 (MQ=255)                                                                 gTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTAGCGAGTTACCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAAc                                                                     >  1:20572/1‑149 (MQ=255)                                                                        gaatgaTCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTCTCCTGCTCTCCAAAGTCAACGGCATgg                                                              >  1:638313/1‑149 (MQ=255)                                                                          atgaTCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGc                                                            >  1:12852/1‑149 (MQ=255)                                                                              tCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAAc                                                        <  1:99620/149‑1 (MQ=255)                                                                                    ccTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCACTGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAAt                                                  <  2:617073/149‑1 (MQ=255)                                                                                     cTGTCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATc                                                 <  2:332232/149‑1 (MQ=255)                                                                                                 ttcCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGt                                     <  1:22149/149‑1 (MQ=255)                                                                                                   ccTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCa                                   >  1:116972/1‑149 (MQ=255)                                                                                                             cTGTGCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGAt                         >  2:152720/1‑149 (MQ=255)                                                                                                              tgtgCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATc                                                 >  2:329619/1‑124 (MQ=255)                                                                                                              tgtgCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATc                                                 <  1:329619/124‑1 (MQ=255)                                                                                                               gtgCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTc                       >  2:218457/1‑149 (MQ=255)                                                                                                                 gCCCAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTg                     <  1:370838/149‑1 (MQ=255)                                                                                                                  cccAATCGTCTTTATCTCTGCGATGATCTGTATTCTCCAGCACAATTGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTgg                    <  2:117036/149‑1 (MQ=255)                                                                                                                  cccAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGCTTCTgg                    >  2:216692/1‑149 (MQ=255)                                                                                                                  cccAAGCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGGTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTgg                    <  2:46945/149‑1 (MQ=255)                                                                                                                    cAATCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGt                  >  1:567806/1‑149 (MQ=255)                                                                                                                     aaTCGTCTTTATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGCTTCTGGGTc                 >  1:43400/1‑149 (MQ=255)                                                                                                                            tttATCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCAGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAGCCTGGAATCCTTTAACGCCGCAACCCCCCTCATTCTGGGTCAGTCTGa          >  2:630503/1‑149 (MQ=255)                                                                                                                               aTCTCTGCACTGATCGCTATTCTCCAGCACATTCGTGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCATGATTCTGGGTCAGTCTGaaaa       <  2:352955/149‑1 (MQ=255)                                                                                                                               aTCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTg                     >  1:421062/1‑135 (MQ=255)                                                                                                                               aTCTCTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTg                     <  2:421062/135‑1 (MQ=255)                                                                                                                                tctcTGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAAAAc      <  1:95469/149‑1 (MQ=255)                                                                                                                                 cttcgCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAAAACt     <  1:47682/145‑1 (MQ=255)                                                                                                                                 ctctGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAATACt     >  2:201321/1‑149 (MQ=255)                                                                                                                                  tctGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAAAACtt    <  1:417684/149‑1 (MQ=255)                                                                                                                                    tGCAATGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAAAACTTTa  >  1:591494/1‑149 (MQ=255)                                                                                                                                        |                                                                                                                                                 GGCTTCTCCGAAATGTTCGAAAAACTGCTCGGATTTGCCAACGAAGGGACTAACTTCGTCTTTGGTAGCATGAATGATCAAGGCCTGGCATTCTTCTTCCTGAAAGTGCTGTGCCCAATCGTCTTTATCTCTGCACTGATCGGTATTCTCCAGCACATTCGCGTGTTGCCGGTGATCATCCGCGCAATTGGTTTCCTGCTCTCCAAAGTCAACGGCATGGGCAAACTGGAATCCTTTAACGCCGTCAGCTCCCTGATTCTGGGTCAGTCTGAAAACTTTA  >  NC_000913_3_KO_pgi/2513216‑2513495

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