Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,027,756 A→G K98R (AAA→AGA)  unknown → Pseudouridine synthase (EC 4.2.1.70), VC1668 type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,027,7560AG100.0% 204.0 / NA 74K98R (AAA→AGA) unknownPseudouridine synthase (EC 4.2.1.70), VC1668 type
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/42);  total (32/42)

ACAAACCGAACGGTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTTCCATGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAAAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAGCGCCGATCGCACCTATCTGCACGCTTATCA  >  1/1027510‑1028002
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aCAAACCGAACGGTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATc                                                                                                                                                                                                                                                    <  1:744313/251‑1 (MQ=255)
         aCGGTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCA‑CGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGata                                                                                                                                                                                                                                          <  1:162711/251‑1 (MQ=255)
            gTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAAt                                                                                                                                                                                                                                                     <  1:183970/238‑1 (MQ=255)
            gTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATg                                                                                                                                                                                                                                        >  1:467919/1‑251 (MQ=255)
            gTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATg                                                                                                                                                                                                                                        >  1:787115/1‑251 (MQ=255)
            gTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATg                                                                                                                                                                                                                                        >  1:391302/1‑251 (MQ=255)
            gTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGGGATATg                                                                                                                                                                                                                                        >  1:598606/1‑251 (MQ=255)
                cAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCaaa                                                                                                                                                                                                                                    >  1:832897/1‑251 (MQ=255)
                       cacaAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAgcgcg                                                                                                                                                                                                                             >  1:272409/1‑251 (MQ=255)
                       cacaAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAgcgcg                                                                                                                                                                                                                             >  1:666593/1‑251 (MQ=255)
                                     cGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGa                                                                                                                                                                                                               <  1:193893/251‑1 (MQ=255)
                                        cATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAAc                                                                                                                                                                                                            <  1:654692/251‑1 (MQ=255)
                                         aTCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAAcg                                                                                                                                                                                                           <  1:170240/251‑2 (MQ=255)
                                           cGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGa                                                                                                                                                                                                         <  1:622489/251‑1 (MQ=255)
                                            ggCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCaa                                                                                                                                                                                                                                     >  1:86722/1‑222 (MQ=255)
                                            ggCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGa                                                                                                                                                                                                         >  1:389970/1‑250 (MQ=255)
                                            ggCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGAGCTTGGAAACTGa                                                                                                                                                                                                         >  1:901806/1‑250 (MQ=255)
                                              cTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGAcc                                                                                                                                                                                                       <  1:701688/250‑1 (MQ=255)
                                                 aaCCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCcaaa                                                                                                                                                                                                   >  1:611783/1‑251 (MQ=255)
                                                   cccGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAAGGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCcaaagc                                                                                                                                                                                                 <  1:434232/251‑1 (MQ=255)
                                                                   aCAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATACAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCg                                                                                                                                                                                 >  1:264110/1‑251 (MQ=255)
                                                                                ccGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAAc                                                                                                                                                                                                            <  1:377387/211‑1 (MQ=255)
                                                                                 cGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTtgtg                                                                                                                                                                   <  1:345946/251‑1 (MQ=255)
                                                                                   cAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTgtgt                                                                                                                                                                  <  1:188303/250‑1 (MQ=255)
                                                                                       tttGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCAt                                                                                                                                                             <  1:463963/251‑1 (MQ=255)
                                                                                              ggTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCg                                                                                                                                                      >  1:886653/1‑251 (MQ=255)
                                                                                                 tCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCAc                                                                                                                                                   <  1:397366/251‑1 (MQ=255)
                                                                                                          ggATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTg                                                                                                                                           >  1:610557/1‑250 (MQ=255)
                                                                                                            aTAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTgcgc                                                                                                                                        <  1:70784/251‑1 (MQ=255)
                                                                                                              aaaGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCAc                                                                                                                                                   >  1:538789/1‑238 (MQ=255)
                                                                                                                      cTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTc                                                                                                                              <  1:466888/251‑1 (MQ=255)
                                                                                                                        cGGGGTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTg                                                                                                                            <  1:718292/251‑1 (MQ=255)
                                                                                                                          gggTTATTAATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAAtt                                                                                                                                             >  1:720596/1‑232 (MQ=255)
                                                                                                                             ttattaATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGc                                                                                                                                                                                                                                       <  1:518461/139‑1 (MQ=255)
                                                                                                                             ttattaATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAgg                                                                                                                                                                                                                       <  1:20452/155‑1 (MQ=255)
                                                                                                                             ttattaATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAAc                                                                                                                        >  1:832837/1‑250 (MQ=255)
                                                                                                                               attaATGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCa                                                                                                                                                                                                                                                   >  1:540357/1‑125 (MQ=255)
                                                                                                                                    tGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAAc                                                                                                                 >  1:803892/1‑250 (MQ=255)
                                                                                                                                    tGTTAGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAAc                                                                                                                 >  1:330863/1‑250 (MQ=255)
                                                                                                                                        aGCATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAAc                                                                                                                                                                                                            >  1:30242/1‑155 (MQ=255)
                                                                                                                                          cATTAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACaaa                                                                                                          >  1:552416/1‑251 (MQ=255)
                                                                                                                                             tAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGata                                                                                                                                                                                                                                          <  1:129766/120‑1 (MQ=255)
                                                                                                                                             tAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCg                                                                                                                      >  1:198196/1‑236 (MQ=255)
                                                                                                                                             tAAATGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAAc                                                                                                        >  1:80352/1‑250 (MQ=255)
                                                                                                                                              aaaTGAACAGGCTGCCGTGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAAt                                                                                                                                                                                                                                                     <  1:594757/108‑1 (MQ=255)
                                                                                                                                                        gctgccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGc                                                                                                <  1:445220/240‑1 (MQ=255)
                                                                                                                                                        gctgccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTc                                                                                              <  1:618997/242‑1 (MQ=255)
                                                                                                                                                           gccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCttattt                                                                                                                                  <  1:623956/206‑1 (MQ=255)
                                                                                                                                                           gccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAACGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCAt                                                                                                                                                             <  1:558887/179‑1 (MQ=255)
                                                                                                                                                            ccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCgtgt                                                                                         <  1:785564/247‑1 (MQ=255)
                                                                                                                                                            ccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCgtgt                                                                                         <  1:412684/247‑1 (MQ=255)
                                                                                                                                                            ccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCgtgt                                                                                         <  1:5982/247‑1 (MQ=255)
                                                                                                                                                            ccgtGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCgtgt                                                                                         >  1:152688/4‑250 (MQ=255)
                                                                                                                                                                      gTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAAc                                                                                                                        >  1:213212/1‑209 (MQ=255)
                                                                                                                                                                      gTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGaaaa                                                                              >  1:317412/1‑251 (MQ=255)
                                                                                                                                                                                tttGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGc                                                                      <  1:880926/249‑1 (MQ=255)
                                                                                                                                                                                tttGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCACATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGc                                                                      <  1:215933/249‑1 (MQ=255)
                                                                                                                                                                                    cATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTcc                                                                 <  1:878745/250‑1 (MQ=255)
                                                                                                                                                                                    cATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTcc                                                                 <  1:277298/250‑1 (MQ=255)
                                                                                                                                                                                                 tCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATg                                                   >  1:891304/1‑251 (MQ=255)
                                                                                                                                                                                                       aaCTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGt                                             >  1:690762/1‑251 (MQ=255)
                                                                                                                                                                                                        aCTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGc                                                                      <  1:396341/225‑1 (MQ=255)
                                                                                                                                                                                                         cTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGTTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTAt                                           <  1:447769/251‑1 (MQ=255)
                                                                                                                                                                                                         cTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTAt                                           <  1:660279/251‑1 (MQ=255)
                                                                                                                                                                                                         cTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTAt                                           <  1:653281/251‑1 (MQ=255)
                                                                                                                                                                                                         cTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTAt                                           <  1:63028/251‑1 (MQ=255)
                                                                                                                                                                                                         cTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTAt                                           <  1:516295/251‑1 (MQ=255)
                                                                                                                                                                                                                  cACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGAt                                  <  1:525337/251‑1 (MQ=255)
                                                                                                                                                                                                                        cAGCTAACAAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAgcgc                             <  1:219825/250‑1 (MQ=255)
                                                                                                                                                                                                                             cacaAACCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAGCGCCGATCg                       >  1:561345/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaCCAAAGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACGAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAGCGCCGATCGCAcc                   <  1:246205/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCgtgt                                                                                         >  1:423548/1‑176 (MQ=255)
                                                                                                                                                                                                                                      aaGAAAAAGCAGGGGAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAGCGCCGATCGCACCTATCt              <  1:265772/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gggAGAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAGCGCCGATCGCACCTATCTGCACGCTTATCa  >  1:103409/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
ACAAACCGAACGGTATCAGCGTACACAAAGACGATGCCGACATCGGCTTAACCCGACTGGTGGCGCAACAGCTTGGCGTGCCGCAGGTTTGGCTGGTTCATCGGCTGGATAAAGTGACCTCGGGGTTATTAATGTTAGCATTAAATGAACAGGCTTCCATGACGTTGTCACAAAAGTTTGCATATCATCAAATTCAGAAAACTTACCTCGCACTGGCAGCTCACAAACCAAAGAAAAAGCAGGGGAAAATCAGCGGTGATATGCAAAAAGCGCGCCGAGGTGCTTGGAAACTGACCCAAAGCAAAGAAAATCCTGCCGTGACAGATTTTGTGTCCCATAGTCTCGCACCAAATTTGCGCTTATTTATTCTGCAACCGAAAACCGACAAAACCCATCAGCTCCGTGTCGCTATGAAAAGTTTAGGCAGTCCGATTTTAGGCGATGCTTTGTATGGCGGCGATAGCGCCGATCGCACCTATCTGCACGCTTATCA  >  1/1027510‑1028002

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