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	87.2%	106.4 / 10.4	39	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base G (12/22);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACG  >  NC_000913/4606068‑4606345                                                                                                                                              |                                                                                                                                                gTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                                  <  1:26061/140‑1 (MQ=255)        tAATTGGTGGGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑ggg‑g‑gg                                                                                                                                        <  1:556013/140‑1 (MQ=255)        tAATTGGTGCGA‑GGGCGGGGACCTGGACCCCCACGTCAGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑ggg‑g‑ggg                                                                                                                                       <  1:309817/140‑1 (MQ=255)               tGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTg                                                                                                                                <  2:25126/140‑1 (MQ=255)                gCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGa                                                                                                                               <  1:960662/140‑1 (MQ=255)                 cGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCGTACTTTTTGTGAGATCTCCTCGTTTGGTGCAAG‑GGG‑G‑GGG‑GACTTGaa                                                                                                                              >  2:270938/1‑140 (MQ=18)                         caggACTTGAACCCCCAAGTCCATAAGGACATTAAAACCAGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCAc                                                                                                                       <  2:307852/138‑1 (MQ=25)                          gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACg                                                                                                                      <  2:664271/140‑1 (MQ=255)                          gggACTTGAACCCCCACGTCCGCAAAGACCCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTGCTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACg                                                                                                                      <  1:270938/140‑1 (MQ=25)                            gACTTGAATCCCCACGTCCGTACGGACACTAGCGCCTGAGGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTc                                                                                                                    <  2:458998/140‑1 (MQ=25)                               ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGt                                                                                                                 <  1:549554/140‑1 (MQ=255)                               ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGATGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGt                                                                                                                 <  2:605575/140‑1 (MQ=255)                                 gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTaa                                                                                                               <  1:248407/140‑1 (MQ=255)                                  aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAg                                                                                                              <  1:895740/140‑1 (MQ=255)                                  aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAg                                                                                                              <  2:728742/140‑1 (MQ=255)                                      cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAac                                                                                                           <  1:810185/140‑1 (MQ=255)                                      cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAac                                                                                                           <  1:967006/140‑1 (MQ=255)                                       ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGG‑‑GGG‑G‑GGGAGACGTGAACCCCCACGTCCGTAAGAactg                                                                                                         >  1:307852/1‑138 (MQ=37)                                          cGTCCGTAAGGACACTACCACCAGAAGCTAGAGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                      <  2:373029/140‑1 (MQ=38)                                          cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCGGTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                      <  2:586188/140‑1 (MQ=255)                                          cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAATACTaa                                                                                                      <  2:920703/140‑1 (MQ=255)                                          cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                      <  1:648008/140‑1 (MQ=255)                                                tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCCCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAAGCCCAACGTCCGTAGGAACACTAACACCTg                                                                                                >  1:373029/1‑140 (MQ=37)                                                         tAACACCTGAAGGTAGCCCGTCAACCAATTCGGTGTCCTTCGCACAGTAATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                       <  2:445345/140‑1 (MQ=25)                                                           acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGcgcg                                                                                     >  1:644673/1‑140 (MQ=255)                                                           acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACCTGAACCCCCACGTCCGAAAGAACCCTAACACCCGTAGCTaacacg                                                                                     >  2:86330/1‑135 (MQ=37)                                                                 gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCGGTAAGAATACTGACACCTGAACCTAGCGCGTCTACCa                                                                              >  1:451320/1‑140 (MQ=255)                                                                   aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                             >  1:973289/1‑140 (MQ=255)                                                                       aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGTG‑GGG‑GACTTCAACCCCCACGTCCGTCAGTACACTAACACCTGAAGCTAGCGCGACTACCAATTcc                                                                         >  2:897240/1‑140 (MQ=25)                                                                        gcgcGTCTACCAATTCCGCCACATTCGCACAATAAACTTCCTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                       <  2:451320/140‑1 (MQ=21)                                                                         gacagcTACCAATTCCGCCACCTTCGCACCGTCATCTCTCTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                       <  1:302301/135‑1 (MQ=21)                                                                         cgcgTCTACCAATTCCGCCACCATCCCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                       <  1:367491/140‑1 (MQ=255)                                                                                     ttCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCa                                                           >  2:543871/1‑140 (MQ=255)                                                                                        cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGGGTCAACTAATTCCCCCACCTTCGCCCAGt                                                         >  1:383378/1‑140 (MQ=21)                                                                                        cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                                       >  2:657169/1‑140 (MQ=255)                                                                                         gCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACTTCCGTAAGAACACTAACACCTGATGCTGGCGGGTCTACCAGTTCCGCCACCTTCGcccctgtc                                                       >  2:302301/1‑135 (MQ=21)                                                                                          gcaccgacGCACAATGATCTTACTTTTTTTGATAGCGAAGCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                      <  1:76506/133‑1 (MQ=17)                                                                                            aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                    <  1:188924/140‑1 (MQ=255)                                                                                               ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                    >  2:29340/1‑137 (MQ=255)                                                                                               ttCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                                    <  1:29340/137‑1 (MQ=255)                                                                                                 cgtacaGTCATCTTACTTTTTCTGCCATCGCCACCTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACGCCCCCGTCCATGAGACCACTAACACCTGAAGCTAGTGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAata                                               <  2:973289/137‑1 (MQ=18)                                                                                                        tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGAGCGTCTACCAATACTGGCACCTTCGTCCACAACGAGCAATATCAACga                                        >  2:450471/1‑139 (MQ=18)                                                                                                                cTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                                >  2:64826/1‑140 (MQ=255)                                                                                                                    ttttGATATCGCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTGCGCCCAGGGCGAGCCATATCAACGTGGTTTTTGGTGc                            >  1:287510/1‑140 (MQ=21)                                                                                                                           atCGCCTCGTTTGGTGCGAGGGGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCACCGCCAACAATATCAACGAGGTTTTTGGAGCGA‑gggg                     >  2:160282/1‑140 (MQ=18)                                                                                                                              gCCTCGTTTGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑gggggggg                 <  1:543871/140‑1 (MQ=255)                                                                                                                                 tCGGTTGGTGCGA‑‑GGG‑GTGGG‑GGATTGAACCCCCCCGTCCGTAAGAACACTAACCCCTTAAGCCAGCGCGTCTACCAATTCCGCCACCCTCGCCCAGTGCGAGTAATATCAACGTGGTTTTTGGTgggg‑gggggggggga              <  2:287510/140‑4 (MQ=18)                                                                                                                                      tGGTGCGAG‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGaa          <  2:644673/140‑1 (MQ=255)                                                                                                                                            gA‑‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg  <  1:657169/140‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                GTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGG‑G‑GGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACG  >  NC_000913/4606068‑4606345

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