Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	704,881	A→C	F31L (TTT→TTG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	704,881	0	A	C	100.0%	78.7 / NA	28	F31L (TTT→TTG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (16/12);  total (16/12)

CGGCAAAAGACAATCAGAAAAGATCCGACTAACGGATGCTGAAAACGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCAAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACACACCAAATGGAACTGTAAGTACCACATTGTTTTCACCCT  >  1/704634‑704985                                                                                                                                                                                                                                                        |                                                                                                         cGGCAAAAGACAATCAGAAAAGATCCGACTAACGGATGCTGAAAACGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCcaaa                                                                                                       <  1:838338/251‑1 (MQ=255)                        tCCGACTAACGGATGCTGAAAACGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTc                                                                                <  1:898334/251‑1 (MQ=255)                                       cTGAAAACGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGc                                                                 <  1:130742/251‑1 (MQ=255)                                          aaaaCGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCa                                                               <  1:859212/250‑1 (MQ=255)                                              cGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGAc                                                          <  1:517810/251‑1 (MQ=255)                                                aaaCAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGa                                                        >  1:324282/1‑251 (MQ=255)                                                         tGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGt                                                 <  1:881643/249‑1 (MQ=255)                                                             gTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGAt                                                                                     >  1:19153/1‑209 (MQ=255)                                                             gTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATca                                           >  1:820650/1‑251 (MQ=255)                                                               cATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATaggag                                                                                        <  1:492080/204‑1 (MQ=255)                                                                aTTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATcacat                                        >  1:25034/1‑250 (MQ=255)                                                                 ttATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATagg                                                                                          >  1:803052/1‑200 (MQ=255)                                                                 ttATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATTGTAAGTAAAGCCAATGACGATTCAAGTCTATcacata                                       >  1:837685/1‑249 (MQ=255)                                                                 ttATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGa                                                           >  1:4469/1‑231 (MQ=255)                                                                 ttATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATcacata                                       >  1:446888/1‑249 (MQ=255)                                                                 ttATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAAAAAAACCGTACCCTATGCATTTGGAATCAAAAAGCTTAAGCGGTCAACCTAAAAAATCCTCACTATATAAAAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGAAAAGCCAATGACGATTCAAGTCTATcccata                                       >  1:829511/1‑246 (MQ=18)                                                                  tataGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATcacatac                                      <  1:503469/251‑4 (MQ=255)                                                                   ataGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATcacatac                                      >  1:153259/1‑247 (MQ=255)                                                                    taGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTc                                                                                >  1:646567/1‑207 (MQ=255)                                                                          gTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGa                              >  1:774111/1‑251 (MQ=255)                                                                            aTTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAAc                            <  1:48873/251‑1 (MQ=255)                                                                            aTTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAAc                            <  1:199794/251‑1 (MQ=255)                                                                               tttCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTTTAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATAAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAAAAAGGCTGAAAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTTTATCACATACCAAATGGaa                             >  1:262171/1‑247 (MQ=17)                                                                                          cGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAACTGTAAGTATCACAt              <  1:681270/251‑1 (MQ=255)                                                                                          cGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAACTGTAAGTATCACAt              <  1:538243/251‑1 (MQ=255)                                                                                                      aTAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAACTGTAAGTATCACAtcgtcttcccccc  <  1:689041/251‑13 (MQ=255)                                                                                                      aTAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAACTGTAAGTATCACATCGTCTTCAccc   >  1:482792/1‑250 (MQ=255)                                                                                                      aTAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAACTGTAAGTATCACATCGTCTTCAccc   >  1:229543/1‑250 (MQ=255)                                                                                                      aTAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCCAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACATACCAAATGGAACTGTAAGTATCACATCGTCTTCAccc   >  1:165471/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                         CGGCAAAAGACAATCAGAAAAGATCCGACTAACGGATGCTGAAAACGAAACAAAAGTGCGGTCATTATAGAGGGTATTTTTCCATTTAGCGAATCATGCCAATAAAAACCCGAGACTGGCTCGGGTTTTTGTTTATAGCTAGACAACGTCTAGCTATAAACAAAACCGTACGCTATGCATATGGAATCAAAGAGCTTAAGCGGTCAACATAAAAAATCCTCACTATATAAGAATAAGGCTGACAACCAAAATTCAAGATATAGGAGGATAAGTCAATGGTAAGTAAAGCCAATGACGATTCAAGTCTATCACACACCAAATGGAACTGTAAGTACCACATTGTTTTCACCCT  >  1/704634‑704985

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