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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,086	1	.	G	94.7%	127.6 / 0.4	38	noncoding (80/87 nt)	leuV	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (14/22);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 5.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

CACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA  >  NC_000913/4605948‑4606207                                                                                                                                             |                                                                                                                          cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑g                                                                                                                          >  2:230617/1‑140 (MQ=255) cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑g                                                                                                                          <  1:572949/140‑1 (MQ=255)   caAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑ggg                                                                                                                        >  2:200170/1‑140 (MQ=255)     aaCAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑ggggg                                                                                                                      >  2:440574/1‑140 (MQ=255)      aCAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑gggggg                                                                                                                     <  2:540526/140‑1 (MQ=255)        aGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑gggggggg                                                                                                                   >  1:468725/1‑140 (MQ=255)         gTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑ggggggggg                                                                                                                  <  2:577754/140‑2 (MQ=255)          ttGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑ggggggggga                                                                                                                 >  1:784626/1‑138 (MQ=255)              ggCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGa                                                                                                             <  2:174019/140‑1 (MQ=255)                    aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAAcccc                                                                                                        <  2:468725/140‑1 (MQ=255)                        tttGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACg                                                                                                    <  2:315647/140‑1 (MQ=255)                          tGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAG‑AAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTcc                                                                                                 >  2:357316/1‑140 (MQ=255)                                cATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                            <  1:440574/140‑1 (MQ=255)                                        cTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAGATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                    <  2:717923/140‑1 (MQ=255)                                         tGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                                                                                     >  1:396622/1‑122 (MQ=255)                                         tGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCAc                                                                                                     <  2:396622/122‑1 (MQ=255)                                            cACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACAcc                                                                                >  2:333271/1‑140 (MQ=255)                                            cACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGAAACTAACAcc                                                                                >  2:526758/1‑140 (MQ=255)                                                 aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAg                                                                           >  2:678960/1‑140 (MQ=255)                                                     cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                       <  2:508982/140‑1 (MQ=255)                                                     cTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                       >  1:609348/1‑140 (MQ=255)                                                         catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                   >  1:141991/1‑140 (MQ=255)                                                         catcatGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                   <  1:200170/140‑1 (MQ=255)                                                                      atAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAG‑AAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     <  1:357316/140‑1 (MQ=255)                                                                       tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     <  1:678960/140‑1 (MQ=255)                                                                       tAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     <  2:384015/140‑1 (MQ=255)                                                                        aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                    >  1:524529/1‑140 (MQ=255)                                                                        aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGcc                                                    >  1:493919/1‑140 (MQ=255)                                                                        aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     >  1:429179/1‑139 (MQ=255)                                                                        aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     <  2:429179/139‑1 (MQ=255)                                                                         aTGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                   >  1:75105/1‑140 (MQ=255)                                                                          tGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                  <  2:542728/140‑1 (MQ=255)                                                                           gAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     >  1:258130/1‑136 (MQ=255)                                                                           gAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     <  2:258130/136‑1 (MQ=255)                                                                                 gAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                          <  2:6662/140‑1 (MQ=255)                                                                                   attaaggATGGAAAATCCAGGGAGAAAACAGGAAAGAAATTAACGATACTTGGGCCGG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                                         <  2:280779/135‑1 (MQ=14)                                                                                   aTCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcaca                                         <  1:118730/140‑1 (MQ=255)                                                                                   aTCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTGCGcaca                                         <  2:342821/140‑1 (MQ=255)                                                                                       gggATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCa                                     <  1:367280/140‑1 (MQ=255)                                                                                            gAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTa                                <  1:792950/140‑1 (MQ=255)                                                                                                  tCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtttct                          >  1:13867/1‑138 (MQ=255)                                                                                                    aaGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttttt                        >  1:591864/1‑140 (MQ=255)                                                                                                        aag‑aaACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGata                    <  2:75105/140‑1 (MQ=255)                                                                                                             aaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                <  2:609348/140‑1 (MQ=255)                                                                                                             aaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                <  2:524529/140‑1 (MQ=255)                                                                                                             aaCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGc                <  1:231234/140‑1 (MQ=255)                                                                                                                           tAAAGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa  <  2:229163/140‑1 (MQ=255)                                                                                                                           caaaGATAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa  <  1:21025/139‑1 (MQ=255)                                                                                                                                             |                                                                                                                          CACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAG‑AAACAAGAAAGGAAGTAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA  >  NC_000913/4605948‑4606207

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