Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 575,748 G→A A299V (GCC→GTC)  unknown ← Probable zinc protease pqqL (EC 3.4.99.‑)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1575,7480GA98.7% 189.0 / NA 75A299V (GCC→GTC) unknownProbable zinc protease pqqL (EC 3.4.99.‑)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/46);  minor base C (0/1);  total (28/47)
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

CTCGTCGGCGGTAAATCCATGTTGTTTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGGCCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATCCGGTCTATACCATT  >  1/575501‑575991
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
cTCGTCGGCGGTAAATCCATGTTGTTTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCg                                                                                                                                                                                                                                                  <  1:576008/251‑1 (MQ=255)
    tcgGCGGTAAATCCATGTTGTTTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAGCCAGTTGTGTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTc                                                                                                                                                                                                                                              <  1:156458/251‑1 (MQ=255)
       gcggTAAATCCATGTTGTTTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCtt                                                                                                                                                                                                                                            <  1:670239/250‑1 (MQ=255)
                     ttgttTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATcttctt                                                                                                                                                                                                                             >  1:39641/1‑251 (MQ=255)
                       gttTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCa                                                                                                                                                                                                                           >  1:353431/1‑251 (MQ=255)
                       gttTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCa                                                                                                                                                                                                                           >  1:738475/1‑251 (MQ=255)
                           aaTTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCaaaaa                                                                                                                                                                                                                       <  1:403888/251‑1 (MQ=255)
                                         aaGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATGGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCttt                                                                                                                                                                                                                                           <  1:427376/217‑1 (MQ=255)
                                            cAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGc                                                                                                                                                                                                       <  1:572995/250‑1 (MQ=255)
                                                  gTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGata                                                                                                                                                                                                >  1:296448/1‑251 (MQ=255)
                                                       gTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTgg                                                                                                                                                                                            <  1:289432/250‑1 (MQ=255)
                                                       gTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTgg                                                                                                                                                                                            <  1:772115/250‑1 (MQ=255)
                                                       gTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGt                                                                                                                                                                                           >  1:350516/1‑251 (MQ=255)
                                                              tttttAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGt                                                                                                                                                                                    >  1:337325/1‑251 (MQ=255)
                                                                   aaaTAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCgg                                                                                                                                                                                >  1:559792/1‑250 (MQ=255)
                                                                       aatTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATa                                                                                                                                                                           >  1:173114/1‑251 (MQ=255)
                                                                              aTCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTa                                                                                                                                                                    <  1:587078/251‑1 (MQ=255)
                                                                                 gATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCaaaaa                                                                                                                                                                                                                       >  1:570439/1‑197 (MQ=255)
                                                                                 gATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTgg                                                                                                                                                                                            <  1:311321/224‑1 (MQ=255)
                                                                                 gATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTAcc                                                                                                                                                                  >  1:773880/1‑250 (MQ=255)
                                                                                 gATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTAcc                                                                                                                                                                  >  1:168419/1‑250 (MQ=255)
                                                                                      aTTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCaa                                                                                                                                                            >  1:672762/1‑251 (MQ=255)
                                                                                              gAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTc                                                                                                                                                                                  <  1:350036/221‑1 (MQ=255)
                                                                                              gAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTGCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGaa                                                                                                                                                     <  1:8862/250‑1 (MQ=255)
                                                                                              gAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACAATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAAt                                                                                                                                                    >  1:527346/1‑251 (MQ=255)
                                                                                               aaaaGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATc                                                                                                                                                   <  1:712856/251‑1 (MQ=255)
                                                                                               aaaaGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATc                                                                                                                                                   <  1:474683/251‑1 (MQ=255)
                                                                                                         ttGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGc                                                                                                                                                                                                       <  1:723696/189‑1 (MQ=255)
                                                                                                         ttGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATg                                                                                                                                          >  1:68483/1‑250 (MQ=255)
                                                                                                         ttGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCAATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAGAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATg                                                                                                                                          >  1:439624/1‑250 (MQ=255)
                                                                                                               ggTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTAtc                                                                                                                                                                                                                                  <  1:635543/156‑1 (MQ=255)
                                                                                                               ggTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTa                                                                                                                                                                    <  1:526313/218‑1 (MQ=255)
                                                                                                                gtttgtttACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAAt                                                                                                                                   >  1:667481/1‑250 (MQ=255)
                                                                                                                gtttgtttACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAAt                                                                                                                                   >  1:788745/1‑250 (MQ=255)
                                                                                                                 tttgtttACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTCCCCGCCAATGATGAATCAGCGGAATGCTAAAATc                                                                                                                                  <  1:759950/250‑1 (MQ=255)
                                                                                                                  ttgtttACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTAtctt                                                                                                                                                                                                                                <  1:409091/155‑1 (MQ=255)
                                                                                                                          cAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCt                                                                                                                        <  1:681188/251‑1 (MQ=255)
                                                                                                                            aGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTaa                                                                                                                      >  1:279679/1‑251 (MQ=255)
                                                                                                                             gATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTaa                                                                                                                      >  1:537159/1‑250 (MQ=255)
                                                                                                                                     cGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTtg                                                                                                             >  1:45262/1‑251 (MQ=255)
                                                                                                                                     cGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTtg                                                                                                             <  1:283523/251‑1 (MQ=255)
                                                                                                                                     cGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTtg                                                                                                             <  1:777114/251‑1 (MQ=255)
                                                                                                                                     cGATAGAAATTCGCCGAATCCAACCAGCTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTAc                                                                                                                                                                   <  1:423296/197‑1 (MQ=255)
                                                                                                                                                               gttgttTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTgg                                                                                                                                                                                            <  1:38567/146‑1 (MQ=255)
                                                                                                                                                               gttgttTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATAc                                                                                                                                                                          <  1:407616/164‑1 (MQ=255)
                                                                                                                                                                ttgttTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGCCTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAgttg                                                                                  <  1:740278/251‑1 (MQ=255)
                                                                                                                                                                ttgttTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAGGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAgttg                                                                                  <  1:265758/251‑1 (MQ=255)
                                                                                                                                                                         cTTTTTTCCCATTTTTGCAGACTTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTCCCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGata                                                                          <  1:320521/250‑1 (MQ=255)
                                                                                                                                                                             tttCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAAc                                                                     <  1:202907/251‑1 (MQ=255)
                                                                                                                                                                             tttCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAAc                                                                     <  1:436943/251‑1 (MQ=255)
                                                                                                                                                                             tttCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAAc                                                                     <  1:623317/251‑1 (MQ=255)
                                                                                                                                                                             tttCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAAc                                                                     <  1:663803/251‑1 (MQ=255)
                                                                                                                                                                               tccaaTTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTTTTTTTGTTTATCTTCTTCAAAAAAGCTAAACTCCAGTGCCGGAATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAgttg                                                                                  <  1:679875/232‑1 (MQ=255)
                                                                                                                                                                                       ttGCAGACGTAAATTGACTAAACGGGTCATAATTTGTGGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAAttt                                                           <  1:536061/251‑1 (MQ=255)
                                                                                                                                                                                                      gACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGatat                                            >  1:400206/1‑251 (MQ=255)
                                                                                                                                                                                                      gACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGatat                                            >  1:879084/1‑251 (MQ=255)
                                                                                                                                                                                                                 tCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGGACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCgacgac                                  <  1:59857/250‑1 (MQ=255)
                                                                                                                                                                                                                    caATTTGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCAc                                         >  1:237718/1‑240 (MQ=255)
                                                                                                                                                                                                                       tttGTTGAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAgg                           <  1:258655/251‑1 (MQ=255)
                                                                                                                                                                                                                        ttgttgAATGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATAGTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAggg                          <  1:794937/251‑1 (MQ=255)
                                                                                                                                                                                                                              aaTGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCC‑ATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACatat                    <  1:409515/250‑1 (MQ=255)
                                                                                                                                                                                                                              aaTGAGATCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACGCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACatat                    <  1:645013/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gagaTCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCATCTCCAGTGCCGGGATATTGGTGCCTTGTTCGGAAATACTAGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAAc                                                                     <  1:690303/199‑1 (MQ=255)
                                                                                                                                                                                                                                 gagaTCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACAtattat                 >  1:397445/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gagaTCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACAtattat                 >  1:404389/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aTCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGTGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATcc               <  1:55083/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aTCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATcc               <  1:578753/250‑1 (MQ=255)
                                                                                                                                                                                                                                    aTCTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATc                <  1:21490/249‑1 (MQ=255)
                                                                                                                                                                                                                                      cTTGTTTGTAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAgttg                                                                                  <  1:33895/181‑1 (MQ=255)
                                                                                                                                                                                                                                        tgtttgtAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTa                                                                                      >  1:507002/1‑175 (MQ=255)
                                                                                                                                                                                                                                         gtttgtAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTtgtg                                                                                                           >  1:266086/1‑153 (MQ=255)
                                                                                                                                                                                                                                            tgtAGCCAATGCCCGTTTCTTTTTTTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCGGGGATATTGGTGCCTTTTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATCCGGTCTATAc      <  1:869256/251‑1 (MQ=255)
                                                                                                                                                                                                                                            tgtAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATCCGGTCTATAc      <  1:207456/251‑1 (MQ=255)
                                                                                                                                                                                                                                              tAGCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATCCGGTCTATAcc     <  1:562032/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gCCAATGACCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATCCGGTCTATACCAtt  >  1:535395/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
CTCGTCGGCGGTAAATCCATGTTGTTTAATTTCGGCGATAAAAGCAAAAAGTGCGGTGATATTTTTTAAATAATTGGTATCGATTAATTGCAAGGAAAAGACACTTTGCAGGGTTTGTTTACCAAGATGAGAACGATAGAAATTCGCCGAATCCAACCAGTTGTTTTGGCTTTTTTCCCATTTTTGCAGACGTAAATTGACTAAACGGGTCACAATTTGTTGAATGAGATCTTGTTTGTAGCCAATGGCCGTTTCTTTTTGTTTATCTTCTTCAAAAAAGCTCAACTCCAGTGCCGGGATATTGGTGCCTTGTTCGGCAATACTTGCTACCCGCCAATGATGAATCAGCGGAATGCTAAAATCAATTTTATCTAACGGGCGTTGTGTGCGAGAACTCGGTTTATCTAGTTGTTGTGATATTAACTGTGCAATTTTGTGCGCATCGATATCACCGACGACAATAAGGGACATATTATCCGGTCTATACCATT  >  1/575501‑575991

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