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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	696,793	0	T	A	10.6%	123.3 / 2.6	44	noncoding (38/75 nt)	glnW	tRNA‑Gln
Reads supporting (aligned to +/- strand):  ref base T (15/24);  new base A (3/2);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGG  >  NC_000913/696644‑696901                                                                                                                                                      |                                                                                                             tatTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                                                              <  2:110554/150‑1 (MQ=255) tatTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                                                              >  1:1338972/1‑150 (MQ=255)  tatggATACATTCAGATTTTGGTTCCTACGAATGGATTCGATCATTTTACCCCTTCATTCTGAGTGTTGTGCTCTAACCAACTGAGCTACGTAGCGAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTa                                                                                                             <  2:487099/146‑1 (MQ=25)   ttCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                                                            >  1:487099/1‑150 (MQ=255)       aTACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTTGGTTACCTGGATTCGAACCAGGGAATGCCGGTATCaaa                                                                                                        <  2:542404/150‑1 (MQ=255)          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                                                    <  2:740773/150‑1 (MQ=255)                                                           ttATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCt                                                    <  1:1030423/150‑1 (MQ=255)                                                           ttATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCt                                                    <  1:1246604/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)                                                                         ctctAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                                                              <  2:1248271/78‑1 (MQ=255)                                                                         ctctAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGt                                                                                                              >  1:1248271/1‑78 (MQ=255)                                                                              aCCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTg                                 <  1:945546/150‑1 (MQ=255)                                                                                cAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTggg                               <  1:540995/150‑1 (MQ=255)                                                                                     gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGTATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTAcc                          <  1:749977/150‑1 (MQ=255)                                                                                     gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                                                            >  1:719140/1‑68 (MQ=255)                                                                                     gAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTAt                                                                                                            <  2:719140/68‑1 (MQ=255)                                                                                                             tcGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCgg  <  1:3766/150‑1 (MQ=255)                                                                                                                             gaggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                                                 >  1:220910/3‑87 (MQ=11)                                                                                                                             gaggATTCGAACCTCGGAATGCCGGAATCAGAATCCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTAc                                                 <  2:220910/85‑1 (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)                                                                                                                                                      |                                                                                                             TATTCGATACATTCAGAATTTGGTGGCTACGACGGGATTCGAACCTGTGACCCCATCATTATGAGTGATGTGCTCTAACCAACTGAGCTACGTAGCCAGATTGTTTCTTCGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGGTATCAAAAACCGGTGCCTTACCGCTTGGCGATACCCCAATAACCGGGCGGTGAACCGCTTACTCGAAGAAGATGGCTGGGGTACCTGGATTCGAACCAGGGAATGCCGG  >  NC_000913/696644‑696901

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