Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,474,818 C→T R420R (CGG→CGA)  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,8180CT98.6% 186.2 / NA 70R420R (CGG→CGA) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (36/33);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTCGGGGGTGTAACCGAAACTGACGGAAAGGAAATCCACCGCACGTTGTTGTTTAACTC‑GTTCGATTAAACGTTGCCCCAGTCCCCGCTGCTGCCAGTTGGGTTGCACGGCGATACGGGAAATGCGCAGTGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTACCACGTGGACGACGTTGTCCGCGCCGAATTTTACGAATTAATGCGTCGTCACCAAATTTGCCCTCAGGCACCGCCCACACACAACCCAACAAGGCCTCCTGCGTCTGCGCAAGATAAAATTGTTGTCGCGGCGCATCCAATAATCGACGTAAATCCAACGGCGAGGTTCGGTAATGCGCCAAGGTCAGCAAGCCATAAACCGATTCGATTTGATAAGGCGGAATCTGTTCTTCGTAAGAAAT  >  1/1474599‑1475036
                                                                                                                                                                                                                            |                                                                                                                                                                                                                          
cagttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaa                                                                                                                                                                                             >  1:305880/5‑240 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAACCACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:446276/3‑238 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAAtttgac                                                                                                                                                                                                                  >  1:143746/3‑225 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:543915/3‑238 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:578174/3‑238 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:149617/3‑238 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:836462/3‑238 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:217141/3‑238 (MQ=255)
  gttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcat                                                                                                                                                                                                      <  1:366912/238‑3 (MQ=255)
   ttCGGGGGTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAACCGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaatc                                                                                                                                                                                          <  1:770770/250‑15 (MQ=255)
       ggggTGTAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcat                                                                                                                                                                                                      >  1:848980/1‑233 (MQ=255)
           tgtAACCGAAACTGACGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaatcctccttc                                                                                                                                                                                   <  1:338311/250‑22 (MQ=255)
                         aCGGAAAGAAAGTCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCGGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATGCGGGAAATGCGCAAAGATCGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcaca                                                                                                                                                                     <  1:85052/250‑1 (MQ=255)
                         aCGGAAAGAAAATCCACCGCAGGTTGTTGTTTAA‑TCGGTTCGATTAAACGTGGCCCCAGCCCTCGCTGCGGCCAGTTGGGCTGCACGGCGATACGGAAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacac                                                                                                                                                                    <  1:16077/251‑1 (MQ=255)
                         aCGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacac                                                                                                                                                                    <  1:598662/251‑1 (MQ=255)
                          cGGAAAGAAAATCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatca                                                                                                                                                                                                    >  1:390200/1‑214 (MQ=255)
                            gaaagaaaaTCCCCCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacacac                                                                                                                                                                  <  1:30271/250‑1 (MQ=255)
                            gaaagaaaaTCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacacac                                                                                                                                                                  <  1:892622/250‑1 (MQ=255)
                                   aaTCCACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaac                                                                                                                                                          <  1:386119/251‑1 (MQ=255)
                                     cccACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacag                                                                                                                                                        <  1:596187/250‑2 (MQ=255)
                                     cccACCGCACGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacag                                                                                                                                                        <  1:308059/250‑2 (MQ=255)
                                             aCGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacacac                                                                                                                                                                  <  1:297852/233‑1 (MQ=255)
                                             aCGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacagggtttcc                                                                                                                                                 >  1:91252/1‑242 (MQ=255)
                                              cGTTGTTGTTTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacac                                                                                                                                                                    <  1:323252/230‑1 (MQ=255)
                                                ttgttgttTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacacac                                                                                                                                                                  >  1:471552/1‑230 (MQ=255)
                                                ttgttgttTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacac                                                                                                                                                                    >  1:437368/1‑228 (MQ=255)
                                                ttgttgttTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacagggtttcccgcg                                                                                                                                             >  1:585054/1‑239 (MQ=255)
                                                  gttgttTAA‑TCTGTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacagggtttcccgcgtt                                                                                                                                           >  1:144890/1‑237 (MQ=255)
                                                     tgtttaatctgTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcat                                                                                                                                                                                                      <  1:325522/180‑3 (MQ=255)
                                                      gtttaatctgTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcc                                                                                                                                                                                                 <  1:487161/185‑8 (MQ=255)
                                                      gtttaatctgTTCGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacacac                                                                                                                                                                  >  1:539036/10‑225 (MQ=255)
                                                                  cGATTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCaaaataaaat                                                                                                                            >  1:445071/1‑251 (MQ=255)
                                                                    aTTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTTCACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacagggtttcccgcgtttgtgcaaaat                                                                                                                                 >  1:9417/1‑220 (MQ=255)
                                                                    aTTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacagggtttcccgcg                                                                                                                                             <  1:790218/232‑13 (MQ=255)
                                                                    aTTAAACGTTGCCCCAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAAttg                                                                                                                          >  1:143521/1‑251 (MQ=255)
                                                                               ccccAGCCCTCGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCgcggtg                                                                                                               >  1:759887/1‑249 (MQ=255)
                                                                                    gccctcGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatca                                                                                                                                                                                                    >  1:227295/6‑157 (MQ=255)
                                                                                      cctcGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCaa                                                                                                        <  1:781539/248‑1 (MQ=255)
                                                                                       ctcGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCaat                                                                                                       <  1:483325/249‑1 (MQ=255)
                                                                                       ctcGCTGCTGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCaat                                                                                                       <  1:232337/249‑1 (MQ=255)
                                                                                          gctgctGCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCaataat                                                                                                    >  1:69320/1‑251 (MQ=255)
                                                                                           ctgctgCCAGTTGGGCTGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCaa                                                                                                        <  1:789812/246‑1 (MQ=255)
                                                                                                       gggctGCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCa                                                                                       >  1:158901/5‑251 (MQ=255)
                                                                                                            gCACGGCGATCCGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGACACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCaaaataaaat                                                                                                                            <  1:615681/209‑1 (MQ=255)
                                                                                                            gCACGGCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcatcatcctcaaat                                                                                                                                                                                           >  1:147719/1‑133 (MQ=255)
                                                                                                                 gCGATACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcac                                                                                                                                                                      <  1:830742/162‑1 (MQ=255)
                                                                                                                     tACGGGAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCcacacac                                                                                                                                                                  <  1:833883/162‑1 (MQ=255)
                                                                                                                         ggAAATGCGCAAAGGACCCAATTTCCACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCTCCGCCCACACTCAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCTCGGTGCATCCAATAATCGTCGTAAATCCAACGGCGAGGtgctgtaa                                                                      <  1:210964/250‑8 (MQ=255)
                                                                                                                         ggAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      <  1:698893/250‑1 (MQ=255)
                                                                                                                         ggAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      <  1:51847/250‑1 (MQ=255)
                                                                                                                         ggAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      <  1:209876/250‑1 (MQ=255)
                                                                                                                         ggAAATGCGCAAAGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      >  1:582402/1‑250 (MQ=255)
                                                                                                                                             aTTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGcat                                                                                                                                                                                                      >  1:473155/1‑100 (MQ=255)
                                                                                                                                               ttACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTcgc                                                                                                                   >  1:283615/1‑183 (MQ=255)
                                                                                                                                               ttACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGt                                                                        <  1:632915/226‑1 (MQ=255)
                                                                                                                                               ttACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATa                                               >  1:733341/1‑251 (MQ=255)
                                                                                                                                                    cGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGCCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGAGGTAATGCGCCAAAGTCAGCAAGCCATAAACcg                                          >  1:506287/1‑251 (MQ=255)
                                                                                                                                                         gTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCcaacag                                                                                                                                                        <  1:701442/136‑2 (MQ=255)
                                                                                                                                                         gTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      <  1:728742/218‑1 (MQ=255)
                                                                                                                                                         gTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAAc                                            >  1:882176/1‑244 (MQ=255)
                                                                                                                                                                             cacaACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      >  1:701921/1‑198 (MQ=255)
                                                                                                                                                                             cacaACATTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAACCGATTCGATTTGATTAGGCGGAATccg                 >  1:390128/1‑249 (MQ=255)
                                                                                                                                                                               caacaTTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAACCGATTCGATTTGATTAGGCGGAATCCGtt               >  1:809807/1‑251 (MQ=255)
                                                                                                                                                                                  caTTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGc                                                       <  1:390954/208‑1 (MQ=255)
                                                                                                                                                                                  caTTTGCGCCACCAAATTGCCGCGCAGGCGACGTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAACCGATTCGATTTGATTAGGCGGAATCCGttctt            >  1:798858/1‑251 (MQ=255)
                                                                                                                                                                                             ccaaattgccgcgcaggcgacgTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAACCGATTCGATTTGATTAGGCGGAATCCGTTCTTTGTAAGAaat  >  1:130677/18‑250 (MQ=255)
                                                                                                                                                                                             ccaaattgccgcgcaggcgacgTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAACCGATTCGATTTGATTAGGCGGAATCCGTTCTTTGTAAGAaat  >  1:28453/18‑250 (MQ=255)
                                                                                                                                                                                                          caggcgacgTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAACcgatt                                       <  1:114989/196‑1 (MQ=255)
                                                                                                                                                                                                               gacgTTGTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTAATGCGCCAAAGTCAGCAAGCCATAAAcc                                           <  1:739850/191‑1 (MQ=255)
                                                                                                                                                                                                                     gTCCGCGTCGAATTTGACGAATTAATGCATCATCCTCAAATCCTCCTTCAGGCACCGCCCACACACAGCCCAACAGGGTTTCCCGCGTTTGTGCAAAATAAAATTGTTGTCGCGGTGCATCCAATAATCGACGTAAATCCAACGGCGAGGTGCGGTaa                                                                      >  1:97440/1‑158 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                                          
CAGCTCGGGGGTGTAACCGAAACTGACGGAAAGGAAATCCACCGCACGTTGTTGTTTAACTC‑GTTCGATTAAACGTTGCCCCAGTCCCCGCTGCTGCCAGTTGGGTTGCACGGCGATACGGGAAATGCGCAGTGAACGCAATTTACACGCCTGTTCCAATCCCGCCTGAAAACACAACATTTGCGCCACCAAATTACCACGTGGACGACGTTGTCCGCGCCGAATTTTACGAATTAATGCGTCGTCACCAAATTTGCCCTCAGGCACCGCCCACACACAACCCAACAAGGCCTCCTGCGTCTGCGCAAGATAAAATTGTTGTCGCGGCGCATCCAATAATCGACGTAAATCCAACGGCGAGGTTCGGTAATGCGCCAAGGTCAGCAAGCCATAAACCGATTCGATTTGATAAGGCGGAATCTGTTCTTCGTAAGAAAT  >  1/1474599‑1475036

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