Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,200,991	A→T	5.9%	intergenic (‑26/‑216)	ygiF ← / → ygiM	inorganic triphosphatase/SH3 domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,200,991	0	A	T	5.9%	243.1 / 2.7	66	intergenic (‑26/‑216)	ygiF/ygiM	inorganic triphosphatase/SH3 domain protein
Reads supporting (aligned to +/- strand):  ref base A (31/31);  new base T (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTCAGCAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCAGTGGCAAGACCTGACACAGAGGATTCCGATTGATGTC  >  NC_000913/3200850‑3201137                                                                                                                                              |                                                                                                                                                   atatTCAGCAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACTTAtgtg                                                                                                                                              <  2:152072/148‑1 (MQ=255)      cAGCAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAg                                                                                                                                                <  3:43289/141‑1 (MQ=255)      cAGCAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAg                                                                                                                                                >  4:43289/1‑141 (MQ=255)         cAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTg                                                                                                                                              >  1:693673/1‑140 (MQ=255)         cAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTg                                                                                                                                              <  2:693673/140‑1 (MQ=255)                 cGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAAc                                                                                                                              >  3:425095/1‑148 (MQ=255)                         tGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACg                                                                                                                                    >  4:169749/1‑134 (MQ=255)                         tGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACg                                                                                                                                    <  3:169749/134‑1 (MQ=255)                                  ccgccCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTg                                                                                                                                              >  2:1122913/1‑115 (MQ=255)                                  ccgccCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTg                                                                                                                                              <  1:1122913/115‑1 (MQ=255)                                        aGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGt                                                                                                       <  4:425095/148‑1 (MQ=255)                                               tGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTaa                                                                                                                          <  4:90927/122‑1 (MQ=255)                                               tGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTaa                                                                                                                          >  3:90927/1‑122 (MQ=255)                                                   aTGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATa                                                                                                                                                  <  3:585518/94‑1 (MQ=255)                                                   aTGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATa                                                                                                                                                  >  4:585518/1‑94 (MQ=255)                                                          cGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAAt                                                                                                                         <  3:741026/112‑1 (MQ=255)                                                          cGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAAt                                                                                                                         >  4:741026/1‑112 (MQ=255)                                                             aaCGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATa                                                                                                                                                  <  2:207569/84‑1 (MQ=255)                                                             aaCGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATa                                                                                                                                                  >  1:207569/1‑84 (MQ=255)                                                              aCGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCtttt                                                                                                               >  3:450079/1‑118 (MQ=255)                                                              aCGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCtttt                                                                                                               <  4:450079/118‑1 (MQ=255)                                                                    cAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCttttt                                                                                                              >  2:1032741/1‑113 (MQ=255)                                                                    cAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCttttt                                                                                                              <  1:1032741/113‑1 (MQ=255)                                                                     aaCGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTc                                                                                                          >  2:1115075/1‑116 (MQ=255)                                                                     aaCGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTc                                                                                                          <  1:1115075/116‑1 (MQ=255)                                                                       cGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCg                                                                                                                      <  1:222577/102‑1 (MQ=255)                                                                       cGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCg                                                                                                                      >  2:222577/1‑102 (MQ=255)                                                                          cACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTaa                                                                                                                          >  3:884921/1‑95 (MQ=255)                                                                          cACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTaa                                                                                                                          <  4:884921/95‑1 (MQ=255)                                                                          cACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTa                                                                                                                           >  1:1025588/1‑94 (MQ=255)                                                                          cACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTa                                                                                                                           <  2:1025588/94‑1 (MQ=255)                                                                            cTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACg                                                                                                                                           >  2:565854/1‑76 (MQ=255)                                                                            cTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACg                                                                                                                                           <  1:565854/76‑1 (MQ=255)                                                                            cTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGa                                                                                                                     <  4:130698/98‑1 (MQ=255)                                                                            cTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGa                                                                                                                     >  3:130698/1‑98 (MQ=255)                                                                                  ttAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCaa                                                                        <  3:534990/137‑1 (MQ=255)                                                                                  ttAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCaa                                                                        >  4:534990/1‑137 (MQ=255)                                                                                        aTAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCaa                                                                                    >  2:430230/1‑119 (MQ=255)                                                                                        aTAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCaa                                                                                    <  1:430230/119‑1 (MQ=255)                                                                                        aTAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTg                                                                                >  4:213160/1‑123 (MQ=255)                                                                                        aTAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTg                                                                                <  3:213160/123‑1 (MQ=255)                                                                                            aCTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTAttt                                                                                            <  2:442649/107‑1 (MQ=255)                                                                                            aCTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTAttt                                                                                            >  1:442649/1‑107 (MQ=255)                                                                                               ttAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTg                                                                                                            >  2:870615/1‑88 (MQ=255)                                                                                               ttAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTg                                                                                                            <  1:870615/88‑1 (MQ=255)                                                                                               ttAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGt                                                                                                       <  4:746381/93‑1 (MQ=255)                                                                                               ttAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGt                                                                                                       >  3:746381/1‑93 (MQ=255)                                                                                               ttAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTAttt                                                                                            <  3:1106843/104‑1 (MQ=255)                                                                                               ttAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTAttt                                                                                            >  4:1106843/1‑104 (MQ=255)                                                                                                   ttCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGcc                                                                              >  2:66904/1‑114 (MQ=255)                                                                                                   ttCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGcc                                                                              <  1:66904/114‑1 (MQ=255)                                                                                                                  cATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGAt                                                                                                   >  1:234811/1‑78 (MQ=255)                                                                                                                  cATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGAt                                                                                                   <  2:234811/78‑1 (MQ=255)                                                                                                                       ccTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAAt                                                                                                                         >  2:261655/1‑51 (MQ=255)                                                                                                                       ccTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAAt                                                                                                                         <  1:261655/51‑1 (MQ=255)                                                                                                                         tttACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGAt                                                                                                   >  3:1084628/1‑71 (MQ=255)                                                                                                                         tttACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGAt                                                                                                   <  4:1084628/71‑1 (MQ=255)                                                                                                                         tttACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCa                                       >  3:975856/1‑131 (MQ=255)                                                                                                                         tttACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCa                                       <  4:975856/131‑1 (MQ=255)                                                                                                                               tATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCAGTGGCaa                                >  2:200005/1‑132 (MQ=255)                                                                                                                               tATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCAGTGGCaa                                <  1:200005/132‑1 (MQ=255)                                                                                                                                       aTGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTaaa                                                  >  2:237276/1‑106 (MQ=255)                                                                                                                                       aTGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTaaa                                                  <  1:237276/106‑1 (MQ=255)                                                                                                                                        cgTCACTTATGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACgg                                                                     >  1:28045/2‑86 (MQ=255)                                                                                                                                        cgTCACTTATGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACgg                                                                     <  2:28045/85‑1 (MQ=255)                                                                                                                                        cgTCACTTATGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCAGTGGCAAGACCTGACACAGAGGATTCCGATTg       >  1:216122/2‑148 (MQ=255)                                                                                                                                             caTAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGAc                                                           <  2:138984/91‑1 (MQ=255)                                                                                                                                             caTAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGAc                                                           >  1:138984/1‑91 (MQ=255)                                                                                                                                             caTAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGt                                               >  1:651168/1‑103 (MQ=255)                                                                                                                                             caTAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGt                                               <  2:651168/103‑1 (MQ=255)                                                                                                                                             caTAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCAGTGGCAAGACCTGACACAGAGGATTCCGATTGATGTc  >  1:93068/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                   ATATTCAGCAACTGCACGGGGTCATGGTGCTCGCCGCCCAGCGTATTGAGATGGTCACGCAACGCCTCAACGGCACTGTGATTAACAATAAACTTTAATTCGATTTCCTGAGCCATGGCCTTTACTTATCGGTTATGTCACATAAGTGACGATGAACGGCGAACTTAATGCGATCTTTTTTGTCAGTAGATAGTATTTTGCGCCAAATTGCCATGCAACGGGCAATTTGACGGGCGTAAAAGTTTGAAGCAGTGGCAAGACCTGACACAGAGGATTCCGATTGATGTC  >  NC_000913/3200850‑3201137

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