Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,294,917	C→A	T137K (ACA→AAA)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,294,917	0	C	A	100.0%	123.9 / NA	36	T137K (ACA→AAA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (22/14);  total (22/14)

AATGAGTTGTGGCATAAGGTACGCCATTATCAGGGCTTTGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAACATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATCCTATTCAAAGAATTGTATTGCAAGATGAAATA  >  1/1294694‑1295083                                                                                                                                                                                                                                |                                                                                                                                                                       aaTGAGTTGTGGCATAAGGTACGCCATTATCAGGGCTTTGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCttt                                                                                                                                              >  1:501545/1‑250 (MQ=255)       ttGTGGCATAAGGTACGCCATTATCAGGGCTTTGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTaag                                                                                                                                        <  1:106563/250‑1 (MQ=255)            gCATAAGGTACGCCATTATCAGGGCTTTGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTg                                                                                                                                  >  1:522455/1‑251 (MQ=255)                            tATCAGGGCTTTGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAattagatt                                                                                                                   <  1:726123/250‑1 (MQ=255)                                     tttGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         <  1:895954/251‑1 (MQ=255)                                     tttGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         <  1:450807/251‑1 (MQ=255)                                                               ctctATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGc                                                                                 <  1:389903/249‑1 (MQ=255)                                                                        tagtagAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAat                                                                                          >  1:832168/1‑231 (MQ=255)                                                                                  aTTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAgg                                                                   >  1:10893/1‑244 (MQ=255)                                                                                                     cTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAg                                                                                              <  1:739141/198‑1 (MQ=255)                                                                                                               gATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAattagatt                                                                                                                   <  1:238341/167‑1 (MQ=255)                                                                                                                    atatAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         >  1:19918/1‑172 (MQ=255)                                                                                                                    atatAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         <  1:22530/172‑1 (MQ=255)                                                                                                                    atatAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         >  1:394549/1‑172 (MQ=255)                                                                                                                    atatAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAatatat                                                                                      >  1:725719/1‑191 (MQ=255)                                                                                                                                      cTTTTAATAAAGATGAATTCGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTc                                                                                                               >  1:707358/1‑148 (MQ=255)                                                                                                                                                                   gTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAattagat                                                                                                                    <  1:184787/114‑1 (MQ=255)                                                                                                                                                                   gTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATAttt                                                                                    >  1:683589/1‑146 (MQ=255)                                                                                                                                                                           tatTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCt                                                                                                              >  1:593617/1‑112 (MQ=255)                                                                                                                                                                                 taAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATCCTATTCAAAGAATTGTATTGCAAGATGaaata  >  1:756252/1‑214 (MQ=255)                                                                                                                                                                                       gCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         >  1:505543/1‑105 (MQ=255)                                                                                                                                                                                       gCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATAttt                                                                                    >  1:103608/1‑126 (MQ=255)                                                                                                                                                                                        cTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGc                                                                                 >  1:321351/1‑128 (MQ=255)                                                                                                                                                                                          aGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         >  1:155938/1‑102 (MQ=255)                                                                                                                                                                                          aGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAgg                                                                   <  1:109332/140‑1 (MQ=255)                                                                                                                                                                                          aGATTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATcc                                 >  1:211705/1‑174 (MQ=255)                                                                                                                                                                                            aTTACATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAAtt                                                        >  1:539977/1‑149 (MQ=255)                                                                                                                                                                                               aCATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCT‑‑‑‑TAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         >  1:279750/1‑93 (MQ=255)                                                                                                                                                                                                cATTCACCCATTCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGa                                              >  1:776770/1‑155 (MQ=255)                                                                                                                                                                                                          ttCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCt                                                                                                              >  1:41221/1‑81 (MQ=255)                                                                                                                                                                                                           tCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         <  1:407439/85‑1 (MQ=255)                                                                                                                                                                                                           tCCCTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATc                                  <  1:328534/156‑1 (MQ=255)                                                                                                                                                                                                            cccTGATGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAagc                                                                                                         >  1:819355/1‑84 (MQ=255)                                                                                                                                                                                                                  tGGGAATAGTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTc                                                                                                               >  1:467689/1‑72 (MQ=255)                                                                                                                                                                                                                          gTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATcc                                 <  1:329351/142‑1 (MQ=255)                                                                                                                                                                                                                          gTCGAAAATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATCCTATTCAAAGa                       <  1:490985/152‑1 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                       AATGAGTTGTGGCATAAGGTACGCCATTATCAGGGCTTTGGGAATATAACCGTTTTTTATAGCTCTATGGATAGTAGAAGTATTAAAGAGGCTGAAGATACTTTAAATGTGATCAATATAAACCAATTAAAAACTTTTAATAAAGATGAATTTGCTAAAGAAGTTGTAGATATTTATAAAAAGCTAGATTACATTCACCCATTCCCTGATGGGAATAGTCGAACATTACGAGAATTTACTCGAATTCTTTCGTAAGAAGTTGGGTTTAAATTAGATTAGTCTAAAGCAGCTCTAACAGAAATATATTTAGCAAGAGATTTTGAGGTAAATTCAATTAGTTTATTGAAAAATTCATATCCTATTCAAAGAATTGTATTGCAAGATGAAATA  >  1/1294694‑1295083

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: ‑