Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,239 G→T A74E (GCG→GAG)  yaaA ← DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,2390GT87.7% 173.6 / ‑1.9 74A74E (GCG→GAG) yaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base T (33/31);  minor base C (0/9);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 3.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTAT  >  NC_000913/5992‑6483
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
gTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCg                                                                                                                                                                                                                                                   <  1:109470/251‑1 (MQ=255)
  tCAGCTTGTTGGTTATGCTACCTCCCCAGAAGTGATACAGATCTTTCCTCCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCAGTAAAGCGAGCGGGCGCAGTTCGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCACAATCGAAATCGTCTTCTCTGAAGGTTTCGGCCTTCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCgg                                                                                                                                                                                                                                                  <  1:125407/250‑1 (MQ=255)
   cAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCg                                                                                                                                                                                                                                                <  1:487481/251‑1 (MQ=255)
    aGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCAGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGt                                                                                                                                                                                                                                               >  1:205916/1‑251 (MQ=255)
     gCTTTTTGGGGATGATATCTCCCCAGAATTGATACAGATCTTGCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCATACGATAATGCTGCATTAAATCGAGCGGGGGCAGTACGCCATACAAGCCGGAAAGCATTGGCAAAGGCTGTG‑GTGCAAAATCGAAATCGTCTTCGCTGAAGGTGGCGGCCTGCAAGCCGGTTTAGACATCACCTTTAAACGCCAGAATCGCCTGGCTGGCAT‑TCTCCGGCGTg                                                                                                                                                                                                                                              <  1:610482/251‑1 (MQ=255)
        tgttgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGaa                                                                                                                                                                                                                                            <  1:61594/250‑1 (MQ=255)
        tgttgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGaa                                                                                                                                                                                                                                            >  1:403872/1‑250 (MQ=255)
            ggTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAACGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTg                                                                                                                                                                                                                                       >  1:237778/1‑251 (MQ=255)
                      tctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTc                                                                                                                                                                                                                             >  2:247606/1‑251 (MQ=255)
                         ctccAGAATTGGTACAGAACTTTTCCTCGTGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCAGTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTTGCGGGCGG‑TCTCCGGCGTTAAAGCTGGCTGACAGTCa                                                                                                                                                                                                                            <  1:501206/247‑1 (MQ=255)
                         ccccAGAATTGATACAGATCTTTCCATCGTGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCAGTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAAGCGCCGGGCGTGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCAt                                                                                                                                                                                                                           <  2:64932/250‑1 (MQ=255)
                           ccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGaa                                                                                                                                                                                                                        <  2:433509/251‑1 (MQ=255)
                              aaaTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAagc                                                                                                                                                                                                                     <  2:173362/250‑1 (MQ=255)
                                       cAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCCGGCTGCCAGTCATGAAAGCGAGCGGCGt                                                                                                                                                                                                            <  1:247757/251‑1 (MQ=255)
                                        aGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCGTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGtt                                                                                                                                                                                                           <  1:290958/251‑1 (MQ=255)
                                             tttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAc                                                                                                                                                                                                      >  1:287054/1‑251 (MQ=255)
                                              tgcccTCGTGCAGTCTCAATGCGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGTAATCGAGCGTCGTTGAtgc                                                                                                                                                                                                      <  2:371235/248‑3 (MQ=255)
                                                cccgcGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAccc                                                                                                                                                                                                    <  1:400483/246‑1 (MQ=255)
                                                cccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCg                                                                                                                                                                                                   >  2:218781/1‑251 (MQ=255)
                                                                 gACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCCCCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATg                                                                                                                                                                                  <  1:179717/251‑1 (MQ=255)
                                                                 gACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGTCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTgctg                                                                                                                                                                                  >  1:40758/1‑248 (MQ=255)
                                                                 gACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGACATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATg                                                                                                                                                                                  >  1:386568/1‑251 (MQ=255)
                                                                  aCGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATgc                                                                                                                                                                                 >  1:497460/1‑251 (MQ=255)
                                                                                ccAGACGATAAGGCTACATTAAGTCGGGAGGGCGCAGTACGCCAAACAAGCCTGAAAGAATCCGCAAA‑GGTGTTTGGGAAAAATCGAAAGCGTCTGGTCTGAGGGTTTCTGCCGGCAAGCCGGTTTAGCCATCTCCTTGAACCCCTAGAATCGTCTGGCGGG‑ATGTCTCCGTCGTAAAATCGGGCTGCCAGTCATGAAAGCGAGCGGCTTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTa                                                                                                                                                                   <  2:33707/251‑1 (MQ=255)
                                                                                   gACGATAAGGCTGCATAAAATCGAGCGGGCGAAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCGGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATc                                                                                                                                                                <  1:502037/251‑1 (MQ=255)
                                                                                       aTAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCg                                                                                                                                                            <  1:287367/251‑1 (MQ=255)
                                                                                         aagtgtccATTAAATCGAGCGGGCGCAGGACTCCATACAAGCCGGAAAGCTTTCGCAAATGCTGTT‑GGGAAAAACCGAAATCGTCTTCGCTGAAGGTTCCGGCCGGCAGGCGGGTGTAGACATCCCCTTTAAACGCCAGAATCGCCGGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCTTGAAAGCGAGCGTCGTTGAGACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGa                                                                                                                                                          <  2:570091/244‑1 (MQ=255)
                                                                                           ggCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAgg                                                                                                                                                        >  1:192069/1‑251 (MQ=255)
                                                                                             cTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTAGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGGGCGGCGTTGATACCCGCCATATTGGCGCTGATGCGCATAAGCGGGCTAATCTGCGGAGGCg                                                                                                                                                      >  2:353811/1‑251 (MQ=255)
                                                                                                 aTTAAACCGAGCGGGCGCAGTACGCCAAACAAGCCGAAAAGCATTCGCAAATGGGGTT‑GGGCAAAATCGAAATCGTCTTCGCGGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCATAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGTGCGGCGTTGATCCCCGCCATTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAg                                                                                                                                                  <  1:30890/251‑1 (MQ=255)
                                                                                                        cGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCcgc                                                                                                                                           >  2:460336/1‑251 (MQ=255)
                                                                                                         gggcgggcgcAGTACCCCATCCAACCCGGAAAGCTTCCGCAAATGCTGTT‑GGGCAAAATCGAAACCGTCTTCGCTGAAGGTTTCGGCCGGCAGGCCGGTGTAGACATCACCTTTAACCTCAAGAATCGCCTGGCGGGCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGATCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCcgcg                                                                                                                                          <  2:403872/249‑1 (MQ=255)
                                                                                                            cgggcgcAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCgcgctt                                                                                                                                       >  1:261807/1‑249 (MQ=255)
                                                                                                                cgcaGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATg                                                                                                                                   >  1:95768/4‑251 (MQ=255)
                                                                                                                   aGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGAt                                                                                                                                >  1:573374/1‑251 (MQ=255)
                                                                                                                        cccATACACGCCGGAAGGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAACCCTCTTCGCTGAATGTTCCGCCCTGCAGGCCGGTTTAGACCTCCCCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGGGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAAct                                                                                                                           <  2:172034/250‑1 (MQ=255)
                                                                                                                                  ccGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAAtt                                                                                                                 >  1:112450/1‑251 (MQ=255)
                                                                                                                                     gAAAGCATTCGCAAATGTTTTT‑GGGAAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCGCCGGCGGGAAATCGGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTc                                                                                                              <  2:210960/251‑1 (MQ=255)
                                                                                                                                     gAAAGCATTCGCAAATGCTGTT‑GGGCAAACTCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTc                                                                                                              >  2:618735/1‑251 (MQ=255)
                                                                                                                                         gCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAAc                                                                                                          >  2:419146/1‑251 (MQ=255)
                                                                                                                                           aTTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGCTCAACTGCTGGGAATTGTCTAAc                                                                                                          >  1:470539/1‑249 (MQ=255)
                                                                                                                                            tccGCAAATGCGGTT‑GGGCAAAATCGATACCGTCTTGGCTGAAGGTGTCGGCCGGCAGGCCGGTGTAGACACCACCTTTAAACGCCAGAACCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCTTTGATTCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAg                                                                                                        <  2:454560/248‑1 (MQ=255)
                                                                                                                                             tCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGGCTAACAGCt                                                                                                      >  2:179280/1‑251 (MQ=255)
                                                                                                                                              cGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAACTCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGCTACCCGCCAGTTTGACGCTGATGCGAATCAGCGTGCTACTCTGCGGAGGCGTCAGTTTCCGCGCTTCATGAATCAACTGCTGGGAATTGTCTAACAGCTc                                                                                                     >  1:585202/1‑251 (MQ=255)
                                                                                                                                               gCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTGGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTcc                                                                                                    >  2:139891/1‑251 (MQ=255)
                                                                                                                                                         gtcggccAAAATCGATATCGTCCTCGCTGAAGCTTTCGGCCTGCAACCCGGTGTAGCCATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCACGGATCACCTGCTGGGAATTGTCTAACAGCTCCGGCAGCGtc                                                                                           >  2:432194/7‑250 (MQ=255)
                                                                                                                                                            gggCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTCAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGCTCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGtat                                                                                          >  2:7221/1‑249 (MQ=255)
                                                                                                                                                            gggCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGtat                                                                                          >  1:212427/1‑249 (MQ=255)
                                                                                                                                                            gggCAAAACCAAAATCGTCTTCGTTGAAGGTTTCGGCCTGAAAGCCGGTTTAGACATCCCCTTTAAACGCCAGAATCGCATGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAg                                                                                        <  2:316412/251‑1 (MQ=255)
                                                                                                                                                              gCAAAATCGAAATCGTCTTCGCTGAAGTTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAgcg                                                                                      <  1:352362/251‑1 (MQ=255)
                                                                                                                                                               caaaatcaaaatCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCCCCTTTAAACGCCAGAAGCGCCGGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCGGCCAGCCATGAAAGCGAGCGGCTTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCt                                                                                                                                                               <  2:500455/177‑1 (MQ=255)
                                                                                                                                                               cAAAATCGCAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCt                                                                                                                                                               >  1:500455/1‑177 (MQ=255)
                                                                                                                                                                     cgaaatcgTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCgtggtg                                                                               >  1:363595/1‑251 (MQ=255)
                                                                                                                                                                      gaaatcgTCATCGCTGACGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGACAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGCAAGCGAGCGGCGTTGATACCCGCAAGTTTGTCGCTGATGCGCATAAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGaggtgg                                                                              >  2:454724/1‑251 (MQ=255)
                                                                                                                                                                          tcgTCTTCGCGGAAGGTTTCGGCCTGCAAGCCGGGGTAGACATCCCCTTTAAACGCCAGAATCGCCGGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTAATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCa                                                                           <  2:13716/250‑1 (MQ=255)
                                                                                                                                                                          ccgTCTTCGCGGAAGGTTTCGGCCGGCAAGCCGGTTTAGACACCACCTTTAAACGCCAGAATCGCCTGGGGGGCAT‑TCTCGGGCGTAAAATCTGGCTGCCAGTCATGAAAGCGGGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCaa                                                                          <  2:497460/250‑1 (MQ=255)
                                                                                                                                                                                tCGCTGAAGTTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACggg                                                                      <  1:218781/249‑1 (MQ=255)
                                                                                                                                                                                 cGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCCGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCtt                                                                   >  2:102138/1‑251 (MQ=255)
                                                                                                                                                                                      aaGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTa                                                              <  1:115106/251‑1 (MQ=255)
                                                                                                                                                                                             cggccggcAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGGGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGc                                                       <  2:173988/251‑1 (MQ=255)
                                                                                                                                                                                             cgcccggCAAGCCGGTGTAGACATCACCTTTAAACGCGAGAATCGCCGGGCGGGCTT‑TCCCCGGCGTAAAATCTGGCTGCCAGTCATGAAAGCGAGCGGGGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTCCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAg                                                        <  2:469461/244‑1 (MQ=255)
                                                                                                                                                                                               gCCTGCAAGCCGGGGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGCGGCGTCAGTTTCAGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCgg                                                     >  2:91775/1‑250 (MQ=255)
                                                                                                                                                                                                  ggaaggcgggTGTAGACACCACCTTAAGACCCAAGAAGCGCCGGGCGGGCAT‑TCTCCGGCGGAAAATCTGCCGGCCAGCCATGAAAGCGAGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTTGTAATCAAGCGtttt                                                  <  2:328308/243‑1 (MQ=255)
                                                                                                                                                                                                   gcaggCCGGTGTAGACATCCCCTTTAACCGCAAGAATCGCCTGGCGGGCTT‑TCTCCGGCGTGAAATCTGGCTGCCGGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTc                                                 <  2:386568/247‑1 (MQ=255)
                                                                                                                                                                                                   gCAAGCCGGTGCCGACATCAACTTTAAACGCCAGAATCGGCTGGCGGGCAT‑TCTCCGGCGAGACATCTGGCTGCCCGTCATGAAAGCGAGCGGCGATGATACCCGACAGTTTGTCGCTGCTGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTCTAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGtttcc                                                 >  2:45673/1‑249 (MQ=255)
                                                                                                                                                                                                      cgcgggTGTAGACACCACCTTTAAAACCAAAAATCGCCTGTCGGGTAT‑TCCCCGGCGTAAAATCGGGCTGCCTCTCATGAAAGCGAGCGGCGTTGATTCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCa                                              <  2:174509/247‑1 (MQ=255)
                                                                                                                                                                                                               agacaccacttgtAAACGCAAGAATCGCCGGGCGGGCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTGGGTAATCAAGCGTTTTCGCAGGTGAAATa                                     <  2:311456/239‑1 (MQ=255)
                                                                                                                                                                                                                  cATCACCTTTAAACCCAAGAATCGCCGGGGGGGCAT‑TCCCCGGCGTGAAATCGGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGa                                  <  2:185474/251‑1 (MQ=255)
                                                                                                                                                                                                                              aCGCCAGAATCGCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCAAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATcc                      <  1:432761/251‑1 (MQ=255)
                                                                                                                                                                                                                                      atagCCTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGACAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTTATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGACTTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGGGGTCACCGGGCTTTGGTCATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTg              >  2:399113/4‑251 (MQ=255)
                                                                                                                                                                                                                                         gCCTGGCGGGCAT‑TCTCCGGCGTGAAATCCGGCTGCAAGTCATGACAGCGAGCGTCGATGATACCCGCCCGTTTGTCGCTGATGCGCATCCGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGACTTGTCTAACAGCTCCGGCAGCGTATAGCGCGAGGTGGTCAACGGGCTTTGGTAATCAAGCGCTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGc             >  2:301475/1‑249 (MQ=255)
                                                                                                                                                                                                                                          ccTGGCGGGCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAgg          <  1:332856/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ggggCAT‑TCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAtt     <  1:179280/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCAT‑TCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTAt  <  1:187362/251‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
GTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTT‑GGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCAT‑TCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTAT  >  NC_000913/5992‑6483

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