Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,215	(G)8→9	noncoding (72/87 nt)	leuP ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,207	1	.	G	95.1%	138.9 / 0.5	41	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (19/20);  total (20/21)
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.98e-01

CGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGG  >  NC_000913/4606084‑4606327                                                                                                                             |                                                                                                                          cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGa                                                                                                            >  1:442772/1‑139 (MQ=255) cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGa                                                                                                            >  1:212573/1‑139 (MQ=255) cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGa                                                                                                            <  2:381908/139‑1 (MQ=255) cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggggggggac                                                                                                               >  1:686997/1‑132 (MQ=255) cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGAAACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggggggggac                                                                                                               <  2:686997/135‑4 (MQ=255)               tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                              <  2:341489/139‑1 (MQ=255)                 caCCCCCACGTGCGTAAAGACACTAACATCTGAAGCTAGCGCGTCTACCAATTCCGGCGACTTCAAACAGTCAACTAACTTTTTTTGATATCGCCTAGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                            <  1:216764/138‑1 (MQ=25)                 aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                            <  2:212573/139‑1 (MQ=255)                 aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGAGAGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTaa                                                                                            <  1:446646/139‑1 (MQ=255)                      ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                           >  1:235037/1‑134 (MQ=255)                      ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                           <  2:235037/134‑1 (MQ=255)                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGCCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                       >  2:16290/1‑139 (MQ=255)                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                   >  1:645786/1‑127 (MQ=255)                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                   <  2:645786/127‑1 (MQ=255)                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                       >  2:404593/1‑139 (MQ=255)                                                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCa                                                           <  2:466731/139‑1 (MQ=255)                                                      aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTc                                                                                                 >  1:562989/1‑97 (MQ=255)                                                      aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTc                                                                                                 <  2:562989/97‑1 (MQ=255)                                                      aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑g                                                                                                                         >  2:672298/1‑72 (MQ=255)                                                      aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑g                                                                                                                         <  1:672298/72‑1 (MQ=255)                                                        cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGCCTGGACCCCCCACTCCCGTAACACTCCTAACGCCTGACAATGGCGAGTCTACCAATTCCg                                                     >  1:130745/1‑139 (MQ=255)                                                         gcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                    <  2:442772/139‑1 (MQ=255)                                                           gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                  >  1:351479/1‑139 (MQ=255)                                                              tACCAATTCCGCCACCTTCGCAAAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                <  2:766219/139‑1 (MQ=255)                                                                  aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACg                                                                                                   >  2:304229/1‑83 (MQ=255)                                                                  aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACg                                                                                                   <  1:304229/83‑1 (MQ=255)                                                                   aTTCCGCCAACTTCGCACAGTCATCTTACTTTTATTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                          <  1:426760/139‑1 (MQ=255)                                                                    ttCCGCCAACGTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                         <  1:774971/139‑1 (MQ=255)                                                                      ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                       >  1:648192/1‑139 (MQ=255)                                                                       cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                      <  2:747280/139‑1 (MQ=255)                                                                        gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                     <  2:765102/139‑1 (MQ=255)                                                                        gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                     <  2:648192/139‑1 (MQ=255)                                                                                cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa                              >  1:803286/1‑139 (MQ=255)                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                      >  2:596637/1‑139 (MQ=255)                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                      >  2:750763/1‑139 (MQ=255)                                                                                         aTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTg                    >  1:363311/1‑139 (MQ=255)                                                                                          tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTTCGTAAGAATACCTACACCTGAAGTTAACGCGTTTACCAATTTCGCCAACTTCGCCACGTGCGACCAATAACCACCTgg                   >  1:758328/1‑139 (MQ=18)                                                                                          tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                      >  1:746177/1‑120 (MQ=255)                                                                                          tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                      <  2:746177/120‑1 (MQ=255)                                                                                          tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                            <  1:755778/130‑1 (MQ=255)                                                                                          tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                            >  2:755778/1‑130 (MQ=255)                                                                                                 ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                   >  2:379492/1‑69 (MQ=21)                                                                                                 ttttttGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                   <  1:379492/69‑1 (MQ=21)                                                                                                         tatCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTAAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                   >  2:604394/1‑108 (MQ=255)                                                                                                         tatCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTAAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                   <  1:604394/108‑1 (MQ=255)                                                                                                           tCGGCTCGTTTTGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAggggg  <  2:130745/139‑1 (MQ=255)                                                                                                              ccTCGTTTGGTGCGAG‑GGGGGGGGACTTGAAccccc                                                                                                      >  1:660847/1‑36 (MQ=17)                                                                                                              ccTCGTTTGGTGCGAG‑GGGGGGGGACTTGAAccccc                                                                                                      <  2:660847/36‑1 (MQ=17)                                                                                                                             |                                                                                                                          CGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGG  >  NC_000913/4606084‑4606327

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