Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 886,434 T→C H43H (CAT→CAC)  unknown → tRNA cytosine(34) acetyltransferase (EC 2.3.1.193)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1886,4340TC100.0% 194.4 / NA 66H43H (CAT→CAC) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (46/20);  total (46/20)

CAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGCGAACCGCTCTATAAAAGCAATAAAAATATCAATGTAGACTGGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCC  >  1/886195‑886668
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          
cAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAg                                                                                                                                                                                                                                 >  1:572754/1‑251 (MQ=255)
  aCGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCa                                                                                                                                                                                                                               >  1:198051/1‑251 (MQ=255)
             ttGAACGCCAAATTCGCGGCGAAAGTGAACCGCTCTATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAAt                                                                                                                                                                                                                    <  1:108532/251‑1 (MQ=255)
                                             ctctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCTCGCGCTTATTTa                                                                                                                                                                                    >  1:467338/1‑251 (MQ=255)
                                               ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCttatttat                                                                                                                                                                                   >  1:154465/1‑250 (MQ=255)
                                               ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatg                                                                                                                                                                                  >  1:11930/1‑251 (MQ=255)
                                               ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatg                                                                                                                                                                                  >  1:590757/1‑251 (MQ=255)
                                               ctATAAAAGCAATAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatg                                                                                                                                                                                  >  1:823761/1‑251 (MQ=255)
                                                          aTAAAAATATTAATGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGgcgtcag                                                                                                                                                                       >  1:160657/1‑247 (MQ=255)
                                                                 tattaaTGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTccc                                                                                                                                                                                                >  1:577465/5‑219 (MQ=255)
                                                                     aaTGTAGATTGGACGGATCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGgcgtcag                                                                                                                                                                       >  1:558484/1‑236 (MQ=255)
                                                                               ggacggatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGctc                                                                                                                                                    >  1:28446/5‑247 (MQ=255)
                                                                                    gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatga                                                                                                                                                                                 <  1:373470/215‑1 (MQ=255)
                                                                                    gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATg                                                                                                                                                       >  1:301001/1‑241 (MQ=255)
                                                                                    gatCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGccat                                                                                                                                             >  1:238826/1‑248 (MQ=255)
                                                                                      tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAGGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatga                                                                                                                                                                                 >  1:437524/1‑213 (MQ=255)
                                                                                      tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTc                                                                                                                                                                                                  >  1:647836/1‑196 (MQ=255)
                                                                                      tCCGTATTCACGGGGTTGGTAAACTATGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTg                                                                                                                                           >  1:179552/1‑251 (MQ=255)
                                                                                                          aaactaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGAt                                                                                                                                                        >  1:492585/6‑218 (MQ=255)
                                                                                                          aaactaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGaa                                                                                                                       <  1:248402/246‑1 (MQ=255)
                                                                                                            actaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGcc                                                                                                                                               >  1:592869/4‑224 (MQ=255)
                                                                                                             ctaTGGCGCGGCAGTTAATTCTTCTCATTGGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCAGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAgg                                                                                                                     >  1:212219/3‑250 (MQ=255)
                                                                                                                        cagttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                          <  1:300262/234‑1 (MQ=255)
                                                                                                                        cagttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                          <  1:193059/234‑1 (MQ=255)
                                                                                                                         agttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgttg                                                                                                        >  1:470535/16‑251 (MQ=255)
                                                                                                                         agttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgttg                                                                                                        >  1:870287/16‑251 (MQ=255)
                                                                                                                          gttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                          <  1:271834/234‑1 (MQ=255)
                                                                                                                           ttaattcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgttg                                                                                                        <  1:56150/236‑1 (MQ=255)
                                                                                                                                tcttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTAt                                                                                                 >  1:726572/9‑251 (MQ=255)
                                                                                                                                  ttctcattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTc                                                                                               <  1:104009/245‑1 (MQ=255)
                                                                                                                                      cattgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTatga                                                                                                                                                                                 >  1:23565/3‑165 (MQ=255)
                                                                                                                                       attgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGa                                                                                          <  1:867890/250‑1 (MQ=255)
                                                                                                                                        ttgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGat                                                                                         >  1:307769/1‑250 (MQ=255)
                                                                                                                                        ttgGTGACGATGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGa                                                                                          <  1:14211/250‑1 (MQ=255)
                                                                                                                                          ggtgacgaTGCCACGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATAttgt                                                                                                          >  1:411686/4‑232 (MQ=255)
                                                                                                                                                      aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggc                                                                           <  1:29227/251‑1 (MQ=255)
                                                                                                                                                      aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTgg                                                                            >  1:59991/1‑250 (MQ=255)
                                                                                                                                                      aCGTTAGCGCACCAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTgg                                                                            >  1:630816/1‑250 (MQ=255)
                                                                                                                                                                 ccAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaac                                                                   >  1:2743/1‑242 (MQ=255)
                                                                                                                                                                 ccAATGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaac                                                                   >  1:756796/1‑242 (MQ=255)
                                                                                                                                                                     tGGGAACATTTTCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatgg                                                            <  1:63206/251‑14 (MQ=255)
                                                                                                                                                                      gggaacattttCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatggg                                                           >  1:253168/9‑237 (MQ=255)
                                                                                                                                                                          acattttCCTGCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacctgatggtgaca                                                       >  1:309866/5‑233 (MQ=255)
                                                                                                                                                                                    gCCCTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGAt                                                                                                                                                        >  1:27031/1‑144 (MQ=255)
                                                                                                                                                                                     cccTGAAAAGTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGt                                                                                                                            >  1:71496/1‑171 (MQ=255)
                                                                                                                                                                                              gTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGa                                                                                                                                                                                                      >  1:575030/1‑88 (MQ=255)
                                                                                                                                                                                              gTGCGGTATGGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTGTTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTggcggtgcaacc                                                                  <  1:893088/220‑8 (MQ=255)
                                                                                                                                                                                                   gtgaaacgccgccaatccataccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGc                                         >  1:335772/22‑240 (MQ=255)
                                                                                                                                                                                                   gtatgGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAaca                              >  1:235089/4‑251 (MQ=255)
                                                                                                                                                                                                     atgGATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGa                                                                                                                        >  1:287215/2‑159 (MQ=255)
                                                                                                                                                                                                       ggATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCg                          >  1:212676/1‑251 (MQ=255)
                                                                                                                                                                                                        gATTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGa                         >  1:210987/1‑251 (MQ=255)
                                                                                                                                                                                                         aTTGGCGGCGTTTCACCGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGAt                        <  1:642971/251‑1 (MQ=255)
                                                                                                                                                                                                              cggcgtttcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAac                    <  1:719897/240‑1 (MQ=255)
                                                                                                                                                                                                                gcgtttcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAac                    >  1:71257/9‑248 (MQ=255)
                                                                                                                                                                                                                 cgtttcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAacacac                <  1:409167/244‑1 (MQ=255)
                                                                                                                                                                                                                     tcaccGCACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCaa                                                                                             >  1:382846/4‑170 (MQ=255)
                                                                                                                                                                                                                           cACTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTAt      <  1:37136/251‑1 (MQ=255)
                                                                                                                                                                                                                             cTTTTCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTg    >  1:436916/1‑251 (MQ=255)
                                                                                                                                                                                                                              ttttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   <  1:322993/251‑1 (MQ=255)
                                                                                                                                                                                                                              ttttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   <  1:25093/251‑1 (MQ=255)
                                                                                                                                                                                                                               tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGATCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   >  1:546230/1‑250 (MQ=255)
                                                                                                                                                                                                                               tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGcc  <  1:163458/251‑1 (MQ=255)
                                                                                                                                                                                                                               tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   >  1:239994/1‑250 (MQ=255)
                                                                                                                                                                                                                               tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   >  1:552108/1‑250 (MQ=255)
                                                                                                                                                                                                                               tttCTGCTCAGCCTCACTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTTATTTATGATGGGCGTCAGGGCATTCATTTAGATGCTCTCGCCATTGCCGCCGGTGCGCTGTCGGAAGGCGGCATATTGTTGCTGTTATTCAATGATTGGCAGGATTTGGCGGTGCAACCTGATGGTGACAGCTTGCGTTGGAGCGGTGAAAAACACCCGATAAACACACCGCACTTTATTGc   >  1:151475/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          
CAACGAGAATCTATTGAACGCCAAATTCGCGGCGAAAGCGAACCGCTCTATAAAAGCAATAAAAATATCAATGTAGACTGGATGGATCCGTATTCGCGGGGTTGGTAAAGCATGGCGCGGCAGCTAATTCTCCTTGTTGGCGACGATGCCACGTTAGCGCGCCAATGGGAATACTTTCCTGCCCTGAAAAGTGCGGTGTGGATTGGCGGTGTTTTGCCGCACTTTTCCGCTCAACTTCATTTCCCTTTTAGCAAAAGTAAAAATCTGCTCGGCAGCGAATTCCCCGCGCTCATTTATGATGGGCGCCAAGGCATTCATTTAGATGCCCTCGCTATTGCCGCCGGTGCGTTGTCGGAAGGCGGCACATTGTTGCTGTTATTCAATGAGTGGCAGGATTTGGCGGCACAACCCGACGGCGACAGCCTGCGTTGGAGCAGTGAAAAACACCCGATAAACACACCGCACTTTATTGCC  >  1/886195‑886668

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