Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,841,383	C→T	G119G (GGG→GGA)	nepI ←	purine ribonucleoside exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,841,383	0	C	T	96.3%	180.6 / ‑0.6	54	coding (357/1191 nt)	nepI	purine ribonucleoside exporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/28);  minor base . (1/1);  total (25/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 = 7.86e-01

CTGCCCAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGCCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAA  >  NC_000913/3841246‑3841516                                                                                                                                          |                                                                                                                                      cTGCCCAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGtcc                                                                                                                                     >  2:710893/1‑140 (MQ=255)   gCCCAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCa                                                                                                                                   >  2:504683/1‑140 (MQ=255)   gCCCAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCa                                                                                                                                   >  2:53230/1‑140 (MQ=255)    cccAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAg                                                                                                                                  <  1:381621/140‑1 (MQ=255)    cccAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAg                                                                                                                                  <  1:516679/140‑1 (MQ=255)     ccAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGa                                                                                                                                 <  2:148642/140‑1 (MQ=255)         cGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAgg                                                                                                                             >  2:150070/1‑140 (MQ=255)            cgcgGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCAc                                                                                                                          >  2:902843/1‑140 (MQ=255)              cgGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGa                                                                                                                        <  2:9581/140‑1 (MQ=255)               ggCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGAc                                                                                                                       <  1:924545/140‑1 (MQ=255)               ggCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGAc                                                                                                                       <  2:440323/140‑1 (MQ=255)                   aTCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGAt                                                                                                                   <  2:227381/140‑1 (MQ=255)                         aGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAg                                                                                                             >  1:632339/1‑140 (MQ=255)                           cgcAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCa                                                                                                           >  1:513684/1‑140 (MQ=255)                           cgcAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCa                                                                                                           <  2:227647/140‑1 (MQ=255)                                 aGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGa                                                                                                     <  1:150070/140‑1 (MQ=255)                                  gAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGaa                                                                                                    <  1:664566/140‑1 (MQ=255)                                  gAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGaa                                                                                                    <  2:168108/140‑1 (MQ=255)                                    aaCCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATg                                                                                                  <  2:753814/140‑1 (MQ=255)                                     aCCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGa                                                                                                 <  2:480907/140‑1 (MQ=255)                                       cgcgCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGt                                                                                               <  1:504683/140‑1 (MQ=255)                                         cgcCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTa                                                                                             <  1:710893/140‑1 (MQ=255)                                          gcCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAg                                                                                            <  2:470504/140‑1 (MQ=255)                                           ccGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGc                                                                                           >  1:953245/1‑140 (MQ=255)                                                 aTCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGa                                                                                     <  1:53230/140‑1 (MQ=255)                                                      cGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCa                                                                                <  1:979825/140‑1 (MQ=255)                                                         cAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCagca                                                                             >  1:703866/1‑140 (MQ=255)                                                          aGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCagcaa                                                                            >  1:272056/1‑140 (MQ=255)                                                           gcgcCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCagcaag                                                                           >  2:671223/1‑140 (MQ=255)                                                            cgcCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCagcaagc                                                                          >  1:83767/1‑140 (MQ=255)                                                                 tCGACACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGTAGgag                                                                     <  2:196515/140‑1 (MQ=255)                                                                         gTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAg                                                             >  2:339482/1‑140 (MQ=255)                                                                               ggcggcACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAG‑‑TCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCa                                                           >  1:440575/1‑134 (MQ=255)                                                                               ggcggcACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAG‑‑TCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCa                                                           <  2:440575/134‑1 (MQ=255)                                                                                 cggcACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACgg                                                     >  1:643003/1‑140 (MQ=255)                                                                                        aGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGcaaacaa                                              <  2:643003/140‑1 (MQ=255)                                                                                         gACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGcaaacaaa                                             >  1:375326/1‑140 (MQ=255)                                                                                              aTGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATaac                                        >  2:974092/1‑140 (MQ=255)                                                                                              aTGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATaac                                        <  2:632339/140‑1 (MQ=255)                                                                                                    aGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTa                                  >  1:836599/1‑140 (MQ=255)                                                                                                        aCGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAgcgg                              >  2:868066/1‑140 (MQ=255)                                                                                                         cGCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAgcggc                             <  1:974092/140‑1 (MQ=255)                                                                                                          gCCGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAgcggcg                            >  1:283638/1‑140 (MQ=255)                                                                                                           ccGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATaa                                         >  2:915774/1‑126 (MQ=255)                                                                                                           ccGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATaa                                         <  1:915774/126‑1 (MQ=255)                                                                                                            cGACATCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGcggcggg                          >  2:234276/1‑139 (MQ=255)                                                                                                                aTCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGt                      <  2:129381/140‑1 (MQ=255)                                                                                                                aTCGCCCAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGt                      <  2:953245/140‑1 (MQ=255)                                                                                                                      cAGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTg                >  2:762713/1‑140 (MQ=255)                                                                                                                       aGAACCCGCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGa               <  2:513684/140‑1 (MQ=255)                                                                                                                              gCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCt        >  1:609918/1‑140 (MQ=255)                                                                                                                              gCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCt        <  1:902843/140‑1 (MQ=255)                                                                                                                              gCCCAGCGCCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCt        <  2:267887/140‑1 (MQ=255)                                                                                                                                   gcgcCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGta   >  1:682767/1‑140 (MQ=255)                                                                                                                                    cgcCAGTCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTCGCCTTAATTGTCTAGGgaa  >  2:179803/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                      CTGCCCAACGGCGCGGCAATCACCAGCGCAATAGAAACCGCGCCGAAGATCACCGACAGCGCCTTCGGCACCGTACGCGGCGGCACCAGACGCATGGTCAGCGACGCCGACATCGCCCAGAACCCGCCCAGCGCCAGCCCCAGACAGGCACGACCGATTAAAAGCAAACTGAATGAGTTAGCAAAGGAAACCAGCAAGCAGGAGAGCGTCAGCAAAACGGCAAACAAAATAACAACGTAGCGGCGGTCAGTAGCCTGAATTGTCTGGGTAA  >  NC_000913/3841246‑3841516

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: ‑