Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,217,252	C→A	L151I (CTC→ATC)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,217,252	0	C	A	100.0%	104.7 / NA	35	L151I (CTC→ATC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (19/16);  total (19/16)

ATTACTGCATGAATCTGTTGCGGACTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATATTCACAAGCACATCGGTACGATTTCTATCAACACGACTCTCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGACAGTTTTTATCTTTATAAAATGAATATTAATTTTTTAT  >  1/1217055‑1217492                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 aTTACTGCATGAATCTGTTGCGGACTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTAt                                                                                                                                                                                             <  1:43700/251‑1 (MQ=255)     cTGCATGAATCTGTTGCGGACTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTAttt                                                                                                                                                                                           >  1:325748/1‑249 (MQ=255)            aaTCTGTTGCGGACTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATcgcg                                                                                                                                                                                                                                             <  1:556309/192‑1 (MQ=255)                         cTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATcg                                                                                                                                                                                                                                               >  1:34594/1‑177 (MQ=255)                         cTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTAc                                                                                                                                                                                                       <  1:259593/217‑1 (MQ=255)                         cTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTAtcgt                                                                                                                                                                      >  1:805659/1‑250 (MQ=255)                         cTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTAtcgt                                                                                                                                                                      >  1:770713/1‑250 (MQ=255)                                   gtgtAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTcc                                                                                                                                                            >  1:846402/1‑250 (MQ=255)                                              tcttTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACAttt                                                                                                                                                 >  1:457845/1‑250 (MQ=255)                                                                  ggTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTg                                                                                                                                                                                         <  1:465678/190‑1 (MQ=255)                                                                       cATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGtt                                                                                                                                                                            <  1:114795/198‑1 (MQ=255)                                                                                 gTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCaa                                                                                                             >  1:76464/1‑251 (MQ=255)                                                                                      aaaaGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGcc                                                                                                         <  1:109376/250‑1 (MQ=255)                                                                                          gTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTcc                                                                                                                                                            >  1:848162/1‑195 (MQ=255)                                                                                           ttAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTa                                                                                                    <  1:422394/250‑1 (MQ=255)                                                                                            tAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACg                                                                                                  <  1:821688/251‑1 (MQ=255)                                                                                                          cAATCTTCTATATTGGCATTACTGTTATCTTTACTCTTTTGTTGAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATGTATGTGATGTTGGCTGCCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGGCTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATAc                                                                                    <  1:346634/251‑1 (MQ=255)                                                                                                            aTCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTAtcg                                                                                                                                                                       >  1:49950/1‑166 (MQ=255)                                                                                                              cttctATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCg                                                                                                                                                                              >  1:584168/1‑157 (MQ=255)                                                                                                              cttctATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATAc                                                                                    <  1:208004/247‑1 (MQ=255)                                                                                                                   atatTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTAtcgtc                                                                                                                                                                     <  1:895919/161‑1 (MQ=255)                                                                                                                   atatTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTAc                                                                                                   >  1:867881/1‑227 (MQ=255)                                                                                                                        ggCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGt                                                                                                 >  1:221855/1‑224 (MQ=255)                                                                                                                                     tCTTTACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGt                                                         >  1:22747/1‑251 (MQ=255)                                                                                                                                        ttACTCTTTTGTTTAAAAATTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATc                                                                     <  1:573650/236‑1 (MQ=255)                                                                                                                                                       aaaaTTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGAc                                       <  1:146566/251‑1 (MQ=255)                                                                                                                                                          aTTATACTCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGacac                                                                            >  1:59425/1‑211 (MQ=255)                                                                                                                                                                ctCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGACAGTTTTTAt                              >  1:665274/2‑251 (MQ=255)                                                                                                                                                                 tCACAAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGACAGTTTTTATc                             <  1:598189/251‑1 (MQ=255)                                                                                                                                                                  cacaAGCACATCGGTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACtt                                                                          >  1:621115/1‑205 (MQ=255)                                                                                                                                                                              ggTACGATTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGACAGTTTTTATCTTTATAAAATGaa                >  1:279100/1‑251 (MQ=255)                                                                                                                                                                                    aTTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTg                                                  >  1:762582/1‑211 (MQ=255)                                                                                                                                                                                    aTTTCTATCAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGcta                                               <  1:304856/214‑1 (MQ=255)                                                                                                                                                                                            cAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAata                                                                                                                                                                                    <  1:662112/73‑1 (MQ=255)                                                                                                                                                                                            cAACACGACTATCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGACAGTTTTTATCTTTATAAAATGAATATTAATTTTTTAt  >  1:849709/1‑251 (MQ=255)                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 ATTACTGCATGAATCTGTTGCGGACTTATATTCTGTGTAATGTCTTCTTTATTCCGGTTATTCTTGGTTACATCCTGACAGTTGGAAAAGTTAATTATCTTTTCGCAATCTTCTATATTGGCATTACTTTTATCTTTACTCTTTTGTTTAAAAATTATATTCACAAGCACATCGGTACGATTTCTATCAACACGACTCTCGCGTTTTTATATAAAAATAAAATTTATGTGATGTTGGCTACCTATTGTTATTTTGTAATATGTGCGTTTATCGTCGGATATTCCCGCCCACATTTCAAAACAGTATTTGACTCATTAGAAATTGAAGGCAAAGCCTATTACGTATTAGCAAAATACAGTGACACTTTTATCTTAGCGAAAAGTACTCGTGCTACTAATGACAGTTTTTATCTTTATAAAATGAATATTAATTTTTTAT  >  1/1217055‑1217492

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