Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,474,214 G→A L622L (CTG→TTG)  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,474,2140GA97.0% 170.7 / NA 66L622L (CTG→TTG) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (33/31);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGGGATCTTCATCAAACTTCACCACCGGTACGCAGCCTCGGCAAGTGCCTTGATCCACGCCTGCTTCTTTGCAGAAATCGGAATACGCATCCCAAGCGCGTTTTAACCAACCTTTTTTGGCTTCTTCAGCCATAGTTGCCTCCTTAGTCG‑TGCCGTTTTATTCTGGGTTGTTATCATAAGAAGGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTTTCAGCCACGTTTTTTTATGTTGAATGTCATCAGGGCGTTTGCAATAGGCTGTCAATAAAGGAAAATCCCGCTCCCCGCTTAATTTCATCAAGCGCCGAATAGCGGGCAAGGCAAACGCCAACGGCAGGAAGAAATCGGCAAAATCCGCCAAGCATTGCCAATCGGTTTCATTGAGCGACCAATCCGTTTCGGTTCCGACAAATTGCGTCGCCAGCGGGTGCAACGAAAGGGGC  >  1/1473979‑1474443
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     
gttcagGATCTTCGTCGAACTTCACCACCGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTAt                                                                                                                                                                                                                         >  1:653584/1‑251 (MQ=255)
     ggATCTTCGTCGAACTTCACCACCGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGa                                                                                                                                                                                                                    >  1:830768/1‑251 (MQ=255)
              tcgaaCTTCACCACCGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGtcatc                                                                                                                                                                                                            <  1:792636/247‑1 (MQ=255)
              tcgaaCTTCACCACCGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGtcatc                                                                                                                                                                                                            <  1:589438/247‑1 (MQ=255)
                        ccaccGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGc                                                                                                                                                                                                 <  1:537030/251‑1 (MQ=255)
                           ccGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAAt                                                                                                                                                                                              <  1:837236/251‑1 (MQ=255)
                           ccGGCACGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAAt                                                                                                                                                                                              <  1:232882/251‑1 (MQ=255)
                              gcaCGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAgg                                                                                                                                                                                           >  1:580848/3‑251 (MQ=255)
                              gcaCGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAg                                                                                                                                                                                            <  1:765256/248‑1 (MQ=255)
                              gcaCGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAg                                                                                                                                                                                            <  1:746828/248‑1 (MQ=255)
                               caCGCAGCCCCGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAgg                                                                                                                                                                                           <  1:728463/249‑1 (MQ=255)
                                     gccccGGCAAGTGCCTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATg                                                                                                                                                                                                                 <  1:222198/218‑1 (MQ=255)
                                                  ccTTGATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGa                                                                                                                                                                             <  1:651044/245‑1 (MQ=255)
                                                      gATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAggac                                                                                                                                                                                                       >  1:385174/1‑213 (MQ=255)
                                                      gATCCACACCGGCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATcctgtt                                                                                                                                                                   >  1:114500/1‑247 (MQ=255)
                                                           acaccggCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCg                                                                                                                                                              >  1:697933/7‑251 (MQ=255)
                                                           acaccggCTTCTTTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCg                                                                                                                                                              >  1:280371/7‑251 (MQ=255)
                                                                  cttcttTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTactt                                                                                                                                                       <  1:766491/251‑4 (MQ=255)
                                                                    tcttTGCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACttt                                                                                                                                                      >  1:674475/1‑250 (MQ=255)
                                                                         gCAGAAATCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAAt                                                                                                                                                                                   >  1:674247/1‑216 (MQ=255)
                                                                                tCGGAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGa                                                                                                                                         >  1:730850/1‑251 (MQ=255)
                                                                                  ggAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAAt                                                                                                                                       >  1:25708/1‑251 (MQ=255)
                                                                                   gAATAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAg                                                                                                                                                                               >  1:155664/1‑210 (MQ=255)
                                                                                    aaTAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAg                                                                                                                                     <  1:243985/251‑1 (MQ=255)
                                                                                    aaTAGGCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAg                                                                                                                                     <  1:474711/251‑1 (MQ=255)
                                                                                      taggCATCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACttt                                                                                                                                                      <  1:115330/229‑1 (MQ=255)
                                                                                        ggctgcccAAGCGCGTTTTAGCCAACCTTTTTTAGCTTCTTCAGTCGTATTTGCCTCCTTAGTCGATGCCGTTTTATGTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAGCCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGc                                                                                                                                                                                          <  1:700849/189‑1 (MQ=255)
                                                                                            tCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCCTCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGggcaag                                                                                                                             >  1:435084/1‑251 (MQ=255)
                                                                                            tCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAgg                                                                                                                                                                              <  1:6763/202‑1 (MQ=255)
                                                                                            tCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGggcaag                                                                                                                             >  1:198214/1‑251 (MQ=255)
                                                                                            tCCCAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGggcaag                                                                                                                             >  1:394941/1‑251 (MQ=255)
                                                                                              ccAAGCGCGTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGggcaagac                                                                                                                           <  1:490986/251‑3 (MQ=255)
                                                                                                   cgcgTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACg                                                                                                                      <  1:716448/251‑1 (MQ=255)
                                                                                                     cgTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGc                                                                                                                     >  1:497927/1‑250 (MQ=255)
                                                                                                     cgTTTTAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTGTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATcctgttc                                                                                                                                                                  <  1:516088/205‑6 (MQ=255)
                                                                                                         ttAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCgg                                                                                                                                  >  1:543942/1‑233 (MQ=255)
                                                                                                          tAACCAACCTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACgg                                                                                                               <  1:848823/251‑1 (MQ=255)
                                                                                                                accTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGgaa                                                                                                         >  1:50446/1‑251 (MQ=255)
                                                                                                                accTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTACTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAAt                                                                                                                                       >  1:118271/1‑221 (MQ=255)
                                                                                                                 ccTTTTTTAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGgaag                                                                                                        <  1:164548/251‑1 (MQ=255)
                                                                                                                    tttttAGCTTCTTCAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGaaa                                                                                                     <  1:599641/251‑1 (MQ=255)
                                                                                                                              cttcAGTCATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGgaaga                                                                                                       >  1:809826/1‑239 (MQ=255)
                                                                                                                                   gtcagatttGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGc                                                                                       >  1:82166/8‑250 (MQ=255)
                                                                                                                                   gtcATATTTGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGc                                                                                       >  1:272769/3‑250 (MQ=255)
                                                                                                                                      atatttGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGc                                                                                                                                                             <  1:872792/173‑1 (MQ=255)
                                                                                                                                        atttGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCATTAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCaaacc                                                                                 <  1:122563/249‑4 (MQ=255)
                                                                                                                                        atttGCCTCCTTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCaaac                                                                                  >  1:261971/3‑248 (MQ=255)
                                                                                                                                                cctTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATcctgttc                                                                                                                                                                  >  1:242297/1‑157 (MQ=255)
                                                                                                                                                cctTAGTCGATGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGccggt                                                                         >  1:26268/1‑248 (MQ=255)
                                                                                                                                                    agtcgatGCCGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGc                                                                                                                     >  1:660204/7‑203 (MQ=255)
                                                                                                                                                            ccGTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTg                                                              <  1:516656/250‑1 (MQ=255)
                                                                                                                                                              gTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTAt                                                                                                                                                                                                                         <  1:801166/93‑1 (MQ=255)
                                                                                                                                                              gTTTTATTTTGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAg                                                            >  1:152514/1‑250 (MQ=255)
                                                                                                                                                                   attttGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAAt                                                                                                                                       <  1:87641/166‑1 (MQ=255)
                                                                                                                                                                   attttGGGTTGTTATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCa                                                      >  1:332150/5‑251 (MQ=255)
                                                                                                                                                                             gttATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACggt                                                                                                              <  1:270628/185‑2 (MQ=255)
                                                                                                                                                                               tATCATAAGAAAGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGc                                                                                                                                                                                          <  1:95856/107‑1 (MQ=255)
                                                                                                                                                                                                            cccTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCa             >  1:785990/1‑251 (MQ=255)
                                                                                                                                                                                                            cccTTTCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCa             >  1:679235/1‑251 (MQ=255)
                                                                                                                                                                                                               tttCTTCACTTCCAACCGACAATTTGCCAGCCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCaatg          <  1:701517/251‑3 (MQ=255)
                                                                                                                                                                                                               tttCTTCACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCaatg          >  1:848411/1‑249 (MQ=255)
                                                                                                                                                                                                                  cttcACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGc              >  1:770914/1‑244 (MQ=255)
                                                                                                                                                                                                                  cttcACTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGaaa       <  1:567562/251‑1 (MQ=255)
                                                                                                                                                                                                                      aCTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAgg     <  1:242533/249‑1 (MQ=255)
                                                                                                                                                                                                                       cTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAgggg   >  1:270092/1‑250 (MQ=255)
                                                                                                                                                                                                                       cTTCCAACCGACAATTTGCCAACCACGTTTTTTTATTTTGAATGTCATCAGGACGCTTGCAATAGGCTGCCAATAAAGGAAAATCCTGTTCCCCGCTTACTTTCATCAAGCGCCGAATAGCGGGCAAGACAAACGCCAACGGTAGGAAGAAATCGGCAAAATCCGCCAAACCTTGCCGGTCGGCTTCATTGAGCGACCAATCCGTTTCGCTTCCGGCAAATTGTGTCGCCAGCGGGTGCAATGAAAGGGGc  <  1:779751/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                     
GTTCGGGATCTTCATCAAACTTCACCACCGGTACGCAGCCTCGGCAAGTGCCTTGATCCACGCCTGCTTCTTTGCAGAAATCGGAATACGCATCCCAAGCGCGTTTTAACCAACCTTTTTTGGCTTCTTCAGCCATAGTTGCCTCCTTAGTCG‑TGCCGTTTTATTCTGGGTTGTTATCATAAGAAGGCGAAATTTTTTTGCAACCCTTTCTTCACTTCCAACCGACAATTTTTCAGCCACGTTTTTTTATGTTGAATGTCATCAGGGCGTTTGCAATAGGCTGTCAATAAAGGAAAATCCCGCTCCCCGCTTAATTTCATCAAGCGCCGAATAGCGGGCAAGGCAAACGCCAACGGCAGGAAGAAATCGGCAAAATCCGCCAAGCATTGCCAATCGGTTTCATTGAGCGACCAATCCGTTTCGGTTCCGACAAATTGCGTCGCCAGCGGGTGCAACGAAAGGGGC  >  1/1473979‑1474443

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