Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	579,194	+CT	intergenic (+2229/‑8935)	unknown → / → unknown	LSU rRNA ## 23S rRNA, large subunit ribosomal RNA/tRNA‑Pro‑TGG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	579,194	1	.	C	100.0%	36.3 / NA	21	intergenic (+2229/‑8935)	unknown/unknown	LSU rRNA ## 23S rRNA, large subunit ribosomal RNA/tRNA‑Pro‑TGG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (7/14);  total (7/14)
*	1	579,194	2	.	T	100.0%	36.0 / NA	21	intergenic (+2229/‑8935)	unknown/unknown	LSU rRNA ## 23S rRNA, large subunit ribosomal RNA/tRNA‑Pro‑TGG
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (7/14);  total (7/14)

CTAATGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAG  >  1/578964‑579265                                                                                                                                                                                                                                        ||                                                                        cTAATGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAag                                                       <  1:490904/251‑1 (MQ=17)    aTGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAgg                                                    >  1:237593/1‑251 (MQ=17)           ttttttCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAAt                                             <  1:624309/251‑1 (MQ=17)             ttttCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATgg                                           <  1:9734/251‑1 (MQ=17)             ttttCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATgg                                           <  1:798772/251‑1 (MQ=17)                           tCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                             >  1:155840/1‑251 (MQ=14)                           tCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                             >  1:853646/1‑251 (MQ=14)                           tCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCtt                             >  1:616262/1‑251 (MQ=14)                            ccTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAagag                                                     <  1:183913/226‑1 (MQ=14)                            ccTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTta                            >  1:323201/1‑251 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:607018/251‑1 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:691095/251‑1 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:898492/251‑1 (MQ=14)                               tAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTtatat                         <  1:110887/251‑1 (MQ=14)                                    tATGCAAATTTTACCGGATTTTTTGATCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAcc                    <  1:548875/251‑1 (MQ=14)                                     aTGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAAccc                   <  1:481439/251‑1 (MQ=14)                                                     aTTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt    <  1:35198/250‑1 (MQ=11)                                                     aTTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt    <  1:170623/250‑1 (MQ=11)                                                      ttttttGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAAt    >  1:25008/1‑249 (MQ=11)                                                      ttttttGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg  <  1:729281/251‑1 (MQ=11)                                                      ttttttGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg  >  1:840211/1‑251 (MQ=11)                                                                                                                                                                                                                                        ||                                                                        CTAATGCTATTTTTTTCATTTTAATGTCCTTAAATTATGCAAATTTTACCGGATTTTTTGACCGCACTTTGCTCTGCACCCCAAAACCTGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAG  >  1/578964‑579265

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