Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,432,994	(C)6→5	coding (567/1125 nt)	smf ←	protein Smf

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,432,989	0	C	.	100.0%	80.9 / NA	20	coding (572/1125 nt)	smf	protein Smf
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (8/12);  total (8/12)

TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATG  >  NC_000913/3432850‑3433104                                                                                                                                            |                                                                                                                    tGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAgcgc                                                                                                                     >  1:86409/1‑140 (MQ=255)     gTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGccccc                                                                                                                 <  1:149477/140‑1 (MQ=255)                    cTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTg                                                                                                 <  2:86409/140‑1 (MQ=255)                                         accaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                            <  1:209171/140‑1 (MQ=255)                                         accaGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGAcgg                                                                            <  2:37707/140‑1 (MQ=255)                                             gTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCggg                                                                        >  1:154828/1‑140 (MQ=255)                                              tACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCaggg                                                                       >  2:393952/1‑140 (MQ=255)                                                cacCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGAt                                                                     <  1:132665/140‑1 (MQ=255)                                                  ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGa                                                                   <  2:348594/140‑1 (MQ=255)                                                  ccTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCAGGGATGa                                                                   <  1:393952/140‑1 (MQ=255)                                                          tAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGGTGAATGGTAtt                                                           <  1:389541/140‑1 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                        >  2:212246/1‑140 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAg                                                        >  1:172186/1‑140 (MQ=255)                                                             aCCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGGATTAAg                                                        >  2:291233/1‑140 (MQ=255)                                                                   gATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCAtt                                                  >  2:402226/1‑140 (MQ=255)                                                                         gCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCaa                                            <  2:172186/140‑1 (MQ=255)                                                                              ttcttcGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAg                                       >  2:320415/1‑140 (MQ=255)                                                                                         aaaTTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACg                            <  2:154828/140‑1 (MQ=255)                                                                                                    aaGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTg                 <  1:402226/140‑1 (MQ=255)                                                                                                            gAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                      <  2:293355/127‑1 (MQ=255)                                                                                                            gAACATCGAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCAttt                      >  1:293355/1‑127 (MQ=255)                                                                                                                   gAGGGGAAATTCCGAGACGAGAGCG‑CCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATg  <  1:212246/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                    TGCTGTCACCAGCGAACCACTACGCAAAGCCGCTTCCACCACCAGTACACCTTTACTTAGACCACTGATAATGCGATTTCTTCGTGGGAAATTGTAAGCAAGGGGTGGAACATCGAGGGGAAATTCCGAGACGAGAGCGCCCCCCTGTTCAAGCAGACTGGCAGCCAGTCGGGCATGACGGCGGGGATGAATGGTATTAAGTCCATTCCCCAATACAGCAATGCTGACGCCATTTACCTGTAAGGCTGCTTTATG  >  NC_000913/3432850‑3433104

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