Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 460,337 G→T P21Q (CCG→CAG)  unknown ← Hydrogenase‑4 maturation protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1460,3370GT95.6% 173.8 / ‑4.0 68P21Q (CCG→CAG) unknownHydrogenase‑4 maturation protease
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (29/36);  total (30/38)
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

AAATGGGCTTTATTTTCGTCCAGCGCTTTCGCCAATAAATTCACATGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCGGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATCCAAACGTAAATTCGGGCGGGC‑ACAGGACCCGCCCCTACGGGAGAATGTAAGCCAATAACTTACCCCGAAAACGGCATACTAATGG  >  1/460097‑460558
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
aaaTGGGCTTTATTTTCGTCCAGCGCTTTCGCCAATAAATTCACATGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGttc                                                                                                                                                                                                                        <  1:776798/251‑1 (MQ=255)
    gggCTTTATTTTCGTCCAGCGCTTTCGCCAATAAATTCACATGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCg                                                                                                                                                                                                                    >  1:485379/1‑251 (MQ=255)
                        gcTTTCGCCAATAAATTCACATGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGcgttcgtt                                                                                                                                                                                               >  1:249557/1‑251 (MQ=255)
                         cTTTCGCCAATAAATTCACATGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTg                                                                                                                                                                                              <  1:896715/251‑1 (MQ=255)
                                         cacaTGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTa                                                                                                                                                                              >  1:333394/1‑251 (MQ=255)
                                            aTGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACgg                                                                                                                                                                           <  1:476386/251‑1 (MQ=255)
                                            aTGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACgg                                                                                                                                                                           <  1:730738/251‑1 (MQ=255)
                                               aCTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACggtgg                                                                                                                                                                        >  1:27547/1‑251 (MQ=255)
                                                  ggagtcgatcaCGATTTCGCCGAATTTCATTAACCCTGCAAATTTACGCTTTGCCTCCCCCTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACCCCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGTTCCGGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCa                                                                                                                                                                     >  1:67599/8‑251 (MQ=255)
                                                   gAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCAc                                                                                                                                                                                                                              <  1:297550/194‑1 (MQ=255)
                                                      tcaatcaCGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGAt                                                                                                                                                                 <  1:421132/251‑1 (MQ=255)
                                                       caatcaCGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATa                                                                                                                                                                <  1:878666/251‑1 (MQ=255)
                                                       caatcaCGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATa                                                                                                                                                                <  1:893594/251‑1 (MQ=255)
                                                        aatcaCGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATa                                                                                                                                                                >  1:790827/1‑250 (MQ=255)
                                                         atcaCGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATg                                                                                                                                                                   <  1:565989/246‑1 (MQ=255)
                                                                tttCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTAc                                                                                                                                                       <  1:25398/251‑1 (MQ=255)
                                                                tttCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTAc                                                                                                                                                       <  1:335917/251‑1 (MQ=255)
                                                                   cGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGc                                                                                                                                                    <  1:657230/251‑1 (MQ=255)
                                                                        aaTTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTAc                                                                                                                                               <  1:397482/251‑1 (MQ=255)
                                                                           ttCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATc                                                                                                                                            <  1:88375/251‑1 (MQ=255)
                                                                            tCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGagaag                                                                                                                                                                                         <  1:324278/205‑1 (MQ=255)
                                                                            tCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGAt                                                                                                                                                                 >  1:747891/1‑229 (MQ=255)
                                                                            tCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTc                                                                                                                                                            >  1:771524/1‑234 (MQ=255)
                                                                             cATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTa                                                                                                                                          >  1:278933/1‑251 (MQ=255)
                                                                             cATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTa                                                                                                                                          >  1:669668/1‑251 (MQ=255)
                                                                                tAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGc                                                                                                                                       <  1:886865/251‑1 (MQ=255)
                                                                                  aCCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACgg                                                                                                                                                                           >  1:663464/1‑213 (MQ=255)
                                                                                   cccTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACg                                                                                                                                                                            >  1:41433/1‑211 (MQ=255)
                                                                                   cccTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCaaa                                                                                                                                    >  1:396323/1‑250 (MQ=255)
                                                                                   cccTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCaa                                                                                                                                     >  1:711538/1‑250 (MQ=255)
                                                                                     cTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACag                                                                                                                                  >  1:503856/1‑251 (MQ=255)
                                                                                       gtagaTTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACagag                                                                                                                                >  1:158218/5‑251 (MQ=255)
                                                                                          aaTTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACAt                                                                                                                                             <  1:752152/235‑1 (MQ=255)
                                                                                                  cTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGc                                                                                                                                       >  1:508988/1‑233 (MQ=255)
                                                                                                   tttGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAg                                                                                                                    <  1:474127/251‑1 (MQ=255)
                                                                                                    ttGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGc                                                                                                                                                    <  1:201036/218‑1 (MQ=255)
                                                                                                               tCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTAc                                                                                                                                                       <  1:777993/204‑1 (MQ=255)
                                                                                                                 cGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCtg‑tg                                                                                                     <  1:361862/251‑1 (MQ=255)
                                                                                                                             gACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCAc                                                                                                                                                                                                                              >  1:365857/1‑120 (MQ=255)
                                                                                                                              aCAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATc                                                                                        <  1:143269/251‑1 (MQ=255)
                                                                                                                                 aaTTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGCTcccttgcccccccctttttca                                                                                       >  1:862138/1‑231 (MQ=255)
                                                                                                                                   ttgttgATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGc                                                                                                                                                    <  1:686028/187‑1 (MQ=255)
                                                                                                                                     gttgATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATcc                                                                                       <  1:482715/245‑1 (MQ=255)
                                                                                                                                     gttgATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGGCCTT‑Tgcccccccccttttcaaacaa                                                                                 >  1:399031/1‑238 (MQ=255)
                                                                                                                                             tCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑Tgcccgcc                                                                                               <  1:757883/229‑1 (MQ=255)
                                                                                                                                             tCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATc                                                                                        <  1:361607/236‑1 (MQ=255)
                                                                                                                                             tCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATCCaaa                                                                                    <  1:679448/240‑1 (MQ=255)
                                                                                                                                              ccTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑Tgcccgcccgaatttacgtttggataatggg                                                                        <  1:633978/251‑23 (MQ=255)
                                                                                                                                                                           aaGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTgaga                                                                                                                                                                                           <  1:254011/108‑1 (MQ=255)
                                                                                                                                                                                atcaatcaCACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTc                                                                                                                                                            <  1:532553/134‑1 (MQ=255)
                                                                                                                                                                                  caatcaCACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGGCCCGTTGCCC‑CCCCTTATCCAAACCTAAAATCAGGCGGGG‑CCAGGACCCGccccccccgggggaagtaaaa                                   >  1:799218/1‑234 (MQ=255)
                                                                                                                                                                                   aaccacaCCCACGTGGGGACCTATTGCAAAAGAGTAATACATTCCTTGCTGCGCCTGCTCTGGCACTTTTTCTTC‑TTTTCGAAAAAAGCTTTTGTAAAGAAAAAAGAATATGAGGGTGGTAATGATCCGCCCTCATGGCTTTACATCAAAGCAAGAGATTTGTTTTGGGGGGGGGGGGTCCTG‑TCCCCCCCCATTACCAAAACGAAATTTCGGG‑GGGCAACAGGCCCCGCCCCTACGGGAGAATGTAAGcc                                  <  1:876055/248‑1 (MQ=255)
                                                                                                                                                                                    atcaCACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCCG‑Tgcccgcc                                                                                               >  1:28545/1‑190 (MQ=255)
                                                                                                                                                                                    atcaCACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGGCCTG‑TGCCCGCCCCTTTTTCAAAAATAAAATCTGGCGGGG‑GCCGGACCCGCCCCCACGGGGGAAAGTaaacc                                  >  1:131283/1‑247 (MQ=255)
                                                                                                                                                                                        cacCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGt                                                                                                            >  1:683627/1‑174 (MQ=255)
                                                                                                                                                                                              caTGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATCCAAACGTAAATTCGGGCGGGC‑ACAGGAc                                                           <  1:565525/215‑1 (MQ=255)
                                                                                                                                                                                               aTGGTGCCCGATTGCAAAAGAGTAATTCATTACTTGCGGCCCCTCCCCCTGCATTTTTTCTTCTTTTTCA‑AAAAAGGGTTCGTTGAGAAGAAAGAAGATGGGGGTGGTAATGATACGTCTTTCCTCTTTTACATCAAAACAACAGATTTGTTTTTGGGGGGGGGGGCCCTG‑CCCCCCCCCTTTATCAAAACTAAAATTGGGGGGGGA‑AAAGGACCCGCCCCTACGGGAGAATGTAAGCCAATAACTTAccc                      <  1:125069/251‑1 (MQ=255)
                                                                                                                                                                                                tGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑Tgc                                                                                                    <  1:70885/173‑1 (MQ=255)
                                                                                                                                                                                                                aaGAGTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTAc                                                                                                                                               >  1:461148/1‑115 (MQ=255)
                                                                                                                                                                                                                  gagTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGagaagatagaag                                                                                                                                                                                  >  1:268016/1‑78 (MQ=255)
                                                                                                                                                                                                                  gagTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATCCAAACGTAAATTCGGGCGGGC‑ACAGGAcc                                                          <  1:272364/196‑1 (MQ=255)
                                                                                                                                                                                                                   aatAATACTTTCCTTGCTGCCCCTCCTCCGCCTTTTTTTCTTGTTTTTCG‑AAAAGGCTTTGGTTGAAAAAAAAGAAGATTACGGTGGTAATGATACGTCCTTCCTTTTTTACATCAAAACAACAGGTTTGTTTTGGGGGGGGGGGGTCCGG‑TCCCCGCCATTTATCAAAACGAAAATTCGGGGGGGC‑ACAGGCCCCGCCCCTACGGGAGAATGTAAGCCAATAACTTACCCCGAAAACGGCATACTAATgg  <  1:189423/249‑1 (MQ=255)
                                                                                                                                                                                                                    gTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACgg                                                                                                                                                                           >  1:250003/1‑83 (MQ=255)
                                                                                                                                                                                                                    gTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGAt                                                                                                                                                                 >  1:640164/1‑93 (MQ=255)
                                                                                                                                                                                                                    gTAATACATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTcc                                                                                                          <  1:173172/148‑1 (MQ=255)
                                                                                                                                                                                                                         aCATTACTTGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCCTTATc                                                                                        >  1:852399/1‑160 (MQ=255)
                                                                                                                                                                                                                                ttGCTGCGCCTGCTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATc                                                                                                                                            >  1:174711/1‑102 (MQ=255)
                                                                                                                                                                                                                                            cTCCTGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACAt                                                                                                                                             <  1:614357/89‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                
AAATGGGCTTTATTTTCGTCCAGCGCTTTCGCCAATAAATTCACATGACTGGAATCAATCACGATTTCGCCGAATTTCATTAACCCTGTAAATTTACGCTTTGCTTCCCCTTCCGGGAAGGCGGCGACAAATTGTTGATAATCCTCGTAGTCGCAAATCAGCAATTTTTTGAAGCAATCAATCACACCCACATGGTGACCGATTGCCAAAGAGTAATACATTACTTGCTGCGCCTGCTCCGGCACTTGTTCTTCGTTTTCG‑ACAAAGCGTTCGTTGAGAAGATAGAAGATTACGGTGGTCATGATACGTCCTTACTGCTTTACATCTAAGCAACAGAGTTGCTTTTGTAGGGGCGGGTCCTG‑TGCCCGCCCATTATCCAAACGTAAATTCGGGCGGGC‑ACAGGACCCGCCCCTACGGGAGAATGTAAGCCAATAACTTACCCCGAAAACGGCATACTAATGG  >  1/460097‑460558

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