Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,028,264	T→C	Y36C (TAC→TGC)	yodC ←	protein YodC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,028,264	0	T	C	95.1%	199.7 / ‑0.1	61	Y36C (TAC→TGC)	yodC	protein YodC
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (33/25);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACT  >  NC_000913/2028130‑2028397                                                                                                                                       |                                                                                                                                      ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCg                                                                                                                                  >  1:426257/1‑140 (MQ=255) ttGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCg                                                                                                                                  >  1:1032003/1‑140 (MQ=255)  tGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGt                                                                                                                                 <  1:744862/140‑1 (MQ=255)  tGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGt                                                                                                                                 <  1:48382/140‑1 (MQ=255)  tGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGt                                                                                                                                 >  1:464343/1‑140 (MQ=255)   gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTc                                                                                                                                <  2:426257/140‑1 (MQ=255)     aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCAt                                                                                                                              >  1:137510/1‑140 (MQ=255)        cATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATAcc                                                                                                                           <  1:827981/140‑1 (MQ=255)              cacaCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGac                                                                                                                     >  2:925676/1‑140 (MQ=255)               acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGaca                                                                                                                    <  1:138937/140‑1 (MQ=255)               acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGaca                                                                                                                    <  2:161224/140‑1 (MQ=255)                cacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGacac                                                                                                                   <  2:486491/140‑1 (MQ=255)                 acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCaa                                                                                                                         <  1:780888/133‑1 (MQ=255)                 acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCaa                                                                                                                         >  2:780888/1‑133 (MQ=255)                 acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACt                                                                                                                  >  1:282555/1‑140 (MQ=255)                 acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACt                                                                                                                  >  2:1057104/1‑140 (MQ=255)                 acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACt                                                                                                                  >  2:253593/1‑140 (MQ=255)                       ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCa                                                                                                            <  1:594759/140‑1 (MQ=255)                            tgtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCg                                                                                                       <  2:309545/140‑1 (MQ=255)                                 tCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGa                                                                                                  <  1:925676/140‑1 (MQ=255)                                     gCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTAt                                                                                              <  2:242257/140‑1 (MQ=255)                                     gCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTAt                                                                                              >  2:751150/1‑140 (MQ=255)                                         gATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAg                                                                                          <  2:282555/140‑1 (MQ=255)                                                    tACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAatcatc                                                                               >  1:332072/1‑140 (MQ=255)                                                     aCCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATcc                                                                              >  2:754997/1‑140 (MQ=255)                                                     aCCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATcc                                                                              >  2:330313/1‑140 (MQ=255)                                                      cccTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCcg                                                                             >  1:96671/1‑140 (MQ=255)                                                       ccTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCcgc                                                                            >  1:521776/1‑140 (MQ=255)                                                       ccTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCcgc                                                                            >  1:932290/1‑140 (MQ=255)                                                         tttCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCgg                                                                          >  2:608381/1‑140 (MQ=255)                                                         tttCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCgg                                                                          >  2:815248/1‑140 (MQ=255)                                                          ttCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCggg                                                                         >  2:115492/1‑140 (MQ=255)                                                             acacTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccg                                                                      <  1:643629/140‑1 (MQ=255)                                                                 ttcttcCGCAGAACGACTCCCCTCCCCCGGAATAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCcgc                                                                            <  2:540408/130‑1 (MQ=255)                                                                 ttcttcCGCAGAACGACTCCCCTCCCCCGGAATAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCcgc                                                                            >  1:540408/1‑130 (MQ=255)                                                                 ttcttcCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccg                                                                      >  2:628091/1‑136 (MQ=255)                                                                 ttcttcCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGgccg                                                                      <  1:628091/136‑1 (MQ=255)                                                                 ttcttcCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCct                                                                  >  1:735796/1‑140 (MQ=255)                                                                    ttcCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCtt                                                               >  2:866741/1‑140 (MQ=255)                                                                    ttcCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCAGgccg                                                                      <  2:1059662/133‑1 (MQ=255)                                                                    ttcCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCAGgccg                                                                      >  1:1059662/1‑133 (MQ=255)                                                                       cGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTaa                                                            <  2:332072/140‑1 (MQ=255)                                                                         cAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAAcc                                                          >  1:221245/1‑140 (MQ=255)                                                                         cAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAAcc                                                          >  2:480688/1‑140 (MQ=255)                                                                            aaCGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTa                                                       >  1:1139649/1‑140 (MQ=255)                                                                                 ctcccctcccCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACATTCATCCGCGGGCCGCCCTCTTTAACCGTAACTTc                                                  <  2:1032003/140‑1 (MQ=255)                                                                                    ccctccccCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTc                                               <  1:754997/140‑1 (MQ=255)                                                                                      ctccccCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCAc                                             >  2:180078/1‑140 (MQ=255)                                                                                      ctccccCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCAc                                             >  1:687335/1‑140 (MQ=255)                                                                                             ggAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAGCGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCAt                                      <  1:1081332/140‑1 (MQ=255)                                                                                              gAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                     <  2:96671/140‑1 (MQ=255)                                                                                              gAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                     <  2:1139649/140‑1 (MQ=255)                                                                                              gAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                     <  1:751150/140‑1 (MQ=255)                                                                                              gAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATa                                     <  1:608381/140‑1 (MQ=255)                                                                                               aaCAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATaa                                    >  2:142829/1‑140 (MQ=255)                                                                                                   aGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCt                                >  2:883556/1‑140 (MQ=255)                                                                                                            tCATGAAAAGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGc                       <  1:330313/140‑1 (MQ=255)                                                                                                                  aaaGCTTCCCGCTTGACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCtt                 >  2:959331/1‑140 (MQ=255)                                                                                                                  aaaGCTTCCCGCTTGACCCC‑CAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCttt                >  2:592541/1‑140 (MQ=255)                                                                                                                                gACCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAAc   <  1:253840/140‑1 (MQ=255)                                                                                                                                 aCCCCGCAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACt  >  2:1026734/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                      TTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTGACCCCGTAACCGTCATACCAACGACACTCAACCATACCGCTGGAGTATCCAGTGACAATCATCCGCGGGCCGCCCTCTTTAACCGTAACTTCCTCACTAACCATAAAGCTCATACTCGCCTCCTTTTTTCGAGTGAAACT  >  NC_000913/2028130‑2028397

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