Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	428,939	T→G	intergenic (‑327/+465)	unknown ← / ← unknown	hypothetical protein/LSU rRNA ## 23S rRNA, large subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	428,939	0	T	G	100.0%	58.3 / NA	23	intergenic (‑327/+465)	unknown/unknown	hypothetical protein/LSU rRNA ## 23S rRNA, large subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (6/17);  total (6/17)

CGTTTATCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTTTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACAC  >  1/428702‑429046                                                                                                                                                                                                                                              |                                                                                                            cctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTacggttttgttt                                                                                                        <  1:551573‑M2/236‑5 (MQ=255)  ttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACaa                                                                                               <  1:398054‑M2/245‑1 (MQ=255)          gTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAg                                                                                        >  1:481458/1‑250 (MQ=255)          gTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGttttgttt                                                                                                        <  1:171621/234‑5 (MQ=11)           ttCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTta                                                                                                          <  1:555390/231‑1 (MQ=17)           ttCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCt                                                                                      <  1:683253/251‑1 (MQ=255)                    ttCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACaaaa                                                                             >  1:423367/1‑251 (MQ=255)                                     cttctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTca                                                             >  1:697884/1‑250 (MQ=255)                                        ctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:383538/251‑1 (MQ=255)                                        ctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:431610/251‑1 (MQ=255)                                           tatTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAgg                                                      <  1:750447/251‑1 (MQ=255)                                                        aGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAgg                                         >  1:605527/1‑251 (MQ=255)                                                          aCGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCTTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGc                                        <  1:447879/250‑1 (MQ=255)                                                          aCGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGc                                        <  1:707975/250‑1 (MQ=255)                                                           cGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:839331/232‑1 (MQ=255)                                                             aTGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:864807/230‑1 (MQ=255)                                                                               cATTTGGTATGTGATAGACTTGAATCGTCAGTGGCTTTACTTACCATTGACTTATCCGCCTGTATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGGCAACGTCTGGCTATAACAAAAAGCCCCGATGCCCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGGCCTACTCTcaca                  <  1:300667/251‑1 (MQ=14)                                                                                     gTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGa            >  1:320201/1‑251 (MQ=37)                                                                                         gtgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGc                                                                                       >  1:214452/1‑172 (MQ=255)                                                                                         gtgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGAcc                            <  1:255498/231‑1 (MQ=25)                                                                                         gtgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATgggg             <  1:101300/246‑1 (MQ=25)                                                                                          tgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGccc       <  1:739401/251‑1 (MQ=21)                                                                                          tgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCCCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGccc       <  1:677476/251‑1 (MQ=21)                                                                                                aCTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCacac  <  1:854493/250‑1 (MQ=17)                                                                                                aCTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCacac  <  1:856424/250‑1 (MQ=17)                                                                                                                                                                                                                                              |                                                                                                            CGTTTATCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTTTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACAC  >  1/428702‑429046

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