Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,724,478	(G)6→7	coding (178/522 nt)	sodC ←	superoxide dismutase (Cu‑Zn)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,724,472	1	.	G	100.0%	93.8 / NA	27	G62R (GGT→CGT)	sodC	superoxide dismutase (Cu‑Zn)
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (14/13);  total (14/13)

CACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑GGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGT  >  NC_000913/1724337‑1724605                                                                                                                                         |                                                                                                                                      cACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggg                                                                                                                                   <  1:289324/140‑1 (MQ=255)   ccTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑gggggg                                                                                                                                 <  2:703011/140‑1 (MQ=255)              ttCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTc                                                                                                                       >  1:643065/1‑140 (MQ=255)                cATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAg                                                                                                                     >  1:95856/1‑140 (MQ=255)                            tATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGaaaa                                                                                                         >  2:365872/1‑140 (MQ=255)                                    ggATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGAc                                                                                                 <  1:969230/140‑1 (MQ=255)                                     gATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGAcc                                                                                                >  1:898415/1‑140 (MQ=255)                                           gATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATc                                                                                          <  2:95856/140‑1 (MQ=255)                                           gATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATc                                                                                          <  2:643065/140‑1 (MQ=255)                                              gCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                                                                                       <  1:1359/140‑1 (MQ=255)                                                ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCt                                                                                               <  1:305213/130‑1 (MQ=255)                                                ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCt                                                                                               >  2:305213/1‑130 (MQ=255)                                                ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGttt                                                                                     >  2:513276/1‑140 (MQ=255)                                                ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                                                                                       >  1:428106/1‑138 (MQ=255)                                                ccGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGt                                                                                       <  2:428106/138‑1 (MQ=255)                                                             ccGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGAc                                                                        <  2:898415/140‑1 (MQ=255)                                                                   cgcTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTAcc                                                                  <  1:365872/140‑1 (MQ=255)                                                                         cTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                            <  1:407245/140‑1 (MQ=255)                                                                         cTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                            <  1:513276/140‑1 (MQ=255)                                                                         cTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGa                                                            <  1:1010189/140‑1 (MQ=255)                                                                            tGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTg                                                         >  2:623124/1‑140 (MQ=255)                                                                                      ggtggCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCt                                               >  2:792673/1‑140 (MQ=255)                                                                                             ggctggcAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTAc                                                   <  1:856154/129‑1 (MQ=255)                                                                                             ggctggcAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTAc                                                   >  2:856154/1‑129 (MQ=255)                                                                                             ggctggcAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACg                                        >  1:148127/1‑140 (MQ=255)                                                                                                                    aaTATGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACtttt                 >  1:735947/1‑140 (MQ=255)                                                                                                                       atGGAAGCCATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCa              <  1:939493/140‑1 (MQ=255)                                                                                                                              ccATGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAg       >  1:425874/1‑140 (MQ=255)                                                                                                                                aTGTTCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCt     >  2:892784/1‑140 (MQ=255)                                                                                                                                   ttCACCGGGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTtgt  >  1:318893/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                      CACCTTCTGGCCCTTCATGTTTACCGGTATTTTGTGGATCAAGATGCCCGCCTGCGGATTCCGCGGCGCTGGCTTTGCCATCTTTGGTGGCTGGCTGGCAGCTTCCTTTGGCATGAATATGGAAGCCATGTTCACC‑GGGGGGTAATGCTTTCAGATCGGGCGAAAACTCCAGACCTTTATCGGTTTCAGTAATGGTGACGCTACCAATTGACTGCCCTACCCCTTGCGACGTGACGAGGTTCATCTCGACTTTTTCACTGGCAGCTTGT  >  NC_000913/1724337‑1724605

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