Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 887,565 G→A 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
*1887,5650GA98.6% 196.7 / NA 70R420R (CGG→CGA) unknowntRNA cytosine(34) acetyltransferase (EC 2.3.1.193)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (33/36);  total (33/37)
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

ATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTG  >  1/887347‑887784
                                                                                                                                                                                                                          |                                                                                                                                                                                                                             
attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                                                                                                                                                                                                <  1:130677/250‑18 (MQ=255)
attTCTTACAAAGAACGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctgcgcggcaatttgg                                                                                                                                                                                                <  1:28453/250‑18 (MQ=255)
          aagaaCGGATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                                                                                     <  1:798858/251‑1 (MQ=255)
             aacggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtgttg                                                                                                                                                                                  <  1:809807/247‑1 (MQ=255)
               cggATTCCGCCTAATCAAATCGAATCGGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                                                                                                <  1:390128/249‑1 (MQ=255)
                                     aatcgGTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtcgcctg                                                                                                                                                                                                             >  1:114989/1‑196 (MQ=255)
                                        cgGTTTATGGCTTGCTGACTTTGGCGCATTACCTCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGGCGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCg                                                                                                                                                       <  1:506287/251‑1 (MQ=255)
                                         ggTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAAcgtc                                                                                                                                                                                                                  >  1:739850/1‑191 (MQ=255)
                                          gTTTATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                                            <  1:882176/244‑1 (MQ=255)
                                             tATGGCTTGCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                                                                                                  <  1:733341/251‑1 (MQ=255)
                                                     gCTGACTTTGGCGCATTACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAAtg                                                                                                                                                                                     >  1:390954/1‑208 (MQ=255)
                                                                    ttacagcaCCTCGCCGTTGGATTTACGACGATTATTGGATGCACCGAGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGAGTGTGGGCGGAGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGGAAATTGGGTCCTTTGCGCATTTcc                                                                                                                            >  1:210964/8‑250 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGAc                                                                                                                                                                                                                        <  1:97440/158‑1 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTtgtg                                                                                                                                                                                <  1:701921/198‑1 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                                            >  1:728742/1‑218 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                            >  1:209876/1‑250 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                            >  1:698893/1‑250 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                            <  1:582402/250‑1 (MQ=255)
                                                                    ttACCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTcc                                                                                                                            >  1:51847/1‑250 (MQ=255)
                                                                      aCCGCACCTCGCCGTTGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                                                                                                  >  1:632915/1‑226 (MQ=255)
                                                                                     tGGATTTACGTCGATTATTGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGccc                                                                                                          <  1:158901/251‑1 (MQ=255)
                                                                                                  attattGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCagcagc                                                                                             <  1:69320/251‑1 (MQ=255)
                                                                                                     attGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgag                                                                                          >  1:232337/1‑249 (MQ=255)
                                                                                                     attGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgag                                                                                          >  1:483325/1‑249 (MQ=255)
                                                                                                      ttGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGcagcag                                                                                              >  1:789812/1‑246 (MQ=255)
                                                                                                      ttGGATGCACCGCGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgagg                                                                                         >  1:781539/1‑248 (MQ=255)
                                                                                                             caccgcGACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTgggg                                                                                  <  1:759887/249‑1 (MQ=255)
                                                                                                                 gcgACAACAATTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTaa                                                                                                                                                  <  1:283615/183‑1 (MQ=255)
                                                                                                                        caaTTTTATTTTGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                       <  1:143521/251‑1 (MQ=255)
                                                                                                                          attttattttGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGTCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGGATCGCCGTGc                                                                                                               >  1:615681/1‑209 (MQ=255)
                                                                                                                          attttattttGCACAAACGCGGGAAACCCTGTTGGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCg                                                                     <  1:445071/251‑1 (MQ=255)
                                                                                                                               attttgcacaaacgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGAAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                       <  1:9417/220‑1 (MQ=255)
                                                                                                                                         aacgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaac                                                    <  1:144890/237‑1 (MQ=255)
                                                                                                                                           cgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAAt                                                                       >  1:790218/13‑232 (MQ=255)
                                                                                                                                           cgcgggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                  <  1:585054/239‑1 (MQ=255)
                                                                                                                                               ggaaaccctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGt                                               <  1:91252/242‑1 (MQ=255)
                                                                                                                                                      ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAAc                                                                                                                                                            >  1:701442/2‑136 (MQ=255)
                                                                                                                                                      ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTggg                                       >  1:308059/2‑250 (MQ=255)
                                                                                                                                                      ctgttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTggg                                       >  1:596187/2‑250 (MQ=255)
                                                                                                                                                        gttgGGCTGTGTGTGGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGAtt                                     >  1:386119/1‑251 (MQ=255)
                                                                                                                                                                gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTa                                                                                                                        >  1:833883/1‑162 (MQ=255)
                                                                                                                                                                gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑cagattaaac                                                        <  1:539036/225‑10 (MQ=255)
                                                                                                                                                                gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                  <  1:471552/230‑1 (MQ=255)
                                                                                                                                                                gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGt                                               >  1:297852/1‑233 (MQ=255)
                                                                                                                                                                gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGAttttctttc                              >  1:892622/1‑250 (MQ=255)
                                                                                                                                                                gtgtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGGGGAttttctttc                              >  1:30271/1‑250 (MQ=255)
                                                                                                                                                                  gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTTCCGTATCGCCGTGCAGCCCAACTGGCCGCAGCGAGGGCTGGGGCCACGTTTAATCGAACCGA‑‑TTAAACAACAACCTGCGGTGGATTTTCTTTCCGt                           >  1:16077/1‑251 (MQ=255)
                                                                                                                                                                  gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAaacaacaa                                                  <  1:437368/228‑1 (MQ=255)
                                                                                                                                                                  gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACg                                                >  1:323252/1‑230 (MQ=255)
                                                                                                                                                                  gtgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGt                           >  1:598662/1‑251 (MQ=255)
                                                                                                                                                                   tgtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGATCTTTGCGCATTTCCCGCATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCCGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGACTTTCTTTCCGt                           >  1:85052/1‑250 (MQ=255)
                                                                                                                                                                    gtgGGCGGTGCCTGAAGGAGGATTTGAGGATGATGCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGc                                                                                                                    >  1:830742/1‑162 (MQ=255)
                                                                                                                                                                                 gaaggaggatttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTaca             >  1:338311/22‑250 (MQ=255)
                                                                                                                                                                                        gatttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGGTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaa     >  1:770770/15‑250 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGc                                                                                                               <  1:147719/133‑1 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:578174/238‑3 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:217141/238‑3 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:149617/238‑3 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:836462/238‑3 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:543915/238‑3 (MQ=255)
                                                                                                                                                                                         atttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGGTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:446276/238‑3 (MQ=255)
                                                                                                                                                                                           ttgaggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGAactg  <  1:305880/240‑1 (MQ=255)
                                                                                                                                                                                               ggatgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑cagattaaac                                                        >  1:487161/8‑185 (MQ=255)
                                                                                                                                                                                                  tgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCgagggc                                                                                       <  1:227295/157‑6 (MQ=255)
                                                                                                                                                                                                  tgatgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCg                            <  1:390200/214‑1 (MQ=255)
                                                                                                                                                                                                    atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAAt                                                                                                                                                <  1:473155/100‑1 (MQ=255)
                                                                                                                                                                                                    atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑cagattaaaca                                                       >  1:325522/3‑180 (MQ=255)
                                                                                                                                                                                                    atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACAcccc         <  1:848980/233‑1 (MQ=255)
                                                                                                                                                                                                    atgCATTAATTCGTCAAATTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    >  1:366912/3‑238 (MQ=255)
                                                                                                                                                                                                                gtcaaaTTCGACGCGGACAACGTCGCCTGCGCGGCAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCTTTGCGCATTTCCCGTATCGCCGTGCAGCCCAACTGGCAGCAGCGAGGGCTGGGGCAACGTTTAATCGAA‑C‑AGATTAAACAACAACGTGCGGTGGATTTTCTTTCCGTCAGTTTCGGTTACACCCCCGaac    <  1:143746/225‑3 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                             
ATTTCTTACGAAGAACAGATTCCGCCTTATCAAATCGAATCGGTTTATGGCTTGCTGACCTTGGCGCATTACCGAACCTCGCCGTTGGATTTACGTCGATTATTGGATGCGCCGCGACAACAATTTTATCTTGCGCAGACGCAGGAGGCCTTGTTGGGTTGTGTGTGGGCTGTGCCTGAGGGCGAATTTGGTGACGACGCATTAATTCGTAAAATTCGGCGCGGACAACGTCGTCCACGTGGTAATTTGGTGGCGCAAATGTTGTGTTTTCAGGCGGGATTGGAACAGGCGTGTAAATTGCGTTCACTGCGCATTTCCCGTATCGCCGTGCAACCCAACTGGCAGCAGCGGGGACTGGGGCAACGTTTAATCGAA‑CGAG‑TTAAACAACAACGTGCGGTGGATTTCCTTTCCGTCAGTTTCGGTTACACCCCCGAGCTG  >  1/887347‑887784

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: ‑