Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	275,065	G→A	T18I (ACC→ATC)	insH1 ←	IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	275,065	0	G	A	94.7%	42.4 / ‑6.7	19	T18I (ACC→ATC)	insH1	IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/10);  total (8/11)
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

AAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTA  >  NC_000913_3_pae_tpiA/275012‑275209                                                      |                                                                                                                                                 gtatgtcnTGTCTGGTCTGATGGTGGTTAGCGGTAGATTCACTGTCGGAGACGATGAGTTGGCGACTCATGATGAACCCTGTTCTATGGTTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTAt                                                   <  2:829676/144‑1 (MQ=11) aaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTGCCGGAGGCGTTAt                                                   >  1:776848/1‑149 (MQ=255) aaaTCCCGTTGCTGGTTGTACGTCGCTTATTGGTGAATTCACGGTCGGGGAAGAGAAGTTGCTGACTCATGATGAACCCTGTTCTATGGTTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTAt                                                   <  2:828295/149‑1 (MQ=11)        ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGATCTgg                                            >  1:153573/1‑149 (MQ=255)        ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTgg                                            >  2:128603/1‑149 (MQ=255)         tttCTGGTCGGACGGCGCTTAGTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggc                                           <  2:186020/149‑1 (MQ=255)          ttCCGGTCTGACGGCGCTTACTGGTGAATTCACTGTCGGGGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGTTGACAAACATGATCTCATATTAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggcg                                          <  2:499447/149‑1 (MQ=14)           tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggcgg                                         <  2:184871/149‑1 (MQ=255)           tCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggcgg                                         <  2:81569/149‑1 (MQ=255)             tGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCACGAACTGATTCGCACCTTCCGGAGGCGGTATGAGCTGGCGGcc                                       >  1:874996/1‑149 (MQ=14)                   gaaggCGCATACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                                 <  2:531939/146‑1 (MQ=255)                                gctgAGTTCACTGTCGGCGAAGATAAGTTGATGACTGATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGattatt                    <  2:528872/149‑1 (MQ=255)                                 ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATAAGCTGGCGGCCCATTTTGTATCTCGttatta                   >  1:353474/1‑149 (MQ=17)                                   gAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAAt                 <  2:211134/149‑1 (MQ=255)                                       tCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATcccc             >  1:62950/1‑149 (MQ=255)                                        cACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCa            <  2:589821/149‑1 (MQ=255)                                                gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                                 >  1:742132/1‑120 (MQ=255)                                                gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                                 <  2:742132/120‑1 (MQ=255)                                                  gAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTACTCCCCACCCGCTATTa  >  1:263376/1‑149 (MQ=255)                                                      |                                                                                                                                                 AAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTA  >  NC_000913_3_pae_tpiA/275012‑275209

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