Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,982,684	(C)8→9	noncoding (79/87 nt)	leuT →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,982,676	1	.	C	96.3%	90.7 / ‑2.6	27	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (11/15);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

TGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGA  >  NC_000913/3982526‑3982818                                                                                                                                                        |                                                                                                                                               tGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGAGAAGGTGGAGGAATTGGTAGACGCGCTAGCTTCAGGCGTTAGTGTCCTTACGGACGTGGGGGTTCAAGt                                                                                                                                                 >  2:661435/1‑151 (MQ=255)             ggATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑ccccccccctc                                                                                                                                     <  2:354533/150‑4 (MQ=255)                 tgtgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGcacc                                                                                                                                 >  1:116983/1‑151 (MQ=255)                 tgtgATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGcacc                                                                                                                                 >  1:213798/1‑151 (MQ=255)                        ccAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACtt                                                                                                                          <  1:674361/151‑1 (MQ=255)                          aGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGCCGTGGGGTTCCAAGTCCCCCCCCCTCGCACCACGACTTTa                                                                                                                        <  2:47328/151‑1 (MQ=255)                          aGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACttt                                                                                                                         >  2:328343/1‑150 (MQ=255)                          aGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACttt                                                                                                                         >  2:213743/1‑150 (MQ=255)                           gTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTa                                                                                                                        <  1:567114/150‑1 (MQ=255)                                cGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTCCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGaa                                                                                                                   <  2:578172/150‑1 (MQ=255)                                     gTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAAc                                                                                                             <  1:256965/151‑1 (MQ=255)                                        cGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTa                                                                                                           <  2:223356/150‑1 (MQ=255)                                             cccATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAtt                                                                                                     <  1:567037/151‑1 (MQ=255)                                             cccATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAt                                                                                                      <  1:357010/150‑1 (MQ=255)                                                   aGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                                                                               <  1:434561/151‑1 (MQ=255)                                                     ccACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                                                                               >  2:435316/1‑149 (MQ=255)                                                               attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                                                                    <  1:537422/150‑1 (MQ=255)                                                                  aGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCg                                                                                 >  2:446793/1‑150 (MQ=255)                                                                    aaGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGt                                                                              >  2:377450/1‑151 (MQ=255)                                                                      gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                                                                            <  1:231433/151‑1 (MQ=255)                                                                         gtgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCa                                                                        <  2:590536/151‑1 (MQ=255)                                                                               aTGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTg                                                                   <  2:450694/151‑1 (MQ=255)                                                                                            ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTgg                                                      >  2:80014/1‑151 (MQ=255)                                                                                                     aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGCGTAGCGCAGCTTGGGAGCGCAACTGGTTTGGGAcc                                             >  1:553201/1‑151 (MQ=255)                                                                                                     aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGAcc                                             >  1:206484/1‑151 (MQ=255)                                                                                                                  tCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAg                                 <  1:344215/150‑1 (MQ=255)                                                                                                                   cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAgg                                <  1:536190/150‑1 (MQ=255)                                                                                                                                  ttACGGACGTGGGGGTTCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCg                >  2:674655/1‑151 (MQ=255)                                                                                                                                                 ttCAAGTCCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGAACCAGTGGGGCGGAGGTTCGAATCCTCTCTCGCCGACCAAGTTTGa  <  1:250181/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               TGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGT‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGA  >  NC_000913/3982526‑3982818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑