Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	696,801	A→T	10.3%	noncoding (30/75 nt)	glnW ←	tRNA‑Gln

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,801	0	A	T	10.3%	109.4 / 3.5	39	noncoding (30/75 nt)	glnW	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base T (2/2);  total (16/23)
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.46e-01

CATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGG  >  NC_000913/696653‑696901                                                                                                                                                     |                                                                                                     cATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAc                                                                                                     >  2:936637/1‑150 (MQ=255)  aTTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                                                                    >  2:1030423/1‑150 (MQ=255)  aTTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAAcc                                                                                                    >  1:1295679/1‑150 (MQ=255)     cAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                                                 >  2:945546/1‑150 (MQ=255)     cAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGt                                                                                                 <  2:1338972/150‑1 (MQ=255)        aaTTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGcc                                                                                              <  1:531710/150‑1 (MQ=255)          tttGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCtt                                                                                            >  2:1218793/1‑150 (MQ=255)            tggtggCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAc                                                                                          <  1:72455/150‑1 (MQ=255)             ggtggCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTAcc                                                                                         <  1:936637/150‑1 (MQ=255)                  cTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTg                                                                                    >  2:3766/1‑150 (MQ=255)                   tACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTgg                                                                                   <  2:66002/150‑1 (MQ=255)                     cgacgGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCg                                                                                 <  1:1218793/150‑1 (MQ=255)                          ggATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATAcc                                                                            >  2:540995/1‑150 (MQ=255)                           gATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGAATCAAAAACCGGTGCATTACCGCTTGGAGATAccc                                                                           >  2:749977/1‑150 (MQ=255)                              tCGAACCTGTGACCCCATCATTATGAGTCATCTGCTCTAACCAACTGAGCTACGTATCCTTATTTTTTCTTAGATGGCTGTGGTACATCGATTCGAACCACGTAATTTCGGTATCAATACGCGGTTCCTTACCGCTTGGCGATACCCCaa                                                                        <  2:113967/150‑1 (MQ=18)                               cGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAAt                                                                       <  1:636173/150‑1 (MQ=255)                                    cTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCaagaggac                                                                  >  1:774827/1‑144 (MQ=255)                                        gACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCATAACGCAAGACCAAACCACAATaagccagcc                                                              >  1:1027077/1‑143 (MQ=25)                                          ccccATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGt                                                            <  1:1199957/150‑1 (MQ=255)                                            ccATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGa                                                          >  1:740773/1‑150 (MQ=255)                                                  ttATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCt                                                    <  1:1246604/150‑1 (MQ=255)                                                  ttATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCt                                                    <  1:1030423/150‑1 (MQ=255)                                                  ttATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCt                                                    <  2:740773/150‑1 (MQ=255)                                                   tATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCtt                                                   <  2:148704/150‑1 (MQ=255)                                                    aTGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTa                                                  >  1:94337/1‑150 (MQ=255)                                                       aGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTc                                               <  2:94337/150‑1 (MQ=255)                                                         tGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATAc                                                                             >  2:266480/1‑118 (MQ=255)                                                         tGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATAc                                                                             <  1:266480/118‑1 (MQ=255)                                                         tGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTATTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCga                                             <  1:1258245/150‑1 (MQ=255)                                                                     aCCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTg                                 <  1:945546/150‑1 (MQ=255)                                                                       cAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTggg                               <  1:540995/150‑1 (MQ=255)                                                                            gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGTATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTAcc                          <  1:749977/150‑1 (MQ=255)                                                                                                    tcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCgg  <  1:3766/150‑1 (MQ=255)                                                                                                                    gaggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                                                 <  2:220910/85‑1 (MQ=11)                                                                                                                    gaggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                                                 >  1:220910/3‑87 (MQ=11)                                                                                                                     tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATa                                                                      <  2:495888/65‑1 (MQ=35)                                                                                                                     tGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATa                                                                      >  1:495888/1‑65 (MQ=35)                                                                                                                                         ccGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc                                                       <  1:341216/60‑1 (MQ=14)                                                                                                                                         ccGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAAcc                                                       >  2:341216/1‑60 (MQ=14)                                                                                                                                                     |                                                                                                     CATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGG  >  NC_000913/696653‑696901

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