Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,331	A→G	noncoding (71/87 nt)	leuQ ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,331	0	A	G	78.5%	126.4 / 35.5	67	noncoding (71/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (6/8);  new base G (19/32);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 7.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAA  >  NC_000913/4606285‑4606361                                                |                               gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCtt                             <  1:84463/50‑1 (MQ=255) gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACtt                             <  1:1602411/50‑1 (MQ=255)  cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            <  1:5182362/50‑1 (MQ=255)  cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            <  1:1546389/50‑1 (MQ=255)  cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTg                            <  1:5903644/50‑1 (MQ=255)   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            >  1:7895814/1‑49 (MQ=255)    cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGaa                          <  1:6108045/50‑1 (MQ=255)     aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAc                         >  1:4234361/1‑50 (MQ=255)      gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcc                        >  1:6807233/1‑50 (MQ=255)      gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcc                        <  1:8226641/50‑1 (MQ=255)       tGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccc                       <  1:1229155/50‑1 (MQ=255)         cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcccc                      <  1:2056524/49‑1 (MQ=255)         cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcccc                      <  1:705605/49‑1 (MQ=255)         cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGG‑CTTGAAccccc                     <  1:8607483/49‑1 (MQ=25)          gAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGCTTGAAccccc                     <  1:2426729/50‑1 (MQ=25)          gAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCa                    <  1:8102413/50‑1 (MQ=255)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGGACCCCCAc                   <  1:7756923/50‑1 (MQ=25)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:4547081/50‑1 (MQ=255)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:7253532/50‑1 (MQ=255)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:2604543/50‑1 (MQ=255)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:190857/50‑1 (MQ=255)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:1655304/50‑1 (MQ=255)           aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:1276424/50‑1 (MQ=255)            gCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCAc                   <  1:1749144/50‑1 (MQ=25)            gCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACg                  >  1:684465/1‑50 (MQ=255)            gCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                  >  1:5924931/1‑50 (MQ=255)            gCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                  >  1:1506881/1‑50 (MQ=255)             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                 >  1:7321508/1‑50 (MQ=255)             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                 <  1:3904005/50‑1 (MQ=255)             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                 <  1:5467464/50‑1 (MQ=255)              aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                <  1:6481906/50‑1 (MQ=255)              aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                >  1:998057/1‑50 (MQ=255)              aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTc                <  1:6628106/50‑1 (MQ=255)               atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTcc               >  1:2040311/1‑50 (MQ=255)               atatCAACATGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTcc               <  1:2924895/50‑1 (MQ=255)                tatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCg              >  1:975314/1‑50 (MQ=255)                tatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCg              <  1:6606403/50‑1 (MQ=255)                  tCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGt             >  1:4201247/1‑49 (MQ=255)                  tCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTa            <  1:6504065/50‑1 (MQ=255)                   cAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGt             >  1:7763023/1‑48 (MQ=255)                    aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTaa           <  1:1968541/49‑1 (MQ=255)                    aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          <  1:1509428/50‑1 (MQ=255)                    aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          <  1:3252618/50‑1 (MQ=255)                    aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  1:3361223/1‑50 (MQ=255)                    aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          <  1:1884241/50‑1 (MQ=255)                    aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          <  1:1553460/50‑1 (MQ=255)                    aaCGTGGTTTTTGGAGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  1:6821053/1‑50 (MQ=17)                     aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAgg         <  1:3547426/50‑1 (MQ=255)                     aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  1:574026/1‑49 (MQ=255)                     aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  1:7279418/1‑49 (MQ=255)                      cGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGa        >  1:842729/1‑50 (MQ=255)                      cGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGa        >  1:7024954/1‑50 (MQ=255)                       gTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGac       <  1:6330020/50‑1 (MQ=255)                       atGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGac       <  1:228480/49‑1 (MQ=255)                        tGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGaca      <  1:815373/50‑1 (MQ=255)                        tGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGaca      >  1:1514898/1‑50 (MQ=255)                        tGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGaca      >  1:3701543/1‑50 (MQ=255)                          gTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac     >  1:5939416/1‑49 (MQ=255)                          gTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACt    <  1:872433/50‑1 (MQ=255)                           tttttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTa   <  1:380777/50‑1 (MQ=255)                           tttttGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt    >  1:8286984/1‑49 (MQ=38)                           tttttGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt    >  1:4539541/1‑49 (MQ=38)                           tttttGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTaa  >  1:3353871/1‑50 (MQ=11)                            ttttGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACt    >  1:5337786/1‑49 (MQ=11)                            ttttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa  >  1:1058691/1‑50 (MQ=255)                             tttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa  <  1:7968542/49‑1 (MQ=12)                             tttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa  <  1:4473974/49‑1 (MQ=12)                                                |                               GCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAA  >  NC_000913/4606285‑4606361

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