Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,610,096	T→C	N11S (AAT→AGT)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,610,096	0	T	C	100.0%	100.8 / NA	36	N11S (AAT→AGT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)

GCATAAACCTGTTCAAAATCCACCGTGCAGTCAGAACCGTCAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACATTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCACCTCCGGTTTAGGTCGTGCCTATGCCGGTGAAGCCGCCAT  >  1/1609884‑1610335                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                gCATAAACCTGTTCAAAATCCACCGTGCAGTCAGAACCGTCAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCcac                                                                                                                                                                                                                  >  1:433672/1‑244 (MQ=255) gCATAAACCTGTTCAAAATCCACCGTGCAGTCAGAACCGTCAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTccc                                                                                                                                                                                                            >  1:891143/1‑250 (MQ=255)              cAAAATCCACCGTGCAGTCAGAACCGTCAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAAt                                                                                                                                                                                              >  1:526849/1‑251 (MQ=255)                                       gTCAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCcac                                                                                                                                                                                                                  <  1:25888/206‑1 (MQ=255)                                         cAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAAcc                                                                                                                                                                   <  1:499769/251‑1 (MQ=255)                                                   aCCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATaa                                                                                                                                                         <  1:416094/251‑1 (MQ=255)                                                      cACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAg                                                                                                                                                      >  1:526507/1‑251 (MQ=255)                                                                          cTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAAcc                                                                                                                                                                   <  1:28104/218‑1 (MQ=255)                                                                          cTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACaaaaa                                                                                                                                             >  1:784968/1‑240 (MQ=255)                                                                                                         gCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTTCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACg                                                                                                       >  1:491888/1‑247 (MQ=255)                                                                                                            cATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACAt                                                                                                                                                                                                                <  1:1351/139‑1 (MQ=255)                                                                                                            cATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGttat                                                                                                >  1:655400/1‑251 (MQ=255)                                                                                                              taataTCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGttt                                                                                                                                                                                                      <  1:39397/147‑1 (MQ=255)                                                                                                              taataTCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGttatta                                                                                              >  1:800749/1‑251 (MQ=255)                                                                                                              taataTCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGttatta                                                                                              >  1:872074/1‑251 (MQ=255)                                                                                                                              aaaGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTg                                                                              <  1:713012/251‑1 (MQ=255)                                                                                                                                              cAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAAc                                                                                            >  1:239594/1‑221 (MQ=255)                                                                                                                                                   aTTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTcc                                                                                                                                                                                                             >  1:54672/1‑103 (MQ=255)                                                                                                                                                   aTTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATcc                                                                                                          >  1:608157/1‑202 (MQ=255)                                                                                                                                                   aTTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAAt                                                         >  1:548521/1‑251 (MQ=255)                                                                                                                                                             tattCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCaa                                                                                                                                                            <  1:28375/142‑1 (MQ=255)                                                                                                                                                              attCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGttt                                                                                                                                                                                     <  1:647399/116‑1 (MQ=255)                                                                                                                                                              attCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTc                                                                          <  1:10758/223‑1 (MQ=255)                                                                                                                                                              attCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAAtt                                              >  1:598828/1‑251 (MQ=255)                                                                                                                                                                   aaTTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCAc                                         <  1:375174/251‑1 (MQ=255)                                                                                                                                                                    aTTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCAc                                         >  1:238971/1‑250 (MQ=255)                                                                                                                                                                                           cTTTTCTTTGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACaaaaa                                                                                                                                             <  1:702825/127‑1 (MQ=255)                                                                                                                                                                                                  ttGTCCACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCACCTCCGGTTTAGGTCGTGCCTATGCCGGTGaa          >  1:881169/1‑251 (MQ=255)                                                                                                                                                                                                      ccACTGGTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCACCTCCGGTTTAGGTCGTGCCTATGCCGGTGAAgccg      >  1:126076/1‑251 (MQ=255)                                                                                                                                                                                                           ggTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACa                                                                                                                                        <  1:98838/116‑1 (MQ=255)                                                                                                                                                                                                           ggTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAAc                                                                                            >  1:596044/1‑160 (MQ=255)                                                                                                                                                                                                           ggTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGCGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCACCTCCGGTTTAGGTCGTGCCTATGCCGGTGAAGCCGCCAt  >  1:437475/1‑250 (MQ=255)                                                                                                                                                                                                           ggTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGc                                                                 >  1:352441/1‑187 (MQ=255)                                                                                                                                                                                                           ggTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAAtt                                                        <  1:25007/196‑1 (MQ=255)                                                                                                                                                                                                           ggTCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCACCTCCGGTTTAGGTCGTGCCTATGCCGGTGAAGCCGCCAt  >  1:632881/1‑250 (MQ=255)                                                                                                                                                                                                             tCTGAACACTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAAc                                                                                            <  1:117997/158‑1 (MQ=255)                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                GCATAAACCTGTTCAAAATCCACCGTGCAGTCAGAACCGTCAATCACCGTACCCACCTCGACACAACAGCAAACTTTTGCTTCTTCCGCCTGACATTTCAATGTGCTCATAATATCTCCTTGTAAAAAGAATTGAAATTTTCAACTATTTTAGCATTATTCAAATTTAATACAAAATCCGACCGCACTTTTCTTTGTCCACTGGTCTGAACATTACGCAAAATATGAAAAGAGGATTGCTCCACATTCCCCTGTTTTTTAGAATGATACGTTTTCCTAAGAAACGATAACCAATTCAATAAAAGGTACAAAAATGACAACCCGCAAATTTACTCAAACCCTTCTATCCACGCTGTTATTAACTGCTGCCGGAGGTGCGTCCGCCGCGGCATTCCAATTAGCAGAAATTTCCACCTCCGGTTTAGGTCGTGCCTATGCCGGTGAAGCCGCCAT  >  1/1609884‑1610335

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