Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 7,724 C→A W211C (TGG→TGT)  unknown ← FIG00031715: Predicted metal‑dependent phosphoesterases (PHP family)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*17,7240CA97.4% 133.3 / NA 39W211C (TGG→TGT) unknownFIG00031715: Predicted metal‑dependent phosphoesterases (PHP family)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (30/8);  total (30/9)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCGGTCGTTTTACAAGGTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCCCAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGACGCGCATAATGCGCTCTGGTGA  >  1/7481‑7965
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
gTGCGGTCGTTTTACAAGGTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACttt                                                                                                                                                                                                                                            >  1:588240/1‑251 (MQ=255)
         ttttACAAGGTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAAtt                                                                                                                                                                                                                                     >  1:86106/1‑249 (MQ=255)
          tttACAAGGTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATtct                                                                                                                                                                                                                                   >  1:161599/1‑250 (MQ=255)
                  gTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGtt                                                                                                                                                                                                                          >  1:68305/1‑251 (MQ=255)
                  gTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGtt                                                                                                                                                                                                                          >  1:620894/1‑251 (MQ=255)
                    aTTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGtttt                                                                                                                                                                                                                        >  1:670470/1‑251 (MQ=255)
                                ctgctgTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCAttt                                                                                                                                                                                                             >  1:21134/1‑250 (MQ=255)
                                     gTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGt                                                                                                                                                                                                       >  1:12966/1‑251 (MQ=255)
                                     gTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGt                                                                                                                                                                                                       >  1:555276/1‑251 (MQ=255)
                                               ggTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGG‑TTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTaa                                                                                                                                                                                             >  1:719103/1‑250 (MQ=255)
                                                            atTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGc                                                                                                                                                                                          >  1:5708/1‑241 (MQ=255)
                                                                  gACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAg                                                                                                                                                                                                                            >  1:738375/1‑201 (MQ=255)
                                                                   aCAAAAGGAGTGTTTCCCGGCGGGGTTTTTATTTTCCCTATTCCCCCCCACGCGAAAGAAAAGAAAAAACCCGACCCCATCGACCCCGCCAATTGGTGTTCCTCAGCCATTTTTCCCAAAAATTGAAATTGATGGGGACG‑TGGCCCGCAACC‑GAAA‑TTTCCCGTGGCACCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAacac                                                                                                                                                                         <  1:887178/251‑1 (MQ=255)
                                                                      aaaGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGc                                                                                                                                                                                          >  1:174095/1‑231 (MQ=255)
                                                                         ggAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAg                                                                                                                                                                                                        >  1:805364/1‑214 (MQ=255)
                                                                                            gTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGaata                                                                                                                                                >  1:804412/1‑251 (MQ=255)
                                                                                                     ttCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAAtg                                                                                                                                       >  1:551225/1‑251 (MQ=255)
                                                                                                      tCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATa                                                                                                                                                                                                         >  1:449233/1‑184 (MQ=255)
                                                                                                      tCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCt                                                                                                                                                         >  1:197474/1‑232 (MQ=255)
                                                                                                      tCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGaataa                                                                                                                                               >  1:683686/1‑242 (MQ=255)
                                                                                                      tCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAAtgg                                                                                                                                      >  1:453402/1‑251 (MQ=255)
                                                                                                      tCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAAtgg                                                                                                                                      >  1:64629/1‑251 (MQ=255)
                                                                                                      tCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAAtgg                                                                                                                                      >  1:713125/1‑251 (MQ=255)
                                                                                                          aGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAAtggtgg                                                                                                                                   >  1:407850/1‑250 (MQ=255)
                                                                                                           gTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAAtggtgg                                                                                                                                   >  1:36369/1‑249 (MQ=255)
                                                                                                                         cAAGAAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCa                                                                                                                   >  1:453056/1‑251 (MQ=255)
                                                                                                                                      aaTCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGaataa                                                                                                                                               >  1:46811/1‑210 (MQ=255)
                                                                                                                                                          gCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACAt                                                                                                     >  1:274525/1‑232 (MQ=255)
                                                                                                                                                               ttGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAAt                                                                             >  1:896627/1‑251 (MQ=255)
                                                                                                                                                                           cccttgtttccccagaatttaaaattgaTGGGGACG‑TTCCCCGAACCCGAAAA‑TTC‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAgg                                                                 <  1:121618/230‑1 (MQ=255)
                                                                                                                                                                            cctattttaacaagaaattaaaattgaCGGGGAGG‑TCGCCCGAACCCGAAAA‑TTT‑CCAGGCCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGc                                                                <  1:394231/236‑1 (MQ=255)
                                                                                                                                                                             cctttttaaccaacatttgaaattgaTCGGG‑CGTTTGCCCGC‑ACCGGAAATTTT‑CCAGGCCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAgg                                                                 <  1:823985/233‑1 (MQ=255)
                                                                                                                                                                                tttttCCCAGCAATAAAATTTGATGGGGACG‑TTCCCCGCAACCGAAAT‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAAt                                                                             <  1:599846/233‑1 (MQ=255)
                                                                                                                                                                                                          ggACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTg                                  >  1:139029/1‑251 (MQ=255)
                                                                                                                                                                                                                          acccgagaatttcgggggcagcgccGCCCCAACTTTTAAAATTCTCCGCTAGTTTTCGCTCCCATTTCATAGGCATGGTGTAACGCAACGGATGCGCCAACTCGGCAGCGCCACCGGCGTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCGTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGAcgc                     <  1:3570/228‑1 (MQ=255)
                                                                                                                                                                                                                                  attcccATGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGACGCGCATAATg             <  1:53743/247‑1 (MQ=255)
                                                                                                                                                                                                                                      ccgtGGCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGACGCGCATAATGCGct         <  1:653368/248‑1 (MQ=255)
                                                                                                                                                                                                                                          ggCATCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGACGCGCATAATGCGCTCTgg     >  1:498949/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ctCACCGCCACAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGACGCGCATAATGCGCTCTGGTGa  <  1:726061/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
GTGCGGTCGTTTTACAAGGTATTTTTAATGGACTGCTGTTCAGGCAAGGTCATTTTTTATATTTAAGACAAAAGGAATGTGTTCCGGCAGCGGTTTTAGTTTTCCTAGTTCCACCCACGCGCAAGGAAAGTGAAAATCCGATCCCATCGAACCCGCCAATTGGTGTTCCTCAGCCAGTTTTACCAGCAATTGAAATTGATCGGGACG‑TTGCCCGCAACCGGAAA‑TTT‑CCATGGCATCACCGCCCCAACTTTTAAAATTCTCCGCTAGTTTTCGCACCCATTTCATAGTCATGGTGTAACGCAACGGATGCGCCAACACGGCAGCGCCACCGGCTTGGTGAATAATGTCAATGGTGGTTGGAATGTCCGCCCATTGCGGCTTCACATAAGCAGACTTGCCCTGTCCCAAATAACGTTTAAAGGCTTGCCCGTCGTTAGAGACTTTGCCGATTTGCACTAAATGACGCGCATAATGCGCTCTGGTGA  >  1/7481‑7965

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