Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,437,890	C→T	Q62* (CAG→TAG)	ytfI →	protein YtfI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,437,890	0	C	T	87.5%	72.8 / 2.3	33	Q62* (CAG→TAG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (14/14);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

ATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGG  >  NC_000913/4437663‑4438126                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              atTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCa                                                                                                                                                                                                                        <  1:430005/250‑1 (MQ=255) atTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTCTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTAGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAg                                                                                                                                                                                                                       <  1:151377/251‑1 (MQ=255)     aaTGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCAAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAAATTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATAAGGTGc                                                                                                                                                                                                                   >  1:467289/1‑251 (MQ=255)      aTGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCt                                                                                                                                                                                                                  >  1:98591/1‑251 (MQ=255)                      aCCAAATAAAATAAGGTTGTCATATGTTTCCCAGGATCAGACACAATAATTGTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCATACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑TTGGTTATCAGGTGCTTTCATCAAAAAATGaa                                                                                                                                                                                                  <  2:350805/251‑1 (MQ=255)                                    accgacgcgtttacgacttggATCAGACAAAAAAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACAATCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTCGAGGGAGTTCTTTCGTTGCGCAACTCAGATGGATAAAA‑GTGGTAATCAGGTGCTTTCATCAAAAAATGAAACTGTTCCCGCCAt                                                                                                                                                                                    <  2:115493/232‑1 (MQ=255)                                         tCATCTGTTACACATGATCAGCCACAATAATTTAATTGGTGCGTTGGAGTTATTTGTAAAGTCTTCGTATACAAACAAACATTCAAACAGTTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAGTTATGCTAACCTATTCACTTTGAGGGATATCTTTCGTTGCTCAACTTAGATTGATAAAA‑GTTGTTATCCGTTGCTTTCTTCAAAAACTGAATCTGTTCACGCCATGGAt                                                                                                                                                                                >  2:25987/1‑250 (MQ=255)                                                       ggATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTCGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCAt                                                                                                                                                                                    >  1:115493/1‑232 (MQ=255)                                                           cAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCa                                                                                                                                                             <  1:636104/251‑1 (MQ=255)                                                                           aTTGGTGCGGTGGAGTTATTTGTAAAGTCTACTTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCTAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTTGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACgg                                                                                                                                             <  2:509758/251‑1 (MQ=255)                                                                                          ttATTTGTAAAGTCTTCGTATCCAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAACAGATTGGATTTCGAATTATGATAGCCTGTTACCTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGATATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAAt                                                                                                                              >  2:670004/1‑251 (MQ=255)                                                                                             tttGTAAAGTCTTCGTATAAAAAAACACATTCAAACACTTTTTTTAACAAACTCCATCATGCATTTAAGAAAAAAGATTGGATTTCCAATTATGATAGCCTGTTACATTTGAGGGGGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTCTCAGGTGCTTTCATCAAAAAATGTAACTGTTCACGCCCTGGATACGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATCTACAATGAc                                                                                                                           >  2:298853/1‑251 (MQ=255)                                                                                                      tCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATGAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGACATGGATAAGTTTCTTATCGGTTTCAGCATTAAAGATAACGGTGCCGAATATACAATGACATTGCGAgg                                                                                                                  >  1:679542/1‑251 (MQ=255)                                                                                                               aCAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTGCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACGGTTCACGCCATGGATAAGTTTCTTAGCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGAt                                                                                                          <  1:426603/250‑1 (MQ=255)                                                                                                                            cAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCAAGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACgaagg                                                                                             >  1:580078/1‑249 (MQ=255)                                                                                                                             aaaCAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACgaagaa                                                                                            >  1:580828/1‑250 (MQ=255)                                                                                                                                  aTTTTTTTAACAACAGCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTGCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATAccc                                                                                      <  1:670004/251‑1 (MQ=255)                                                                                                                                      ttttAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATa                                                                                   <  2:123187/250‑1 (MQ=255)                                                                                                                                         tAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACa                                                                                <  1:404346/250‑1 (MQ=255)                                                                                                                                          aacaacATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAAt                                                                              >  1:332259/1‑251 (MQ=255)                                                                                                                                               catccatcATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACaa                                                                                                                               <  1:5290/197‑1 (MQ=255)                                                                                                                                               catccatcATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACaa                                                                                                                               >  2:5290/1‑197 (MQ=255)                                                                                                                                                atccatcATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTTTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTTCGAGGTAGTGGATTTGAATCCGAAGAAATACCCATATCAATAAATGa                                                                        <  1:298853/251‑1 (MQ=255)                                                                                                                                                           aTTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACATTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GGGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGAAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACtaat                                                             >  1:573564/1‑248 (MQ=255)                                                                                                                                                                   aaaaaGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGTAGTTCTTTCGGTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTTCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTTATGGAt                                                     <  2:81496/251‑1 (MQ=255)                                                                                                                                                                        gATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATGAAA‑ATGGTTATACGGAGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGGTTCAGCCTTAAAGCTAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGCTTTGAATACGAAGAAATACCCATAACAATAAATGCATATAACTCATTCATGGATTTCa                                                 >  1:435999/1‑250 (MQ=255)                                                                                                                                                                                  cGAATTATGATAGCCTGTTAACTTTGAGGGATGTCTTTCGGTGCGCAACTCAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTGCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGa                                       <  2:332259/250‑1 (MQ=255)                                                                                                                                                                                                       cTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGAGAATTTCCGTTAGAACAAAAc                    >  2:502761/1‑248 (MQ=255)                                                                                                                                                                                                          tGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACg               >  2:510170/1‑250 (MQ=255)                                                                                                                                                                                                           gAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGAt             >  1:584353/1‑251 (MQ=255)                                                                                                                                                                                                                aGTTCTTTCGTTGCGCAAGTGAAATTGATAAAAAGTG‑TTATCAGTTGTCTTCATCAAAAAAAGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTGCAGACCGAAAGATAACGGTGCCGAATAAACAATTACATTGCAAGTTAGAGGATTAGAATAGGAAAAAAAACCCATTACATTAAATGAATATATCACATGCATGGATTTCAAAAACCGAGAATTGCGGTTGGAACAAAACACACGATTGTAt        <  2:435999/251‑1 (MQ=255)                                                                                                                                                                                                                      ttctttGCGCAACTTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTg   <  2:98591/246‑1 (MQ=255)                                                                                                                                                                                                                       tCGTTGCGCAACGTAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTCCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTgg  <  2:580828/250‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                              ATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAA‑GTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAGACGATTGTATGCCTGG  >  NC_000913/4437663‑4438126

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