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	93.9%	102.8 / ‑0.7	34	Q62* (CAG→TAG)	ytfI	protein YtfI
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (19/12);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

TATAGTGAGATATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAG  >  NC_000913/4437652‑4438110                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             tataGTGAGATATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATaaa                                                                                                                                                                                                                   >  2:601810/1‑250 (MQ=255)  ataGTGAGATATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACACTAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCCTCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATaaa                                                                                                                                                                                                                   >  1:475778/1‑249 (MQ=255)  ataGTGAGATATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATaaa                                                                                                                                                                                                                   >  2:145003/1‑249 (MQ=255)            atTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAg                                                                                                                                                                                                       <  1:384801/251‑1 (MQ=255)                   ggCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCttt                                                                                                                                                                                                >  1:704293/1‑251 (MQ=255)                        cATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTtcatc                                                                                                                                                                                            >  2:414941/1‑250 (MQ=255)                         aTAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTtcatca                                                                                                                                                                                           >  1:409573/1‑250 (MQ=255)                           aaaGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCaaa                                                                                                                                                                                         <  2:601672/250‑1 (MQ=255)                            aaGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCaaaa                                                                                                                                                                                        >  1:245593/1‑250 (MQ=255)                                        aaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGAt                                                                                                                                                                                                                      >  1:614736/1‑208 (MQ=255)                                        aaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGAt                                                                                                                                                                                                                      <  2:614736/208‑1 (MQ=255)                                         aaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAAATTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAACCAATTTTTTTAACAGCATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAAGTATGATAGCCTGTTAACTTGGAGGGAGTTCTTTCGGTGCGAAACTTAGACTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAAcggctca                                                                                                                                                                          <  2:261229/251‑7 (MQ=255)                                                              cccAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAATTGAAACTGTTCACGCCATGGATAAGTTTCTTAt                                                                                                                                                     >  2:572244/1‑251 (MQ=255)                                                                cAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATc                                                                                                                                                    >  1:757873/1‑250 (MQ=255)                                                                             aataatTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTa                                                                                                                                       >  1:668817/1‑250 (MQ=255)                                                                                      aTTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACgg                                                                                                                             >  1:456541/1‑251 (MQ=255)                                                                                                          tGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACa                                                                                                          <  1:414941/250‑1 (MQ=255)                                                                                                                    tCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAACATGTTAACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGt                                                                                                 >  2:307819/1‑249 (MQ=255)                                                                                                                           aaaaaaaACATTCAAACAATTTTTTTAACAACATCCAACATGCATTTAATAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGTGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTTCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTATTTGAttt                                                                                        <  2:661775/250‑1 (MQ=255)                                                                                                                              aaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGa                                                                                      >  2:19528/1‑250 (MQ=255)                                                                                                                               aaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAAt                                                                                    <  1:134099/251‑1 (MQ=255)                                                                                                                                       cAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCTAAAGATGATAGCCTGTTAAATTTGAGGGAGTTCTTTCTATGCGCAAATTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAAGGAAACTGTTCACGACATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATGGCGAGGTAGAGGATTTGAATACgaaga                                                                             <  2:32807/250‑1 (MQ=255)                                                                                                                                             aTTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATAccc                                                                      >  2:703123/1‑251 (MQ=255)                                                                                                                                               ttttttAACAACATCCATCATGCATTTACGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGAGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGTTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCAt                                                                    >  1:598197/1‑251 (MQ=255)                                                                                                                                                                  atGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAAc                                                 >  1:810454/1‑251 (MQ=255)                                                                                                                                                                      aTTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTACCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACtcat                                             <  2:765583/251‑1 (MQ=255)                                                                                                                                                                           agtaaaaaGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATg                                        <  2:395531/248‑1 (MQ=255)                                                                                                                                                                                     ttGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCaaa                               >  2:575571/1‑250 (MQ=255)                                                                                                                                                                                                   atgatgggcTGTTAAATTTGATGGAGTTTTTTCGGTGTGCAACTGAGATTGATAAACGGGGGTATCAGGTGCTTTCAGTAAAAAAAGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGAGAGGTAGTGTATTTGAATACTAAGAAATACCTATAACAATAAACGAATACAACTCATTCATTGCTTTCAAAAACCGGGAATGTcc                 <  2:475778/242‑1 (MQ=255)                                                                                                                                                                                                              ttAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGATTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACtc                                               >  2:391723/1‑209 (MQ=255)                                                                                                                                                                                                              ttAACTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGATTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACtc                                               <  1:391723/209‑1 (MQ=255)                                                                                                                                                                                                                 aCTTTGAGGGAGTTCTTTCGTTGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAg  <  1:575571/251‑1 (MQ=255)                                                                                                                                                                                                                                      tGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAatat                                                                                                                   >  2:352018/1‑117 (MQ=255)                                                                                                                                                                                                                                      tGCGCAACTTAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAatat                                                                                                                   <  1:352018/117‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                             TATAGTGAGATATTTAATGGCGACATAAAGTAACCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCCTTAAAGATAACGGTGCCGAATATACAATGACATTGCGAGGTAGTGGATTTGAATACGAAGAAATACCCATAACAATAAATGAATATAACTCATTCATGGATTTCAAAAACCGTGAATTTCCGTTAGAACAAAACAG  >  NC_000913/4437652‑4438110

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