Predicted mutation
evidence seq id position mutation annotation gene description
RA 5 15,914 T→C I274M (ATA→ATG)  unknown ← 1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase (EC 1.1.1.267)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*515,9140TC100.0% 156.2 / NA 54I274M (ATA→ATG) unknown1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase (EC 1.1.1.267)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/45);  total (9/45)

AATTTAATTTGGCGATCTAAAAATGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCTATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTCCGGCGTAATGCCGTCAAATTCGCTTAACGCCGTAT  >  5/15669‑16152
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
aatttaatttGGCGATCTAAAAATGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTg                                                                                                                                                                                                                                           <  1:81322/251‑1 (MQ=255)
        ttGGCGATCTAAAAATGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGa                                                                                                                                                                                                                                   <  1:835840/251‑1 (MQ=255)
          ggCGATCTAAAAATGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACg                                                                                                                                                                                                                                 >  1:534936/1‑251 (MQ=255)
                     aaTGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTg                                                                                                                                                                                                                                           >  1:185316/1‑230 (MQ=255)
                     aaTGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTCTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGat                                                                                                                                                                                                                      <  1:450079/251‑1 (MQ=255)
                             cacaGCGATTTCATAGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCGAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGc                                                                                                                                                                                                               <  1:636572/250‑1 (MQ=255)
                                  cGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCAt                                                                                                                                                                                                          <  1:296346/250‑1 (MQ=255)
                                      ttCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAgaat                                                                                                                                                                                                     <  1:349221/251‑1 (MQ=255)
                                             gcggcATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGc                                                                                                                                                                                                               <  1:435782/234‑1 (MQ=255)
                                                gcATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCTGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATgg                                                                                                                                                                                            <  1:428208/250‑1 (MQ=255)
                                                gcATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATgg                                                                                                                                                                                            <  1:610809/250‑1 (MQ=255)
                                                gcATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATgg                                                                                                                                                                                            <  1:324166/250‑1 (MQ=255)
                                                 cattcattGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGAt                                                                                                                                                                                          <  1:344521/251‑1 (MQ=255)
                                                  attcattGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCAcc                                                                                                                                                                                                            <  1:529642/232‑1 (MQ=255)
                                                  attcattGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGAtt                                                                                                                                                                                         >  1:299209/1‑251 (MQ=255)
                                                   ttcattGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGaatg                                                                                                                                                                                        <  1:778252/251‑4 (MQ=255)
                                                   ttcattGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTg                                                                                                                                                                                        <  1:778266/251‑1 (MQ=255)
                                                            gtcgtAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAgaa                                                                                                                                                                                                      >  1:59744/1‑228 (MQ=255)
                                                             tcgtAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGaat                                                                                                                                                                               <  1:396350/250‑1 (MQ=255)
                                                             tcgtAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGGGTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGaat                                                                                                                                                                               <  1:770111/250‑1 (MQ=255)
                                                                   cATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGaat                                                                                                                                                                               >  1:223095/1‑244 (MQ=255)
                                                                    atatTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATg                                                                                                                                                                                             <  1:158206/229‑1 (MQ=255)
                                                                         gACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCAtt                                                                                                                                                                  >  1:95939/1‑251 (MQ=255)
                                                                                       gCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGat                                                                                                                                                                                                                      <  1:626216/185‑1 (MQ=255)
                                                                                                gCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAaaca                                                                                                                                            <  1:578763/250‑1 (MQ=255)
                                                                                                 cGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAaacaac                                                                                                                                          <  1:670003/251‑1 (MQ=255)
                                                                                                    gTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGc                                                                                                                                                                                                               >  1:648580/1‑179 (MQ=255)
                                                                                                               ggATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTg                                                                                                                                                                                                                                           <  1:270900/140‑1 (MQ=255)
                                                                                                               ggATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCg                                                                                                                                                                                                                           <  1:21466/156‑1 (MQ=255)
                                                                                                               ggATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAAc                                                                                                                                                                                                                 <  1:399139/166‑1 (MQ=255)
                                                                                                                 aTAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAatat                                                                                                                          <  1:524920/251‑1 (MQ=255)
                                                                                                                 aTAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAatat                                                                                                                          <  1:517753/251‑1 (MQ=255)
                                                                                                                     cGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTc                                                                                                                      <  1:283215/251‑1 (MQ=255)
                                                                                                                       ttGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGAt                                                                                                                                                                                          <  1:34369/181‑1 (MQ=255)
                                                                                                                       ttGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAatat                                                                                                                          <  1:681978/245‑1 (MQ=255)
                                                                                                                        tGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGatat                                                                                                                                                                                                                    >  1:643726/1‑154 (MQ=255)
                                                                                                                           tAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGcc                                                                                                                <  1:607007/251‑1 (MQ=255)
                                                                                                                           tAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGcc                                                                                                                <  1:373130/251‑1 (MQ=255)
                                                                                                                                cGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGt                                                                                                           <  1:624140/251‑1 (MQ=255)
                                                                                                                                    gCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCAt                                                                                                                                                                                                          <  1:761840/152‑1 (MQ=255)
                                                                                                                                       aaaaaaGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTcatcat                                                                                                    <  1:621152/251‑1 (MQ=255)
                                                                                                                                         aaaaGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCAtgg                                                                                                  <  1:658579/251‑1 (MQ=255)
                                                                                                                                                     gTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGAt                                                                                                                                                                                          <  1:694118/151‑1 (MQ=255)
                                                                                                                                                                                 aCGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGcc                                                                                                                >  1:677874/1‑197 (MQ=255)
                                                                                                                                                                                       gAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTc                                                                                                                      <  1:618960/185‑1 (MQ=255)
                                                                                                                                                                                             gtgcggtgcgGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGAcc                                                              <  1:512598/235‑1 (MQ=255)
                                                                                                                                                                                                  gtgcggATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATg                                                    <  1:807257/240‑1 (MQ=255)
                                                                                                                                                                                                      ggATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGc                                           <  1:223116/245‑1 (MQ=255)
                                                                                                                                                                                                       gATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTc                                     <  1:736114/250‑1 (MQ=255)
                                                                                                                                                                                                                  gTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTCCGGCGTAATGc                          <  1:788969/250‑1 (MQ=255)
                                                                                                                                                                                                                                    cgcATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTCCGGCGTAATGCCGTCAAATTCGCTTAACGc       <  1:635290/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cgcATATCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTCCGGCGTAATGCCGTCAAATTCGCTTAACGc       <  1:781964/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tatCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAaaca                                                                                                                                            <  1:477412/114‑1 (MQ=255)
                                                                                                                                                                                                                                         atCAGGATTGCCCATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTCCGGCGTAATGCCGTCAAATTCGCTTAACGCCGtat  <  1:683935/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              
AATTTAATTTGGCGATCTAAAAATGCCTGCACAGCGATTTCATTGGCGGCATTCATTGCCGTCGTAGCATATTGACCGGCGGCGAACGCCTCGATGGCGAGTTTCAAGCAAGGATAACGTTGATAGTCCGGTGCCAAAAAAGTCAGCCCGTTTAATTGATAGAAATCTAGCGGCGCGACGCCGGAAAACGTGCGGTGCGGATACGCCATGGTTTCGGCGATCGGCGTGCGCATATCAGGATTGCCTATTTGCGCGATGACGCTGCCGTCGATATAACGCACCATAGAATGAATAATGGATTGCGGGTGAATAATCACTTCCATTTCGTCGGCGGACACGTTAAACAACCAGCGGGCTTCAATATATTCCAAGCCTTTGTTCATCATGGTGGCGGAATCCACGGAAATTTTCTTGCCCATAGACCAGTTCGGATGCGCTACCGCTTGTTCCGGCGTAATGCCGTCAAATTCGCTTAACGCCGTAT  >  5/15669‑16152

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