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.3%	170.6 / 50.6	85	noncoding (71/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (8/10);  new base G (27/38);  total (36/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

CGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAA  >  NC_000913/4606284‑4606361                                                 |                               cGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGa                                <  1:1563006/48‑1 (MQ=255)  gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCt                              <  1:2876793/49‑1 (MQ=255)  gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACtt                             >  1:4280827/1‑50 (MQ=255)  gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACtt                             >  1:5676747/1‑50 (MQ=255)  gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACtt                             >  2:4111372/1‑50 (MQ=255)   cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCtt                             >  2:4332167/1‑49 (MQ=255)   cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            <  1:4667429/50‑1 (MQ=255)    ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            >  1:491868/1‑49 (MQ=255)    ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            >  2:2307600/1‑49 (MQ=255)    ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTg                            >  2:665613/1‑49 (MQ=255)    ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGa                           <  1:2323576/50‑1 (MQ=255)    ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGa                           <  2:5085288/50‑1 (MQ=255)     cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGaa                          >  2:5098419/1‑50 (MQ=255)     cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGaa                          >  1:294964/1‑50 (MQ=255)     cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGaa                          >  2:597531/1‑50 (MQ=255)     cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGaa                          >  2:954396/1‑50 (MQ=255)     cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGaa                          <  1:4862999/50‑1 (MQ=255)     cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGaa                          >  1:3149981/1‑50 (MQ=255)      aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAc                         >  1:3440499/1‑50 (MQ=255)      aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAc                         >  1:808617/1‑50 (MQ=255)      aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAc                         <  1:4017458/50‑1 (MQ=255)      aGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAc                         <  2:560014/50‑1 (MQ=255)       gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcc                        <  2:250125/50‑1 (MQ=255)       gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcc                        <  1:4474820/50‑1 (MQ=255)       gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcc                        <  2:1414937/50‑1 (MQ=255)       gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcc                        <  1:2353698/50‑1 (MQ=255)       gTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGG‑CTTGAAccc                       >  1:71826/1‑50 (MQ=37)        tGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccc                       <  2:3150817/50‑1 (MQ=255)        tGCGAGCAATATAAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccc                       >  2:3619867/1‑50 (MQ=25)         gCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcccc                      >  1:416610/1‑50 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccccc                     <  2:2736777/50‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccccc                     <  2:2835582/50‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccccc                     <  1:4872652/50‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccccc                     <  2:3295600/50‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAccccc                     <  1:3343886/50‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcccc                      <  2:2155511/49‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcccc                      >  2:96591/1‑49 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAAcccc                      <  1:4556693/49‑1 (MQ=255)          cGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAAccccc                     <  1:62286/50‑1 (MQ=255)           gAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCa                    >  2:1203320/1‑50 (MQ=255)           gAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCa                    <  2:136612/50‑1 (MQ=255)            aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  1:150744/50‑1 (MQ=255)            aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   >  2:4161789/1‑50 (MQ=255)            aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCAc                   <  2:4114288/50‑1 (MQ=255)            aGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCAc                   >  1:4724818/1‑50 (MQ=255)             gCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACg                  <  1:3421196/50‑1 (MQ=255)             gCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGACTTGAACCCCCACGt                 <  2:1282208/50‑1 (MQ=25)              cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                 >  2:4712362/1‑50 (MQ=255)              cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGt                 >  1:4977880/1‑50 (MQ=255)               aaTATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGt                 >  2:4586469/1‑50 (MQ=21)               aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                >  1:65620/1‑50 (MQ=255)               aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                >  2:5124319/1‑50 (MQ=255)               aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                <  2:5547818/50‑1 (MQ=255)               aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                <  1:2321232/50‑1 (MQ=255)               aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTc                >  2:2941230/1‑50 (MQ=255)               aaTATCAACGTGGTTTTTGGTGCGA‑GGGGGGGG‑CTTGAACCCCCACGTcc               <  1:959272/50‑1 (MQ=21)                atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTcc               <  2:4740679/50‑1 (MQ=255)                atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTcc               <  2:3920942/50‑1 (MQ=255)                atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTcc               <  2:2356235/50‑1 (MQ=255)                atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTcc               <  1:1763567/50‑1 (MQ=255)                atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTcc               <  2:3834079/50‑1 (MQ=255)                   tCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGt             <  1:2432655/49‑1 (MQ=255)                   tCAACGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTa            <  1:1739552/50‑1 (MQ=255)                   tCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGt             <  1:5440364/49‑1 (MQ=255)                     aaCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  1:3358954/1‑50 (MQ=255)                     aaCGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg          <  1:4195614/50‑1 (MQ=255)                     aaCGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg          <  2:2642771/50‑1 (MQ=255)                     aaCGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg          >  2:1709712/1‑50 (MQ=255)                      tcGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg         <  2:3180667/49‑1 (MQ=255)                      aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAgg         <  2:5287715/50‑1 (MQ=255)                      aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAgg         <  1:5491285/50‑1 (MQ=255)                      aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  2:4308542/1‑49 (MQ=255)                      aCGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAg          >  2:5361815/1‑49 (MQ=255)                       cGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGa        <  2:3548648/50‑1 (MQ=255)                       cGTGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGa        <  2:5521332/50‑1 (MQ=255)                         tGGTTTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGac       <  1:4767667/50‑1 (MQ=14)                         tGGTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGaca      >  1:2363607/1‑50 (MQ=255)                          ggTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac     >  2:310632/1‑50 (MQ=255)                          ggTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac     >  2:1330703/1‑50 (MQ=255)                           gTTTTTGGTGCGAGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGaca      <  2:1048764/49‑1 (MQ=14)                           gTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac     >  1:3314660/1‑49 (MQ=255)                           gTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGacac     <  1:3894471/49‑1 (MQ=255)                           gTTTTTGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACt    <  2:516552/50‑1 (MQ=255)                            tttttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACCCt    >  2:4942440/1‑49 (MQ=11)                              tttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa  <  1:1927697/49‑1 (MQ=12)                              tttGGTGCGA‑GGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa  <  2:2943099/49‑1 (MQ=12)                                                 |                               CGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAA  >  NC_000913/4606284‑4606361

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