Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 886,901 2 bp→AA coding (596‑597/1923 nt) 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,9010CA100.0% 123.6 / NA 43A199E (GCG→GAA) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/22);  total (21/22)
*1886,9020GA100.0% 123.8 / NA 43A199E (GCG→GAA) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/22);  total (21/22)

CTTTCTGCAAGAAAAAATCCGGAAATACAGCTTCCCCATTTATCAACATTTGCCGTCTGCGCTCGTTTTCCCTAAGATTGAAAAGTCGCTTTCAGAACAACACCAACCGACCTTGGAACAAGCCCATTTATTGCGGCAAATGGCGACATCTGACGCGGAGATGTTAATCGTTACGGCAAAGCGCGGGCGGGGAAAATCCGCTTTGGCGGGCTTGTTCGCCAAACAGCAATTAGCGGCAAATAAGCCCGTTCTATTGACGGCGCCGAACAAAAGTGCGGTCAAGATTTTCAATGAATTTTCCGGTGTCGAAATCCCCTTTATGCCACCTGATGAACTGCGCCAAAAGGTTGGTGAAAATCCACAGCAATTTGCGGATTCGTGGTTATTTGTGGACGAAGCGGCGATGAT  >  1/886668‑887075
                                                                                                                                                                                                                                         ||                                                                                                                                                                             
cTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTc                                                                                                                                                               <  1:204357/251‑1 (MQ=255)
cTTTCTGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTc                                                                                                                                                               <  1:515335/251‑1 (MQ=255)
     tGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          <  1:578806/251‑1 (MQ=255)
     tGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          <  1:65902/251‑1 (MQ=255)
     tGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          <  1:396196/251‑1 (MQ=255)
     tGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          <  1:471226/251‑1 (MQ=255)
     tGCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          <  1:660618/251‑1 (MQ=255)
      gCAAGAAAAAATCCAAAAATACGGCTTTCCCATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGa                                                                                                                                                         >  1:324372/1‑251 (MQ=255)
                                  cccATTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          >  1:782676/1‑222 (MQ=255)
                                     aTTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTAt                                                                                                                                                            <  1:89370/217‑1 (MQ=255)
                                     aTTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATAtttt                                                                                                                          <  1:128232/251‑1 (MQ=255)
                                     aTTTATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATAtttt                                                                                                                          <  1:246325/251‑1 (MQ=255)
                                      tttATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTc                                                                                                                         <  1:116385/251‑1 (MQ=255)
                                       ttATCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCa                                                                                                                        >  1:795066/1‑251 (MQ=255)
                                         aTCAATATCTGCCGCCCGCGCTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACgg                                                                                                                                                      >  1:629920/1‑219 (MQ=255)
                                                tctgccgcccgcgcTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAAtttg                                                                                                               >  1:791685/11‑250 (MQ=255)
                                                          gcgcTGGTCTTCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCgctgcgga                                                                                                     >  1:853880/1‑244 (MQ=255)
                                                                gtcttCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATcc                                                                                                <  1:641658/247‑1 (MQ=255)
                                                                    tCCCTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGttt                                                                                           <  1:891812/251‑1 (MQ=255)
                                                                     cccTAGGATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          >  1:471769/1‑187 (MQ=255)
                                                                          ggATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGAc                                                                                                                                                        <  1:613195/183‑1 (MQ=255)
                                                                          ggATTGAAACGTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATcc                                                                                                <  1:837591/239‑1 (MQ=255)
                                                                                    gTCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGaa                                                                            >  1:77978/1‑250 (MQ=255)
                                                                                     tCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAAc                                                                           <  1:860768/250‑1 (MQ=255)
                                                                                     tCGCTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAAc                                                                           <  1:177452/250‑1 (MQ=255)
                                                                                        cTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTgcg                                                                       >  1:579153/1‑251 (MQ=255)
                                                                                        cTTTCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGACAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          >  1:71132/1‑168 (MQ=255)
                                                                                          ttCAGAACTACATCAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTgcgt                                                                      >  1:832531/1‑249 (MQ=255)
                                                                                               aactacatcAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTgcgtcgaaa                                                                 <  1:536029/242‑7 (MQ=255)
                                                                                                   acatcAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTgcgtcgaaagattg                                                            >  1:529655/5‑240 (MQ=255)
                                                                                                    catcAACCGACCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTgg                                                           >  1:536541/4‑251 (MQ=255)
                                                                                                             aCCTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTg                                                                                                                                                          >  1:645381/1‑147 (MQ=255)
                                                                                                              ccTTGGAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATccg                                                 >  1:450363/1‑250 (MQ=255)
                                                                                                                   gAACAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTAtat                                                                                       >  1:418508/1‑206 (MQ=255)
                                                                                                                    aaCAAACCCATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCaa                                           <  1:14588/251‑1 (MQ=255)
                                                                                                                        aacccATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAg                                                                                                                                         <  1:724398/151‑1 (MQ=255)
                                                                                                                          cccATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCg                                     >  1:897422/1‑251 (MQ=255)
                                                                                                                          cccATTTATTGCGGCAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCg                                     >  1:272634/1‑251 (MQ=255)
                                                                                                                                        cAAATGGCGACATCCGGTGCAGAGATGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCGGATTTATGGTTAtt                       <  1:834540/251‑1 (MQ=255)
                                                                                                                                             ggcgacatccggtgcagagaTGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCGGATTTATGGTTATTTGTgg                  >  1:68449/17‑251 (MQ=255)
                                                                                                                                                    tccggtgcagagaTGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCGGATTTATGGTTATTTGTGGATGAAgc           <  1:278547/242‑1 (MQ=255)
                                                                                                                                                     ccggtgcagagaTGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCGGATTTATGGTTATTTGTGGATGAAgcg          >  1:443875/9‑251 (MQ=255)
                                                                                                                                                             gagaTGTTAATCGTTACGGCAAAGCGTGGGCGGGGCAAATCCGCTTTGGCAGGCTTGTTCGCCAAGCAGCAATTGGAAGCAAATCAGCCTGTTCTATTGACGGCGCTGAACAAAAGTGCGGTCAATATTTTCAATGAATTTGCCGCTGCGGAAATCCCGTTTATATCGCCTGATGAACTGCGTCGAAAGATTGGTGAAAATCCGCAGCAATTTGCGGATTTATGGTTATTTGTGGATGAAGCGGCgatgat  >  1:694667/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ||                                                                                                                                                                             
CTTTCTGCAAGAAAAAATCCGGAAATACAGCTTCCCCATTTATCAACATTTGCCGTCTGCGCTCGTTTTCCCTAAGATTGAAAAGTCGCTTTCAGAACAACACCAACCGACCTTGGAACAAGCCCATTTATTGCGGCAAATGGCGACATCTGACGCGGAGATGTTAATCGTTACGGCAAAGCGCGGGCGGGGAAAATCCGCTTTGGCGGGCTTGTTCGCCAAACAGCAATTAGCGGCAAATAAGCCCGTTCTATTGACGGCGCCGAACAAAAGTGCGGTCAAGATTTTCAATGAATTTTCCGGTGTCGAAATCCCCTTTATGCCACCTGATGAACTGCGCCAAAAGGTTGGTGAAAATCCACAGCAATTTGCGGATTCGTGGTTATTTGTGGACGAAGCGGCGATGAT  >  1/886668‑887075

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