Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,475,617 T→A Q154L (CAA→CTA)  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
*11,475,6170TA100.0% 196.8 / NA 74Q154L (CAA→CTA) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (39/35);  total (39/35)

TTTTGGCGCAGTTCATCAGGTGGCATAAAGGGGATTTCGACACCGGAAAATTCATTGAAAATCTTGACCGCACTTTTGTTCGGCGCCGTCAATAGAACGGGCTTATTTGCCGCTAATTGCTGTTTGGCGAACAAGCCCGCCAAAGCGGATTTTCCCCGCCCGCGCTTTGCCGTAACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGACTTTTCAATCTTAGGGAAAACGAGCGCAGACGGCAAATGTTGATAAATGGGGAAGCTGTATTTCCGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACTGCTCCAACGCAGGCTGTCGCCGTCGGGTTGTGCCGCCAAATCCTGCCACTCATTGAATAACAGCAACAATGTGCCGCCTTCCGACAACGCACCGGCGGCAATAGC  >  1/1475371‑1475858
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
tttcgacgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATgtagtt                                                                                                                                                                                                                                                >  1:536029/7‑250 (MQ=255)
     acgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGaa                                                                                                                                                                                                                                           <  1:832531/249‑1 (MQ=255)
      cgcAGTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAg                                                                                                                                                                                                                                         <  1:579153/251‑1 (MQ=255)
          gTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGa                                                                                                                                                                                                                                      >  1:860768/1‑250 (MQ=255)
          gTTCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGa                                                                                                                                                                                                                                      >  1:177452/1‑250 (MQ=255)
           ttCATCAGGCGATATAAACGGGATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGAc                                                                                                                                                                                                                                     <  1:77978/250‑1 (MQ=255)
                          aaacgggATTTCCGCAGCGGCAAATTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGa                                                                                                                                                                                                                     >  1:891812/5‑251 (MQ=255)
                               ggatttccgcagcggcaaaTTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATcc                                                                                                                                                                                                                           >  1:837591/17‑239 (MQ=255)
                               ggatttccgcagcggcaaaTTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGaagac                                                                                                                                                                                                                 >  1:641658/17‑247 (MQ=255)
                                    tccgcagcggcaaaTTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAgcgc                                                                                                                                                                                                           <  1:853880/240‑1 (MQ=255)
                                              caaaTTCATTGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCaga                                                                                                                                                                                                 <  1:791685/250‑3 (MQ=255)
                                                       tGAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATaa                                                                                                                                                                                        <  1:795066/251‑1 (MQ=255)
                                                        gAAAATATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATaaa                                                                                                                                                                                       >  1:116385/1‑251 (MQ=255)
                                                         aaaaTATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAAt                                                                                                                                                                                      >  1:128232/1‑251 (MQ=255)
                                                         aaaaTATTGACCGCACTTTTGTTCAGCGCCGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAAt                                                                                                                                                                                      >  1:246325/1‑251 (MQ=255)
                                                                                     ccGTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGAt                                                                                                                                                                                          <  1:629920/219‑1 (MQ=255)
                                                                                       gTCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATcc                                                                                                                                                                                                                           >  1:613195/1‑183 (MQ=255)
                                                                                        tCAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGc                                                                                                                                                       <  1:324372/251‑1 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGTCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAg                                                                                                                                                                                                                                         <  1:71132/168‑1 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAggg                                                                                                                                                                                                                      <  1:471769/187‑1 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATggg                                                                                                                                                                                   <  1:782676/222‑1 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCa                                                                                                                                                      >  1:65902/1‑251 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCa                                                                                                                                                      >  1:660618/1‑251 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCa                                                                                                                                                      >  1:578806/1‑251 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCa                                                                                                                                                      >  1:396196/1‑251 (MQ=255)
                                                                                         cAATAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCa                                                                                                                                                      >  1:471226/1‑251 (MQ=255)
                                                                                           aTAGAACAGGCTGATTTGCTTCCAATTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAAt                                                                                                                                                                                      >  1:89370/1‑217 (MQ=255)
                                                                                              gaacaggctgatttgcttccaaTTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAg                                                                                                                                                 >  1:204357/21‑251 (MQ=255)
                                                                                              gaacaggctgatttgcttccaaTTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAg                                                                                                                                                 >  1:515335/21‑251 (MQ=255)
                                                                                                   ggctgatttgcttccaaTTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAggg                                                                                                                                                                                                                      >  1:48300/16‑177 (MQ=255)
                                                                                                            gcttccaaTTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATcc                                                                                                                                                                                                                           >  1:335158/7‑162 (MQ=255)
                                                                                                                   aTTGCTGCTTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGttt                                                                                                                            <  1:738563/251‑1 (MQ=255)
                                                                                                                      gctgctTGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTAt                                                                                                                          >  1:669075/1‑250 (MQ=255)
                                                                                                                       ctgcttGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTAAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATc                                                                                                                         <  1:463983/246‑1 (MQ=255)
                                                                                                                          cttGGCGAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCggg                                                                                                                      <  1:518565/249‑1 (MQ=255)
                                                                                                                                gAACAAGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGttttt                                                                                                               >  1:104312/1‑251 (MQ=255)
                                                                                                                                     aGCCTGCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCACTCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCAccg                                                                                                          >  1:718452/1‑249 (MQ=255)
                                                                                                                                          gCCAAAGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGAGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCa                                                                                                     >  1:582660/1‑251 (MQ=255)
                                                                                                                                              aaGCGGATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACg                                                                                                  >  1:478547/1‑250 (MQ=255)
                                                                                                                                                   gATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCaagc                                                                                             <  1:639360/250‑4 (MQ=255)
                                                                                                                                                   gATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCaag                                                                                              >  1:296060/1‑247 (MQ=255)
                                                                                                                                                   gATTTGCCCCGCCCACGCTTTGCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCaag                                                                                              >  1:241184/1‑247 (MQ=255)
                                                                                                                                                                        gCCGTAACGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATgg                                                                                                                                                                                    >  1:31774/1‑142 (MQ=255)
                                                                                                                                                                              aCGATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATc                                                                 <  1:780503/251‑1 (MQ=255)
                                                                                                                                                                                gATTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATc                                                                 >  1:707442/1‑249 (MQ=255)
                                                                                                                                                                                 aTTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGAt                                                                                                                                                                                          <  1:295264/127‑1 (MQ=255)
                                                                                                                                                                                 aTTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGc                                                                                                 >  1:171207/1‑216 (MQ=255)
                                                                                                                                                                                 aTTAACATCTCTGCACCGGATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTcaccatcaggt                                                                               <  1:539034/234‑11 (MQ=255)
                                                                                                                                                                                        tctctgcaccggATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGc                                                                                                 >  1:253815/12‑209 (MQ=255)
                                                                                                                                                                                                 cggATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAAc                                              <  1:775271/249‑1 (MQ=255)
                                                                                                                                                                                                  ggATGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGAt                                                                                                                                                                                          >  1:143867/2‑110 (MQ=255)
                                                                                                                                                                                                    aTGTCGCCATTTGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAg                                            <  1:661606/250‑1 (MQ=255)
                                                                                                                                                                                                               tGCCGCAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATgcc                                <  1:195416/251‑1 (MQ=255)
                                                                                                                                                                                                                    cAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATc                                                                                                                         >  1:26160/1‑157 (MQ=255)
                                                                                                                                                                                                                    cAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTc                                                                                                              <  1:101027/168‑1 (MQ=255)
                                                                                                                                                                                                                    cAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTcacc                                                                                      <  1:187747/192‑4 (MQ=255)
                                                                                                                                                                                                                    cAATAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCtt                           >  1:581539/1‑251 (MQ=255)
                                                                                                                                                                                                                      aTAAATGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTc                          >  1:168614/1‑250 (MQ=255)
                                                                                                                                                                                                                        aaaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGcaacaa                                      >  1:181389/1‑236 (MQ=255)
                                                                                                                                                                                                                        aaaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATgc                                 >  1:755147/1‑241 (MQ=255)
                                                                                                                                                                                                                        aaaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGa                       >  1:58554/1‑251 (MQ=255)
                                                                                                                                                                                                                         aaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                      <  1:408960/251‑1 (MQ=255)
                                                                                                                                                                                                                         aaTGGGTTTGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                      <  1:380027/251‑1 (MQ=255)
                                                                                                                                                                                                                             ggtttGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTAt                                                                                                                          <  1:520638/144‑1 (MQ=255)
                                                                                                                                                                                                                             ggtttGTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACagcg                  >  1:200952/4‑248 (MQ=255)
                                                                                                                                                                                                                                  gTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGc                                                                                                                                     <  1:693442/131‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCAc                                                                                                            <  1:830534/156‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAAc                                              <  1:49114/218‑1 (MQ=255)
                                                                                                                                                                                                                                  gTTCCAAGGTCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACcg             >  1:840272/1‑251 (MQ=255)
                                                                                                                                                                                                                                           tCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAAtg    <  1:49567/251‑2 (MQ=255)
                                                                                                                                                                                                                                           tCGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAAt     <  1:578240/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTTTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCAAtgg   <  1:739493/251‑3 (MQ=255)
                                                                                                                                                                                                                                              gTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGAc                      <  1:120324/230‑1 (MQ=255)
                                                                                                                                                                                                                                              gTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACcggcg          <  1:182502/242‑1 (MQ=255)
                                                                                                                                                                                                                                              gTTGATGTAGTTCTGAAAGCGACGTTTCAATCCTAGGGAAGACCAGCGCGGGCGGCAGATATTGATAAATGGGAAAGCCGTATTTTTGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACCGCTCCAACGCAAGCTGTCACCATCAGGTTGCACCGCCAAATCCTGCCAATCATTGAATAACAGCAACAATATGCCGCCTTCCGACAGCGCACCGGCGGCatattg  <  1:650412/250‑6 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
TTTTGGCGCAGTTCATCAGGTGGCATAAAGGGGATTTCGACACCGGAAAATTCATTGAAAATCTTGACCGCACTTTTGTTCGGCGCCGTCAATAGAACGGGCTTATTTGCCGCTAATTGCTGTTTGGCGAACAAGCCCGCCAAAGCGGATTTTCCCCGCCCGCGCTTTGCCGTAACGATTAACATCTCCGCGTCAGATGTCGCCATTTGCCGCAATAAATGGGCTTGTTCCAAGGTCGGTTGGTGTTGTTCTGAAAGCGACTTTTCAATCTTAGGGAAAACGAGCGCAGACGGCAAATGTTGATAAATGGGGAAGCTGTATTTCCGGATTTTTTCTTGCAGAAAGGCAATAAAGTGCGGTGTGTTTATCGGGTGTTTTTCACTGCTCCAACGCAGGCTGTCGCCGTCGGGTTGTGCCGCCAAATCCTGCCACTCATTGAATAACAGCAACAATGTGCCGCCTTCCGACAACGCACCGGCGGCAATAGC  >  1/1475371‑1475858

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