Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,206,807	A→G	Q371Q (CAA→CAG)	yehM →	putative YehN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,206,807	0	A	G	96.1%	171.3 / NA	51	Q371Q (CAA→CAG)	yehM	putative YehN
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (21/28);  minor base C (1/1);  total (22/29)
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

CCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAG  >  NC_000913/2206680‑2206942                                                                                                                                |                                                                                                                                        cccTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCa                                                                                                                             >  2:797122/1‑140 (MQ=255) cccTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCa                                                                                                                             >  2:746760/1‑140 (MQ=255) cccTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCa                                                                                                                             >  1:443888/1‑140 (MQ=255)     gAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAg                                                                                                                         >  2:680082/1‑140 (MQ=255)      aaCGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGa                                                                                                                        >  2:817979/1‑140 (MQ=255)       aCGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGAt                                                                                                                       <  2:4356/140‑1 (MQ=255)           gcCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCa                                                                                                                   >  2:624220/1‑140 (MQ=255)           gcCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTAAGCCACACCCGCCTGTGCCAAAAGATGTCa                                                                                                                   >  1:270983/1‑140 (MQ=255)              gTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCt                                                                                                                <  1:797122/140‑1 (MQ=255)              gTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCt                                                                                                                <  1:746760/140‑1 (MQ=255)                       gAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTgg                                                                                                       >  2:158050/1‑140 (MQ=255)                         gCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACCTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGaa                                                                                                     >  2:502497/1‑140 (MQ=255)                             tCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTa                                                                                                 <  1:158050/140‑1 (MQ=255)                                               cGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCcgg                                                                               <  1:472172/140‑1 (MQ=255)                                                     cgatccgatTCTGCTTACGTTGCTACACAAATTAGCGTGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGc                                                                         <  1:247829/140‑1 (MQ=255)                                                     cgatccgatTCTGCTAACGTTGATAGACACATTAGCGTGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGc                                                                         <  1:273182/140‑1 (MQ=255)                                                     cgatccgatTCTGCAAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGc                                                                         <  2:146800/140‑1 (MQ=255)                                                     agatccgatTCTGCTAACGTTGGTAGACACATTAGGGGGGGACGGATTGGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGc                                                                         <  2:340970/139‑1 (MQ=255)                                                       atccgatTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCtt                                                                       >  1:696435/1‑140 (MQ=255)                                                       atccgatTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCtt                                                                       >  1:586045/1‑140 (MQ=255)                                                                 gCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATc                                                             <  1:768180/140‑1 (MQ=255)                                                                  cTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTCCCGCAATCTCCATTCCGGCGGCGCTTACCTTAAATct                                                            >  1:340970/1‑139 (MQ=255)                                                                    aaCGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCt                                                          <  1:680082/140‑1 (MQ=255)                                                                    aaCGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCt                                                          <  1:624220/140‑1 (MQ=255)                                                                      cGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTTGAACGTACAGCACTCTTACTTCCGGAGGAGCCTACCTTAAATCGCttt                                                        >  2:394421/1‑140 (MQ=255)                                                                       gTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTa                                                       <  1:86462/140‑1 (MQ=255)                                                                       gTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTa                                                       >  2:642155/1‑140 (MQ=255)                                                                       gTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTa                                                       <  1:179270/140‑1 (MQ=255)                                                                            tAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGCACCGCAATCTCCCTTCTGGCGGAGCATACCTTAAATCGCTTTAccccc                                                  >  2:289244/1‑140 (MQ=255)                                                                                   aTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGc                                           <  1:565155/140‑1 (MQ=255)                                                                                     tAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTg                                         >  2:120604/1‑140 (MQ=255)                                                                                       gCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGgc                                       <  2:586045/140‑1 (MQ=255)                                                                                          ggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCa                                    >  1:163185/1‑140 (MQ=255)                                                                                          ggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCa                                    >  1:410183/1‑140 (MQ=255)                                                                                          ggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCa                                    >  1:745986/1‑140 (MQ=255)                                                                                          ggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCa                                    >  1:727462/1‑140 (MQ=255)                                                                                          ggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACATACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGcaca                                    >  1:547286/1‑140 (MQ=255)                                                                                             gACGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAACCGCTTTACCCCCGATGGGCTGGCGCAAAg                                 >  1:290339/1‑140 (MQ=255)                                                                                             gACGGATTCGGTAAACTTGCCCCTTCTACGCAACACCCGACTCTGCCAGAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAg                                 <  2:270983/140‑1 (MQ=255)                                                                                               cGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTc                               <  1:817979/140‑1 (MQ=255)                                                                                               cGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTc                               <  2:745986/140‑1 (MQ=255)                                                                                               cGGATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTc                               <  1:642155/140‑1 (MQ=255)                                                                                                ggATTCGGTAAACTTGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCa                              <  2:163185/140‑1 (MQ=255)                                                                                                            cTTGCCCCTTCTACGTCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCg                  <  1:289244/140‑1 (MQ=255)                                                                                                             ttGCCCCTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCgg                 <  2:443888/140‑1 (MQ=255)                                                                                                                 cccTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTgg             <  2:460756/140‑1 (MQ=255)                                                                                                                  ccTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGc            <  1:120604/140‑1 (MQ=255)                                                                                                                  ccTTCTACGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGc            <  2:696435/140‑1 (MQ=255)                                                                                                                        aCGCCACAGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACt      <  2:290339/140‑1 (MQ=255)                                                                                                                           ccACAGCCGCCTCTGCCCAAAGCTGTCACCTTCGACCTGTAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGga   <  2:245929/140‑1 (MQ=255)                                                                                                                            cacaGCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGgag  <  1:502497/140‑1 (MQ=255)                                                                                                                                |                                                                                                                                        CCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAG  >  NC_000913/2206680‑2206942

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