Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	819,874	T→G	G138G (GGA→GGC)	unknown ←	Molybdenum ABC transporter permease protein ModB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	819,874	0	T	G	100.0%	183.6 / NA	67	G138G (GGA→GGC)	unknown	Molybdenum ABC transporter permease protein ModB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (42/25);  total (42/25)

AAACCAGGGATAAAATTATAGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCTCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATACCAAAGGCAAATGGATGATGC  >  1/819630‑820118                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        aaaCCAGGGATAAAATTATAGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAg                                                                                                                                                                                                                                                   <  1:634705/251‑1 (MQ=255)  acccAGGGATAAAATTATAGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCCCCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAgg                                                                                                                                                                                                                                                  <  1:614196/249‑1 (MQ=255)  aaCCAGGGATAAAATTATAGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCCCCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAgg                                                                                                                                                                                                                                                  <  1:531245/251‑1 (MQ=255)                    aGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCa                                                                                                                                                                                                                                >  1:770311/1‑251 (MQ=255)                    aGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCa                                                                                                                                                                                                                                >  1:262238/1‑251 (MQ=255)                    aGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCa                                                                                                                                                                                                                                >  1:181822/1‑251 (MQ=255)                                tAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGa                                                                                                                                                                                                                    >  1:571279/1‑251 (MQ=255)                                        cggcTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTc                                                                                                                                                                                                             >  1:258711/1‑250 (MQ=255)                                            ttGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCg                                                                                                                                                                                                        <  1:326735/251‑1 (MQ=255)                                              gTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGaa                                                                                                                                                                                                      >  1:554319/1‑251 (MQ=255)                                                                   gaatgaaAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAGAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCa                                                                                                                                                                                   >  1:62633/1‑249 (MQ=255)                                                                        aaaGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCCCCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATg                                                                                                                                                                             <  1:654200/250‑1 (MQ=255)                                                                        aaaGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATgg                                                                                                                                                                            >  1:202414/1‑251 (MQ=255)                                                                          aGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGaa                                                                                                                                                                          >  1:538365/1‑251 (MQ=255)                                                                             aTACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATg                                                                                                                                                                       >  1:635356/1‑251 (MQ=255)                                                                               aaaTTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGcc                                                                                                                                                                     <  1:194987/249‑1 (MQ=255)                                                                                  ttGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGccac                                                                                                                                                                   >  1:151283/1‑250 (MQ=255)                                                                                    gCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAgg                                                                                                                                                                                                                                                  <  1:849218/169‑1 (MQ=255)                                                                                    gCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGt                                                                                                                                                                                                                                                 <  1:370149/170‑1 (MQ=255)                                                                                    gCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGccacca                                                                                                                                                                 >  1:486152/1‑250 (MQ=255)                                                                                       aaaGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACgg                                                                                                                                                              >  1:399011/1‑250 (MQ=255)                                                                                        aaGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGGAATGGTTGCGCCGAATTCCCCCAAGGAGCGGGCAAATCCCAATACGTTGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATg                                                                                                                                                                             <  1:323336/234‑1 (MQ=255)                                                                                          ggAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAg                                                                                                                                                           >  1:801548/1‑250 (MQ=255)                                                                                               aGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTa                                                                                                                                                     >  1:213397/1‑251 (MQ=255)                                                                                                   tgtgTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATAc                                                                                                                                                 <  1:141565/251‑1 (MQ=255)                                                                                                     tgtTTCACCGGGGATATTGGAAACGAAAGAAATGGTTGCGCCGAATTCACCAAGGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGGGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATAcc                                                                                                                                                <  1:464965/250‑1 (MQ=255)                                                                                                             cGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATAc                                                                                                                                                 >  1:458955/1‑241 (MQ=255)                                                                                                               gCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCa                                                                                                                                     >  1:170724/1‑251 (MQ=255)                                                                                                                                 gTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAg                                                                                                                   >  1:827306/1‑251 (MQ=255)                                                                                                                                                      cccAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAata                                                                                              >  1:322321/1‑251 (MQ=255)                                                                                                                                                       ccAAGGGGCGGGCAAATCCCAATACGAGGCCCGCTAAAATCCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAatat                                                                                             <  1:720689/251‑1 (MQ=255)                                                                                                                                                                 ggCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTggc                                                                                                                                                                                                                                           <  1:620383/99‑1 (MQ=255)                                                                                                                                                                   cAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGCCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGCCAGCCCGAACCATGGGAATAAATATTTGCCGATGaaa                                                                                 >  1:668713/1‑251 (MQ=255)                                                                                                                                                                      aTCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGTCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACTGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGCCAGCCCGAACCATGGGGATAAATATTTGCCGATGAAACCg                                                                              >  1:670542/1‑251 (MQ=255)                                                                                                                                                                      aTCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCg                                                                              >  1:267927/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATAc                                                                                                                                                                                                                 >  1:486739/1‑119 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAata                                                                                              >  1:491682/1‑234 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGt                                                                             >  1:529763/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGt                                                                             >  1:813467/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGt                                                                             >  1:770941/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGt                                                                             >  1:669814/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGt                                                                             >  1:568865/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGt                                                                             >  1:229401/1‑251 (MQ=255)                                                                                                                                                                       tCCCAATACGATGCCCGCAAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCTCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCg                                                                              <  1:428268/250‑1 (MQ=255)                                                                                                                                                                        cccAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAata                                                                                              >  1:671502/1‑233 (MQ=255)                                                                                                                                                                         ccAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGaa                                                                                                                                                                          <  1:760011/156‑1 (MQ=255)                                                                                                                                                                         ccAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTg                                                                           >  1:811669/1‑251 (MQ=255)                                                                                                                                                                                 aTGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCAt                                                                   <  1:317507/251‑1 (MQ=255)                                                                                                                                                                                                    gCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCt                                                                                                                                                      >  1:596235/1‑149 (MQ=255)                                                                                                                                                                                                         gATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAAt                                        >  1:273153/1‑251 (MQ=255)                                                                                                                                                                                                         gATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAAt                                        >  1:288865/1‑251 (MQ=255)                                                                                                                                                                                                         gATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAAt                                        >  1:320512/1‑251 (MQ=255)                                                                                                                                                                                                          aTAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAAt                                        >  1:614541/1‑250 (MQ=255)                                                                                                                                                                                                              gggCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAat                                                                                               <  1:147022/194‑1 (MQ=255)                                                                                                                                                                                                              gggCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGGAACCAATcacca                                   >  1:770452/1‑251 (MQ=255)                                                                                                                                                                                                                gCAAGGTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATcaccac                                  <  1:267299/250‑1 (MQ=255)                                                                                                                                                                                                                 aagggTGATGGTAAAAAATCCCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAg                                <  1:47190/248‑1 (MQ=255)                                                                                                                                                                                                                    ggTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATCCCACAACGCCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGAATAAat                                                                                               >  1:64894/1‑188 (MQ=255)                                                                                                                                                                                                                     gTGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAgggg                             <  1:725385/250‑1 (MQ=255)                                                                                                                                                                                                                      tGATGGTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAggggg                            <  1:114633/250‑1 (MQ=255)                                                                                                                                                                                                                           gTAAAAAATACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATa                       >  1:792527/1‑250 (MQ=255)                                                                                                                                                                                                                                 aaTACCCGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATACCAAAgg                <  1:306530/251‑1 (MQ=255)                                                                                                                                                                                                                                       cGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCACCAACAACAACAGGTAACCAATCACCACAGGGGGCAATACCAAAGGCaa             >  1:139631/1‑251 (MQ=255)                                                                                                                                                                                                                                       cGCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATAc                      <  1:308220/239‑1 (MQ=255)                                                                                                                                                                                                                                        gCCAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGc                           >  1:574124/1‑233 (MQ=255)                                                                                                                                                                                                                                         ccAAGGGGAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATGTTAAGTCGATACTGTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATACCAAAGGCAAATgg         <  1:745447/250‑1 (MQ=255)                                                                                                                                                                                                                                                gAAGCGCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATACCAAAGGCAAATGGATGATGc  <  1:723136/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        AAACCAGGGATAAAATTATAGCAAAAACACATAAACGGGCGGCTTGTTCTTCCGCGCCGGGCGTTTGAATGAAAGAATACATTGCTAAAGGAATAGTTTGTGTTTCACCGGCGATATTGGAAACGAAAGTAATGGTTGCGCCGAATTCACCCAAGGAGCGGGCAAATCCCAATACGATGCCCGCTAAAATTCCCGGCAAAGATAGGGGCAAGGTGATGGTAAAAAATACCCGCCAAGGGGAAGCTCCCAAGGTTCTGGCGGCTTGTTCCAATTTTAAGTCGATACTTTCCAGCGAAAGACGAATTGCCCGAACCACTAATGGAAATGCCACCACGGCAGCGGCTAATACCGCCCCGCGCCAGCTGAAGCCGAAAGTCAGCCCGAACCATTGGTATAAATATTTGCCGATGAAACCGTTGCGTCCCATTGCCAC‑‑‑CAACAACAGGTAACCAATCACCACAGGGGGCAATACCAAAGGCAAATGGATGATGC  >  1/819630‑820118

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