Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,518,092	G→T	noncoding (25/76 nt)	alaX ←	tRNA‑Ala

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,518,092	0	G	T	93.3%	65.9 / NA	30	noncoding (25/76 nt)	alaX	tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base T (12/16);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCA  >  NC_000913/2517851‑2518332                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  gATCGGGGTTTTAAGGAACACAAATACTAAGTTTATTCGCATCGAGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGCTCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTATccc                                                                                                                                                                                                                                            >  2:36458/1‑249 (MQ=255)               ggCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGccc                                                                                                                                                                                                                           <  1:160490/251‑1 (MQ=255)                     caAATAATAAGTTTATTCGCCTCGCGAAGCAAACAGGGGGCCAGCCTGGTGCCGATTTTTCATCCAATGAGGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCtata                                                                                                                                                                                                                               <  2:370815/241‑1 (MQ=255)                     caAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGAt                                                                                                                                                                                                                     <  1:2561/251‑1 (MQ=255)                     caAATAATAAGATTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCCCCAGGCGGCACTAACCCGGACGTAGCCTGTATTACAAACGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCtata                                                                                                                                                                                                                               >  1:370815/1‑241 (MQ=255)                      aaaTAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATg                                                                                                                                                                                                                    >  1:99599/1‑251 (MQ=255)                                 ttATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTa                                                                                                                                                                                                          <  1:107784/250‑1 (MQ=255)                                 ttATTCGCAGCGCGAAGCAAAAAAGGCGCCAGCCTCTTGACGCTTTTACATCCACTGATGCTTTCAACAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAACGCCGAGAAAAAAAACTCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGACGCTAAGCGGGATCGAACCGCTGACAACTTGCATGCCATGCAAGCTATCTCCGCG‑CTGAGCTATAATCCAACGATGCGTTTCCGTAc                                                                                                                                                                                                         >  2:262937/1‑251 (MQ=255)                                       ccATCGCGAATCAAACAGGGGGCCAGCCTGTTGGCGATTGTTCATCCAATGATGCTTTCGCCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATGTGTTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCTTGCAAGCTCTCTCCCAG‑CGGAGCTATAGCCCCACGAAGCTTTTACGTACCAAGtt                                                                                                                                                                                                   <  2:23050/250‑1 (MQ=38)                                           cgcgAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCt                                                                                                                                                                                               <  1:393649/251‑1 (MQ=255)                                            gcgAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTg                                                                                                                                                                                              <  2:142706/251‑1 (MQ=255)                                                   aaaCAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGGACCAAGTTTGCTGGGTGCa                                                                                                                                                                                        >  2:164466/1‑250 (MQ=255)                                                      cAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACAAGGCGGCACTAACCCGGACGAAGACTGTATTACAAAAGAGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGCGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAt                                                                                                                                                                                    >  1:163964/1‑251 (MQ=255)                                                                                                         ggcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGATGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGAGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCa                                                                                                                                   >  1:178630/1‑249 (MQ=255)                                                                                                         ggcACTAACCAGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGATGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGAAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCa                                                                                                                                   >  2:378773/1‑249 (MQ=255)                                                                                                             cTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAATCTCTCTCCCAG‑CTg                                                                                                                                                                                                                                      <  2:377294/146‑1 (MQ=255)                                                                                                             cTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTg                                                                                                                                                                                                                                      >  1:377294/1‑146 (MQ=255)                                                                                                                  ccGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTccc                                                                                                                        <  1:189560/251‑1 (MQ=255)                                                                                                                     gACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGc                                                                                                                     <  1:53923/251‑1 (MQ=255)                                                                                                                                        aagcatcAAAAATCAGTGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTTGTGTAGCTAAGCGTGATCGAACCTCTTCCCTCTTGCATGCCATGCTAGCTCTCTCCCAG‑CTCAGCTATAACCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTTGAACTAAGCGGGATCGAACCGCTTACCACTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCAc                                                                                                  <  1:178473/245‑1 (MQ=25)                                                                                                                                                 aaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAGGTTTGCTGGGTGCAAAATTTGGTGGAGCTGAGAGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAACGCtt                                                                                         >  2:44462/1‑251 (MQ=255)                                                                                                                                                            aaaaaaACCCAGCTTTGCAGCAAGGGATGGAAATTTGGTGGAGATAAGCCGGATCGAACGGCTGACCTCTTGCCTGCCATGCCAGCTCTATCCCAG‑CTGAGCTATAGCCCCCCGATGCGTTTACGTACCCAGTTTGCAGCGCGCAAACTTTGGTGCAGCTACGAGGGATCGAACCGCTGACCTCTTGCATGCCTTGCAAGCGCACTACAAGCTGAACTATAACACCGTCACGTAAAGATTGTCGAGTTgcc                                                                              >  2:359511/1‑249 (MQ=255)                                                                                                                                                                                     gTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGccc                                                     <  1:130750/251‑1 (MQ=255)                                                                                                                                                                                                 gtggAGCTAAGGGGGATCGAACCGCTGACCTCGTGCATGCCAGGCAAGCGCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAAc                                         <  2:172220/251‑1 (MQ=255)                                                                                                                                                                                                              ggATCGAACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                            <  1:164466/251‑1 (MQ=255)                                                                                                                                                                                                                agcGAACCGCTGACCTCATGCCTGCCAAGCCAGCTCTCGTCCAG‑CTGGTCCATAGACCCACGATTCGTTTACGGCCCAAGTTTGCTGGGTGCCAAATTTGGTTGAGCGAAGCGGGATCGAACCGCTGACCTCTTGCATGCCCTGCAAGCGCTCTCCCAGCTGATCTATAGCCCCGGCACGTAAAGCTTGTCGAGTTGACGGTCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGTAAAGTGCAACgt                          <  2:158294/249‑1 (MQ=21)                                                                                                                                                                                                                aTCGCACCGCTGACCTCTTGCATGCCATGCAAGCTCTCTCCCAG‑CGGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACgt                          <  2:6468/251‑1 (MQ=255)                                                                                                                                                                                                                          tGACCTCTTGCATGCCCTGCAAGCTCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGCAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCCTCACGTAACACTTGTCAAATTGACGGGCGGCATCATATGAATTCCGCTCGCATGTGTCACCGGCACATTGCAACgtgtcttttca                 >  2:330274/1‑243 (MQ=37)                                                                                                                                                                                                                                        cctttcAAGCTCTCTACCAG‑TTGAGCCATAGCCCCACGATGTGTTTACGTACTAAGTTTGCTGGTTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTTCCAGCTGTGCTATTGCCCCTTCACGTATAGCGGGTCGAGTTGACGGTCGGCATCTTATGAATTCCGCCCGAATGTGCCAACGGCAAATTGCAACTTGTAGTGTCAATCGCTGAATAATCa  <  2:391944/246‑1 (MQ=21)                                                                                                                                                                                                                                        ccagtcAAGCGCTCTCCCAGAGTG‑GCTATAGCCCCAGTAAGCGTTTACGTACCAAGTTTGCTGGGTGGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCAGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCa  <  2:15966/246‑1 (MQ=37)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  GATCGGGGTTTTCAGGCACACAAATAATAAGTTTATTCGCATCGCGAAGCAAACAGGGGGCCAGCCTGTTGCCGATTTTTCATCCAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAG‑CTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCA  >  NC_000913/2517851‑2518332

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