Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 887,489 A→G Q395R (CAG→CGG)  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
*1887,4890AG100.0% 176.6 / NA 61Q395R (CAG→CGG) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/35);  total (26/35)

TGACAAACTGGAGCGCTTTATTGATGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGACGCAACCGCGTTACGATGCCAACCTTGCCAAACAGATTCAGATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCG  >  1/887247‑887719
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
tGACAAACTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCgggaaacc                                                                                                                                                                                                                                 >  1:428935/1‑246 (MQ=255)
     aaCTGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCgggaaaccctgttg                                                                                                                                                                                                                           <  1:604455/251‑11 (MQ=255)
        tGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                                                                                                                                        <  1:527982/251‑2 (MQ=255)
        tGGAGCGCTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                                                                                                                                        <  1:367894/251‑2 (MQ=255)
               cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                                                                                                                                        <  1:644225/244‑2 (MQ=255)
               cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCtgtgtgt                                                                                                                                                                                                                 <  1:39237/251‑1 (MQ=255)
               cTTTATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCtgtgtg                                                                                                                                                                                                                  <  1:420895/250‑1 (MQ=255)
                  tATTGACGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTggg                                                                                                                                                                                                              <  1:278758/251‑1 (MQ=255)
                     tgacgaATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGcgg                                                                                                                                                                                                           <  1:607949/247‑3 (MQ=255)
                      gacgaATTGTTGTTGTGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggctg                                                                                                                                                                                                                      >  1:52133/4‑236 (MQ=255)
                             tgttgttgtGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaa                                                                                                                                                                                                   <  1:175357/251‑2 (MQ=255)
                               ttgttgtGTGATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaag                                                                                                                                                                                                  <  1:619874/250‑3 (MQ=255)
                                    gtgtgATTGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggat                                                                                                                                                                                            >  1:406502/1‑243 (MQ=255)
                                          ttGCGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaag                                                                                                                                                                                                  >  1:342096/1‑237 (MQ=255)
                                             cGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCTTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggatttgaggat                                                                                                                                                                                    <  1:443662/250‑17 (MQ=255)
                                             cGAAGATCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggatttgaggat                                                                                                                                                                                    <  1:674652/250‑17 (MQ=255)
                                                  aTCGTCTGCCGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaag                                                                                                                                                                                                  >  1:823891/1‑229 (MQ=255)
                                                        tgccGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                        <  1:819051/248‑1 (MQ=255)
                                                         gccGCAACCGTGTTATGACGCTAACCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                        >  1:408787/3‑250 (MQ=255)
                                                                  gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGCTTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAAtt                                                                                                                                                               >  1:890573/14‑250 (MQ=255)
                                                                  gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaa                                                                                                                                                                                                   <  1:252535/201‑2 (MQ=255)
                                                                  gtgttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAAtt                                                                                                                                                               >  1:283929/14‑250 (MQ=255)
                                                                    gttatgacgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCga                                                                                                                                                            >  1:177839/12‑250 (MQ=255)
                                                                          acgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCtgtgtg                                                                                                                                                                                                                  >  1:290264/6‑191 (MQ=255)
                                                                          acgctaaCCTTGCCGAACAGATTCAGATTTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaaggaggatttg                                                                                                                                                                                         <  1:695858/211‑12 (MQ=255)
                                                                                             gattcagattTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggctg                                                                                                                                                                                                                      <  1:735718/168‑4 (MQ=255)
                                                                                                   gattTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                        >  1:77889/1‑208 (MQ=255)
                                                                                                    attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                                                                                                                             <  1:130677/250‑18 (MQ=255)
                                                                                                    attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                                                                                                                             <  1:28453/250‑18 (MQ=255)
                                                                                                           acaaagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTgggc                                                                                                                                                                                                                        >  1:81350/4‑151 (MQ=255)
                                                                                                            caaagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTc                                                                                                                                                              >  1:803115/3‑209 (MQ=255)
                                                                                                              aagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                  <  1:798858/251‑1 (MQ=255)
                                                                                                                gaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAAtt                                                                                                                                                                        <  1:775515/195‑1 (MQ=255)
                                                                                                                 aacggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtgttg                                                                                                               <  1:809807/247‑1 (MQ=255)
                                                                                                                   cggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                             <  1:390128/249‑1 (MQ=255)
                                                                                                                                         aatcgGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctg                                                                                                                                          >  1:114989/1‑196 (MQ=255)
                                                                                                                                            cgGTTTATGGCTTGCTGACTTTGGCGCATTACCTCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGGCGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCg                                                                                    <  1:506287/251‑1 (MQ=255)
                                                                                                                                             ggTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtc                                                                                                                                               >  1:739850/1‑191 (MQ=255)
                                                                                                                                              gTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCgggaaaccctgttg                                                                                                                                                                                                                           >  1:488911/1‑104 (MQ=255)
                                                                                                                                              gTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGaagg                                                                                                                                                                                                 <  1:229874/140‑4 (MQ=255)
                                                                                                                                              gTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                         <  1:882176/244‑1 (MQ=255)
                                                                                                                                                 tATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                               <  1:733341/251‑1 (MQ=255)
                                                                                                                                                       ttGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGcgg                                                                                                                                                                                                           <  1:685081/121‑3 (MQ=255)
                                                                                                                                                         gCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                  >  1:390954/1‑208 (MQ=255)
                                                                                                                                                                        ttacagcaCCTCGCCGTTGGATTTACGACGATTATTGGATGCACCGAGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGAGTGTGGGCGGAGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGGAAATTGGGTCCTTTGCGCATTTcc                                                         >  1:210964/8‑250 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCgggaaaccc                                                                                                                                                                                                                                <  1:468169/83‑6 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGAc                                                                                                                                                     <  1:97440/158‑1 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                             <  1:701921/198‑1 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                         >  1:728742/1‑218 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                         >  1:698893/1‑250 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                         >  1:51847/1‑250 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                         >  1:209876/1‑250 (MQ=255)
                                                                                                                                                                        ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                         <  1:582402/250‑1 (MQ=255)
                                                                                                                                                                          aCCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                               >  1:632915/1‑226 (MQ=255)
                                                                                                                                                                                         tGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGccc                                       <  1:158901/251‑1 (MQ=255)
                                                                                                                                                                                                      attattGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCagcagc                          <  1:69320/251‑1 (MQ=255)
                                                                                                                                                                                                         attGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgag                       >  1:483325/1‑249 (MQ=255)
                                                                                                                                                                                                         attGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgag                       >  1:232337/1‑249 (MQ=255)
                                                                                                                                                                                                          ttGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGcagcag                           >  1:789812/1‑246 (MQ=255)
                                                                                                                                                                                                          ttGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgagg                      >  1:781539/1‑248 (MQ=255)
                                                                                                                                                                                                                 caccgcGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTgggg               <  1:759887/249‑1 (MQ=255)
                                                                                                                                                                                                                     gcgACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                               <  1:283615/183‑1 (MQ=255)
                                                                                                                                                                                                                            caaTTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt    <  1:143521/251‑1 (MQ=255)
                                                                                                                                                                                                                              attttattttGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGTCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGGATCGCCGTGc                                            >  1:615681/1‑209 (MQ=255)
                                                                                                                                                                                                                              attttattttGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCg  <  1:445071/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
TGACAAACTGGAGCGCTTTATTGATGAATTGTTGTTGTGTGATTGCGAAGATCGTCTGACGCAACCGCGTTACGATGCCAACCTTGCCAAACAGATTCAGATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCG  >  1/887247‑887719

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