Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,771,316	T→A	13.0%	Y82* (TAT→TAA)	ydiM →	putative MFS transporter, membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,771,316	0	T	A	13.0%	97.2 / 5.4	38	Y82* (TAT→TAA)	ydiM	putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base A (2/3);  total (18/20)
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

GAGACACTTTGGCAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGA  >  NC_000913/1771170‑1771458                                                                                                                                                   |                                                                                                                                               gagaCTCTTTGGCAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGGCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCtata                                                                                                                                               <  1:293949/148‑1 (MQ=255) gagaCACTTTGGCAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCtat                                                                                                                                                <  7:64282/147‑1 (MQ=255)         ttGGCAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGCCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCcttc                                                                                                                                       >  3:110152/1‑148 (MQ=255)          tGGCAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCcttc                                                                                                                                       <  3:143038/147‑1 (MQ=255)             cAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGCCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCtttt                                                                                                                                   <  4:110152/148‑1 (MQ=255)              aGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCttttt                                                                                                                                  <  3:241800/148‑1 (MQ=255)               gACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTTTCATGCTCGGGAAGTGCTGCTATATGGCCTTCtttttt                                                                                                                                 >  4:297789/1‑148 (MQ=255)                 cTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTgg                                                                                                                               >  3:98737/1‑148 (MQ=255)                       ccGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGGTATATGGCCTTCTTTTTTGGCATCCt                                                                                                                         >  6:233241/1‑148 (MQ=255)                              tgtCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCAGATCGCTTTGGTCTCGGCCCTTTTATCATGCTCGGGATGTGCTGCtata                                                                                                                                               <  5:246111/119‑1 (MQ=255)                               gtCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATAATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGATGCTatat                                                                                                                                              >  6:246111/1‑119 (MQ=255)                                   cGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCtat                                                                                                                                                >  3:71318/1‑113 (MQ=255)                                   cGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCtat                                                                                                                                                <  4:71318/113‑1 (MQ=255)                                   cGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTAAATGGCCTTCTTTTTTGGCAt                                                                                                                            <  1:90151/133‑1 (MQ=255)                                   cGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTAAATGGCCTTCTTTTTTGGCAt                                                                                                                            <  1:85733/133‑1 (MQ=255)                                   cGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTAAATGGCCTTCTTTTTTGGCAt                                                                                                                            >  2:85733/1‑133 (MQ=255)                                   cGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTAAATGGCCTTCTTTTTTGGCAt                                                                                                                            >  2:90151/1‑133 (MQ=255)                                    gATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATaa                                                                                                             <  1:14697/147‑1 (MQ=255)                                        gTTATCTCATCGCTGGGCCTTGTTCGATTAAGTGTCTTGATTTATGCAGGATTTTTCTCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTTCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATTAcatca                                                                                                        <  1:312003/148‑1 (MQ=255)                                         ttATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAAcatcat                                                                                                       <  2:175889/148‑1 (MQ=255)                                             ctcATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATc                                                                                                   >  1:249269/1‑148 (MQ=255)                                               cATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGc                                                                                                 >  7:297741/1‑148 (MQ=255)                                                           ttGGTCGATTAAGTGTCTTGCTTTTTGAAGGATTTTTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCtat                                                                                                                                                <  6:120696/89‑1 (MQ=255)                                                            tGGTCGATTAAGTGTCTTGCTTTATGCAGGATTATTATCCGATAGCTTTGGTCTCAGCACTATTATCATGCTCTGGATGTGCTGCtata                                                                                                                                               >  5:120696/1‑89 (MQ=255)                                                              gTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATc                                                                                                                           <  2:304481/107‑1 (MQ=255)                                                              gTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATc                                                                                                                           >  1:304481/1‑107 (MQ=255)                                                               tCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTgg                                                                                                                               >  3:90093/1‑102 (MQ=255)                                                               tCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTgg                                                                                                                               <  4:90093/102‑1 (MQ=255)                                                                  aTTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTg                                                                              >  3:35977/1‑148 (MQ=255)                                                                        tgtCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGa                                                                        <  4:98737/148‑1 (MQ=255)                                                                           cTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTTTCATGCTCGGGATGTGATGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATg                                                                     >  3:209507/1‑148 (MQ=255)                                                                               cTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCa                                                                  >  1:227179/1‑147 (MQ=255)                                                                                      gaGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTAAATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGtt                                                           <  6:181566/146‑1 (MQ=255)                                                                                      cAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATAATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGtt                                                           <  4:81246/147‑1 (MQ=255)                                                                                           ttattaTCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCg                                                     <  4:35977/148‑1 (MQ=255)                                                                                            tattatCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGa                                                    >  6:171695/1‑148 (MQ=255)                                                                                             attatCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCg                                                     <  7:244029/146‑1 (MQ=255)                                                                                             attatCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCg                                                     >  8:244029/1‑146 (MQ=255)                                                                                                atCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGc                                                 <  2:249269/147‑1 (MQ=255)                                                                                                atCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGc                                                 <  2:227179/147‑1 (MQ=255)                                                                                                                    gccCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGAt                            <  7:93479/148‑1 (MQ=255)                                                                                                                       cTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATgg                          >  2:259369/1‑147 (MQ=255)                                                                                                                       cTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATgg                          >  2:228702/1‑147 (MQ=255)                                                                                                                                tGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTcc                <  5:171695/148‑1 (MQ=255)                                                                                                                                           tgctgcTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCt      <  1:259369/147‑1 (MQ=255)                                                                                                                                                ctATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGa  >  5:42798/1‑146 (MQ=255)                                                                                                                                                   |                                                                                                                                               GAGACACTTTGGCAGACTAATGCCGCGGGTGTCTCGATAGTTATCTCATCGCTGGGCATTGGTCGATTAAGTGTCTTGCTTTTTGCAGGATTATTATCCGATCGCTTTGGTCGCCGCCCTTTTATCATGCTCGGGATGTGCTGCTATATGGCCTTCTTTTTTGGCATCCTGCAGACCAATAACATCATTATCGCTTATGTTTTTGGCTTTCTGGCGGGAATGGCAAACAGTTTTCTCGATGCAGGCACTTATCCCAGTTTGATGGAAGCTTTTCCACGCTCACCTGGGA  >  NC_000913/1771170‑1771458

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