Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,819,986	T→C	intergenic (‑51/+6)	unknown ← / ← unknown	F200/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,819,986	0	T	C	100.0%	110.2 / NA	39	intergenic (‑51/+6)	unknown/unknown	F200/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/17);  total (22/17)

ACTTAGCCATGTTGGAGATATTAAAGGTTTTAAGGCCTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCTGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTACGCCGGTTAATAAATTAATTTTTTTCATGATAAACTCCTTA  >  1/1819780‑1820214                                                                                                                                                                                                               |                                                                                                                                                                                                                                     aCTTAGCCATGTTGGAGATATTAAAGGTTTTAAGGCCTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAAc                                                                                                                                                                                          <  1:197760/251‑1 (MQ=255)   ttAGCCATGTTGGAGATATTAAAGGTTTTAAGGCCTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATAAGTCTATCATACTCTCTAACAg                                                                                                                                                                                        <  1:552418/251‑1 (MQ=255)                           gTTTTAAGGCCTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAg                                                                                                                                                                                        >  1:659854/1‑227 (MQ=255)                           gTTTTAAGGCCTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACtt                                                                                                                                                                 >  1:126463/1‑250 (MQ=255)                                    ccTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATActc                                                                                                                                                                                                <  1:711470/210‑1 (MQ=255)                                                  tCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAg                                                                                                                                         <  1:311369/251‑1 (MQ=255)                                                                                ataataAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATg                                                                                                            <  1:163159/250‑1 (MQ=255)                                                                                       aaaCAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTacac                                                                                                    <  1:281962/251‑1 (MQ=255)                                                                                          caacTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACActtt                                                                                                 >  1:655918/1‑251 (MQ=255)                                                                                                   aTCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGttt                                                                                                                           <  1:211857/216‑1 (MQ=255)                                                                                                    tCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATAc                                                                                                                                                                                                  <  1:654051/144‑1 (MQ=255)                                                                                                    tCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAAc                                                                                                                                                                                          >  1:460138/1‑152 (MQ=255)                                                                                                                                                                    ctcCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCt                                                                                                                                       >  1:829881/1‑139 (MQ=255)                                                                                                                                                                        aTGTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCt                                                                          <  1:681152/196‑1 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCa                                                                                                                                                 >  1:800304/1‑123 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATg                                                                                                                                               <  1:124966/125‑1 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGc                                                                                                                                        <  1:555637/132‑1 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTa                                                                                                                                      <  1:893055/134‑1 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGc                                                >  1:552674/1‑220 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTa                                           >  1:160377/1‑225 (MQ=255)                                                                                                                                                                          gTATAATAATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTACGCCgg                                     >  1:315440/1‑231 (MQ=255)                                                                                                                                                                          gTATAATAATGACTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTACGCCGGTTAATAAATTAAttttttt                  >  1:240282/1‑250 (MQ=255)                                                                                                                                                                            ataataATGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAatat                                                                                                                      <  1:272971/148‑1 (MQ=255)                                                                                                                                                                             taataaTGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCAt                                                                                                                                                                                                    >  1:233182/1‑69 (MQ=255)                                                                                                                                                                             taataaTGGCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTACGCCGGTTAATAAATTAATTTTTTTCATg              >  1:561101/1‑251 (MQ=255)                                                                                                                                                                               ataatGGCTTATCTTATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATc                                                                                                                                                         >  1:811352/1‑110 (MQ=255)                                                                                                                                                                                    ggCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATAc                                                                                                                                                                                                  >  1:348474/1‑64 (MQ=255)                                                                                                                                                                                     gCTTATCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCt                                                                                                                                       >  1:642039/1‑122 (MQ=255)                                                                                                                                                                                         aTCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACActttct                                                                                               <  1:205458/158‑1 (MQ=255)                                                                                                                                                                                         aTCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCatat                                                                                        >  1:219276/1‑165 (MQ=255)                                                                                                                                                                                         aTCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGc                                                <  1:366818/205‑1 (MQ=255)                                                                                                                                                                                         aTCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTa                                           <  1:716997/210‑1 (MQ=255)                                                                                                                                                                                         aTCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTACGCCGGTTAATAAATTAATTTTTTTCATGATAAACTCCTTa  >  1:675482/1‑251 (MQ=255)                                                                                                                                                                                          tCATATTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCa                                                                                                                                                                                                     >  1:552139/1‑55 (MQ=255)                                                                                                                                                                                              atTCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTa                                           >  1:455515/1‑205 (MQ=255)                                                                                                                                                                                               ttCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGAt                                                                                                                                                                            >  1:135757/1‑75 (MQ=255)                                                                                                                                                                                               ttCTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATAc                                                                                                                                                                   >  1:239221/1‑84 (MQ=255)                                                                                                                                                                                                 cTATGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGc                                                                                                                                        <  1:336256/109‑1 (MQ=255)                                                                                                                                                                                                   aTGATAAGCCTCCGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTcc                                                                                                                                                   >  1:380147/1‑96 (MQ=255)                                                                                                                                                                                                               |                                                                                                                                                                                                                                     ACTTAGCCATGTTGGAGATATTAAAGGTTTTAAGGCCTCTACCTCAAAATCAAACCATTTATAGTTACCAAATAATGCGATAATAAAAACAACTGAAAATCTTAATACGGTAATATCAATATTAGAAAGGCAAAACTTTCTAAGAACTTCTTTCATATAAATTCTCCATGTATAATAATGGCTTATCATATTCTATGATAAGCCTCTGAAAATTATTTATCCCCTTTCATACGTCTATCATACTCTCTAACAGTCGGTTTATGATAATTCATACTTCCCAGATCATTTCCCATGTTATAGCTAATCTCTTGTTTAATATCGTCTGAATGCATTACACTTTCTTTCATATTTGTTGTGCTAGCTGATGCTAATGATGCAATAGTAAGAGCAAGTACGCCGGTTAATAAATTAATTTTTTTCATGATAAACTCCTTA  >  1/1819780‑1820214

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