Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,331	Δ1 bp	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	.	98.2%	219.5 / ‑2.1	54	noncoding (71/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (30/23);  minor base G (0/1);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

TTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGC  >  NC_000913/4606198‑4606460                                                                                                                                       |                                                                                                                                  tttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGaa                                                                                                                             >  2:1138968/1‑140 (MQ=255) tttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGaa                                                                                                                             >  1:1144316/1‑140 (MQ=255) tttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGaa                                                                                                                             >  1:814760/1‑140 (MQ=255)   tGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAAcc                                                                                                                           >  2:1332360/1‑140 (MQ=255)     gTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAAcccc                                                                                                                         >  2:57601/1‑140 (MQ=255)       gcgagggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCa                                                                                                                       <  1:171572/135‑1 (MQ=255)         gagggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACg                                                                                                                     <  2:1245318/137‑1 (MQ=255)          agggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGt                                                                                                                    >  1:1145473/3‑140 (MQ=255)           gggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTc                                                                                                                   <  1:1138968/139‑1 (MQ=255)                 gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAg                                                                                                             <  1:394719/140‑1 (MQ=255)                  ggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAgg                                                                                                            <  1:590504/140‑1 (MQ=255)                     cTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGaca                                                                                                         <  1:376588/140‑1 (MQ=255)                         aaCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                     >  2:424285/1‑140 (MQ=255)                          aCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGCTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                                     <  2:156433/140‑1 (MQ=255)                            ccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTCCCAAGGCCGCCACCTTCGCCCAGTACGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAacac                                                                                                  <  2:1144316/140‑1 (MQ=37)                                 cGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGaa                                                                                             >  1:594856/1‑140 (MQ=255)                                        aaGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                                                                      >  2:561210/1‑140 (MQ=255)                                          gAACACTAACACCTGAAGCTAGCTCGTCTACCAATTCCGACACCTACGACCAATGCGGGCAATAACAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                                    <  2:289500/140‑1 (MQ=25)                                            acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                                                                  >  2:481533/1‑140 (MQ=255)                                            acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGCGGTTTTTGGAGCGAGGGGCGGG‑CTTGAACCCCCACGTCCGAAAGAACACTAACACCTGCAGCTAGCGCGTc                                                                                  >  2:588447/1‑140 (MQ=25)                                                tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTCCGCCCAGTGCGAGCATTATCAACGCGGTTTTTGGTGCGAGGGCGGGG‑CTTGAACCCCCACATCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                              >  1:297905/1‑140 (MQ=25)                                                 aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTACGAGCACTATCAACGTGGTTTTTGGTCCGAGGGGCAGG‑CTTGAACCCCCACGCCCGTAAGGACACTAACACCTGAGGCTAGCGCGTCTACCa                                                                             >  2:401137/1‑140 (MQ=25)                                                  acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                            >  2:623799/1‑140 (MQ=255)                                                  acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGAGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                            >  1:288679/1‑140 (MQ=255)                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                      <  1:57601/140‑1 (MQ=255)                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                      >  2:701189/1‑140 (MQ=255)                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                      <  2:1145473/140‑1 (MQ=255)                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                      >  1:19513/1‑140 (MQ=255)                                                          aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                                    <  1:1150068/140‑1 (MQ=255)                                                             tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                                                 <  2:497664/140‑1 (MQ=255)                                                              aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATGAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGAGACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                                >  2:970354/1‑140 (MQ=255)                                                                cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                              >  2:995251/1‑140 (MQ=255)                                                                  cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                            <  2:288679/140‑1 (MQ=255)                                                                   gTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCAt                                                           >  1:431711/1‑140 (MQ=255)                                                                           aTTCCGCCACCTTCGCCCAGTGCGTGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                                   <  1:501877/140‑1 (MQ=255)                                                                              ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTc                                                >  1:1145653/1‑140 (MQ=255)                                                                              ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGCGGGG‑CTTGAACCTCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTc                                                <  2:19513/140‑1 (MQ=255)                                                                               cGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCt                                               >  1:848797/1‑140 (MQ=255)                                                                                 ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTa                                             >  1:538969/1‑140 (MQ=255)                                                                                  cACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaa                                            >  1:405353/1‑140 (MQ=255)                                                                                         gCCCAGTGCGGGCAATATCAATGTGGTTTTTGGCGCGA‑GGGGGGGACTTGAACCCCAAAGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTCGCATACCATCAATTCTTAAAAAGAAt                                     <  2:297905/140‑1 (MQ=21)                                                                                                   aGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                           <  1:603097/140‑1 (MQ=255)                                                                                                        tatCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgc                      >  1:35882/1‑140 (MQ=255)                                                                                                         atCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg                          <  2:457956/135‑1 (MQ=255)                                                                                                         atCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg                          >  1:457956/1‑135 (MQ=255)                                                                                                         atCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                           <  2:928839/134‑1 (MQ=255)                                                                                                         atCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                           >  1:928839/1‑134 (MQ=255)                                                                                                         atCAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                     <  2:1145653/140‑1 (MQ=255)                                                                                                           cAACGTGGTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCa                   <  1:588447/140‑1 (MQ=255)                                                                                                               gTGGTTTTTGGTGCGAGTGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCt               >  1:1239401/1‑140 (MQ=255)                                                                                                                tGGTTTTGGGTGGGAGGGGGGGG‑CTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa              <  2:838848/140‑1 (MQ=21)                                                                                                                 ggTTTGTGGGGAGAGGGGGGGG‑CTTGAACACCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg             <  2:1269557/140‑1 (MQ=21)                                                                                                                  gTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                     <  1:1019790/131‑1 (MQ=255)                                                                                                                  gTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                     >  2:1019790/1‑131 (MQ=255)                                                                                                                  gTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt            >  1:1293373/1‑140 (MQ=255)                                                                                                                  gTTTTTGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt            >  1:561322/1‑140 (MQ=255)                                                                                                                   tttttGGTGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTg           <  1:701189/140‑1 (MQ=255)                                                                                                                          tGCGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCa    >  2:104752/1‑140 (MQ=255)                                                                                                                            cGAGGGGGGGG‑CTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGc  <  1:801404/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                  TTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTTTCAGC  >  NC_000913/4606198‑4606460

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