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	93.3%	89.4 / 0.4	30	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (10/18);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGATAATTGGTGC‑GA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGA  >  NC_000913/4606069‑4606336                                                                                                                                              |                                                                                                                                    tAAAGATAATTGGTGC‑GAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                      <  1:662133/140‑1 (MQ=255)       tAATTGGTGCGGG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggggg                                                                                                                              <  1:182037/140‑1 (MQ=255)       tAATTGGTGC‑GAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggggg                                                                                                                              <  1:439978/140‑1 (MQ=255)         aTTGGTGC‑GAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggggggg                                                                                                                            <  1:321255/140‑1 (MQ=255)         aTTGGTGC‑GAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggggggg                                                                                                                            >  1:689547/1‑140 (MQ=255)           tGGTGC‑GAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑ggggggggga                                                                                                                          >  1:281386/1‑138 (MQ=255)                 cgagggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACTAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGa                                                                                                                      <  1:767350/136‑1 (MQ=255)                 cgagggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                      >  2:702450/5‑125 (MQ=255)                 cgagggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                      <  1:702450/121‑1 (MQ=255)                   agggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAAc                                                                                                                    >  2:586132/3‑140 (MQ=255)                   agggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGTGGGGACTTGAAcc                                                                                                                   >  1:350294/3‑140 (MQ=255)                      gggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGCGAG‑GGGGGGGGACTTGAAcccc                                                                                                                 >  2:612108/1‑140 (MQ=255)                         ggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCa                                                                                                               >  1:530151/1‑140 (MQ=255)                                tGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                       >  2:455471/1‑140 (MQ=255)                                  aaCCCCCCCGTGCGTAAAGTCATTAACACCTGAAGATAGCGCGTCTTCCAAATCCGGCACCTAAACACAGTGAACTTACTATTATTGAAATCGACTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                     <  2:30286/140‑1 (MQ=19)                                      cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                 >  1:82353/1‑140 (MQ=255)                                       ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                        <  2:801117/132‑1 (MQ=255)                                       ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                        >  1:801117/1‑132 (MQ=255)                                       ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGTGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                               <  2:350294/140‑1 (MQ=255)                                        cccGTCCGTAAGGACACTAACACCTGAAGCTAGCACGTATTCCAATTCCGGCACCATCGCACAGGCATCTTACATTTATTGATATAGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                               <  2:134292/138‑1 (MQ=25)                                         aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                               <  1:76818/140‑1 (MQ=255)                                          cGTCCGTAAGGACACTACCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGGCATCTAACATTTATTGAGATCGCTTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                             <  2:737628/140‑1 (MQ=25)                                          cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                             <  2:281386/140‑1 (MQ=255)                                                      cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                 >  2:744475/1‑140 (MQ=255)                                                      cacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTa                                                                                 <  1:586132/140‑1 (MQ=255)                                                                         cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              >  2:445307/1‑140 (MQ=255)                                                                         cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              >  2:596/1‑140 (MQ=255)                                                                          gcgTCTACCAAATCCGGCACCTTCGAACAGTGATCTTACTTTTTTTGAAATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                              <  2:64095/140‑1 (MQ=21)                                                                                   aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                    <  1:445307/140‑1 (MQ=255)                                                                                   aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                    <  1:596/140‑1 (MQ=255)                                                                                   aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                    <  1:744475/140‑1 (MQ=255)                                                                                    aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGTCACCTTCGcc                                                    <  2:530151/140‑1 (MQ=21)                                                                                          ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAT‑GGGGGGGGACTTGAATTCCCACGACCGTAAGCACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAgaacc                                             >  2:320491/1‑136 (MQ=21)                                                                                          ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                          <  2:337355/127‑1 (MQ=255)                                                                                          ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                          >  1:337355/1‑127 (MQ=255)                                                                                             ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                             <  1:510086/137‑1 (MQ=255)                                                                                             ccTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                             >  2:510086/1‑137 (MQ=255)                                                                                                      aGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                 <  1:455471/140‑1 (MQ=255)                                                                                                                                 tCGTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCTGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑gggggggggac     <  1:612108/140‑4 (MQ=255)                                                                                                                                   gTTTGGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTg   <  2:177708/140‑1 (MQ=255)                                                                                                                                     ttGGTGCGAG‑GGGGGGGGGCCAGAACCACCAGGTCAGAAAGAACACTAGCACCTGAAGCTAGCGCGTATTCCAAGTCCGCCACCCTCGCCAAGTGCGAGAAATAAAATAGTGGTTTTTGGTGCGAGGGGGGGGGACTTGa  <  2:584720/140‑1 (MQ=255)                                                                                                                                        gTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑ggggggggg       >  2:535849/1‑130 (MQ=255)                                                                                                                                        gTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑ggggggggg       <  1:535849/131‑2 (MQ=255)                                                                                                                                              |                                                                                                                                    TAAAGATAATTGGTGC‑GA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGA  >  NC_000913/4606069‑4606336

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: ‑