Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,534 T→C S86S (AGT→AGC)  ispH → 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091326,5340TC95.9% 203.5 / NA 75S86S (AGT→AGC) ispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  major base C (31/39);  minor base A (3/0);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 9.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAC  >  NC_000913/26289‑26779
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
ttGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCcgcg                                                                                                                                                                                                                                                    <  2:387209/250‑1 (MQ=255)
      aaCCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGAc                                                                                                                                                                                                                                             >  1:356409/1‑251 (MQ=255)
        cctgcttgggttttGTGCCGGTGTAGACCGCGATATCAACATTGTTTAAATCGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCTAGAAGTGGTACATAACCGCTATGTGGTCGATATCTTGCGTGAGCGTGGGGCTATCTTGAGTGAGCAGATTAGCGAAGTACCGTGCGGCGCGATCCTCATTGTCTCCGCACACGGTGGTTCTCAGGCGGTACGGAACGAAGCAAAAAGCCGCGATTTGACgg                                                                                                                                                                                                                                           <  2:46346/241‑1 (MQ=255)
            ggTGGTTTTGGGGCCGGGGTAGATCGCGCTATAAGCATTGTTGAAAACGCGCTGGCCATTTACGGCTCACCGATATATGTACGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGgtggt                                                                                                                                                                                                                                       <  2:74769/250‑3 (MQ=255)
                            gggTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCg                                                                                                                                                                                                                       >  1:305220/1‑251 (MQ=255)
                             ggTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATACATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGTGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAATCCGCGATTTGACGGTGTTTGATGCCACCt                                                                                                                                                                                                                            <  1:486065/245‑1 (MQ=255)
                             ggTAGACCGCGCTATCAGCATTGTTGAAAAAGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCCCGTAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGTGCAGATTAGCGAAGTACCGAACGGCGCGATCATGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCCAAAAGCCGCGTTTTGACGGTGTTTGATGCCAcat                                                                                                                                                                                                                            >  2:486065/1‑243 (MQ=255)
                                aGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGTCGGTGTTTGATGCCACCTGTCCGCTgg                                                                                                                                                                                                                   <  1:370852/251‑1 (MQ=255)
                                     gcgcTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGAcc                                                                                                                                                                                                              >  1:43184/1‑251 (MQ=255)
                                       gcTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCa                                                                                                                                                                                                             >  1:496297/1‑250 (MQ=255)
                                        cTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCaa                                                                                                                                                                                                            >  1:558736/1‑250 (MQ=255)
                                             aGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGACGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAACCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAgggc                                                                                                                                                                                                       >  1:356068/1‑247 (MQ=255)
                                              gCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGAGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCat                                                                                                                                                                                                     >  1:122761/1‑251 (MQ=255)
                                                aTTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCatat                                                                                                                                                                                                   <  1:322902/251‑1 (MQ=255)
                                                aTTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCatat                                                                                                                                                                                                   <  1:236397/251‑1 (MQ=255)
                                                  tgttgAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATg                                                                                                                                                                                                  <  1:627362/250‑1 (MQ=255)
                                                  tgttgAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATg                                                                                                                                                                                                  <  1:499591/250‑1 (MQ=255)
                                                     tgcaaaCGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGTTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGAACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGCCGGTGTTTGTTGCCACCTGTCAGCTGGTGACCAAAGTGCATATGGa                                                                                                                                                                                                >  2:641691/4‑249 (MQ=255)
                                                     tgAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGa                                                                                                                                                                                                >  1:555783/1‑249 (MQ=255)
                                                      gAAAACGCGCTGGCCAGTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAg                                                                                                                                                                                              <  2:124569/250‑1 (MQ=255)
                                                              gcTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCcgcg                                                                                                                                                                                     >  1:287409/1‑251 (MQ=255)
                                                               cTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGACGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGAGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACAAAAGTGCATATGGAAGTCGCCCgcgc                                                                                                                                                                                    >  1:666689/1‑251 (MQ=255)
                                                                  gCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCACGATCCTGATTTTCTCCGCACACGGAGTTTATCCGGCGGTACGTAACGAAGCAAAAAGACGAGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCACGCGcc                                                                                                                                                                                   >  2:593580/1‑249 (MQ=255)
                                                                   ccATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGTTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAgc                                                                                                                                                                                <  1:260272/251‑2 (MQ=255)
                                                                   ccATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAgc                                                                                                                                                                                >  2:146492/1‑250 (MQ=255)
                                                                                  cGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCgaagag                                                                                                                                                                <  1:31484/251‑2 (MQ=255)
                                                                                  cGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCgaagag                                                                                                                                                                <  2:305220/251‑2 (MQ=255)
                                                                                     tatatGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCAcc                                                                                                                                                                                                                             >  1:573248/1‑188 (MQ=255)
                                                                                     tatatGCCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATATCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCAcc                                                                                                                                                                                                                             <  2:573248/188‑1 (MQ=255)
                                                                                        atGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCATGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCgaagagtcgatt                                                                                                                                                          >  1:56617/1‑244 (MQ=255)
                                                                                                 aCGAAGTGGTACATAACCGCTATGTGGGCGATAGGTTGCGTGAGCGAGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGACTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCCTATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCAtaggc                                                                                                                                                 >  1:272532/1‑247 (MQ=255)
                                                                                                          tACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCGCACGGTGGTTCTCAGGCGGTACTTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCcgcg                                                                                                                                                                                     <  1:629386/207‑1 (MQ=255)
                                                                                                          tACATAAACGCTATGTGGTCGATAGCTAGCGTGAGCGTGGGGCTATCTTTATTGCGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCAGCTGGTGACCAACGTGCATATGGAAGTCGCCcgcg                                                                                                                                                                                     >  2:629386/1‑207 (MQ=255)
                                                                                                                 ccGCTATGAGGACGATAGCTCGCGTGAGCGTGGGGCTATCTTTCTTGAGCAGATTAGCGAAGTACAGGACGGCGCGATCCTGATTGTCTTCGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGAGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGGTGCGAAAAGTCGATTCTCATCGGCCACGCCGGGCATCCgg                                                                                                                                 >  2:137683/1‑251 (MQ=255)
                                                                                                                      attgggTCGATAGCTTTCGTGAGCGTGGGGTTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTGCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCCCCCGGAAg                                                                                                                              <  2:479571/245‑1 (MQ=255)
                                                                                                                        gtgGTCGAGATCTTGCTTTAGCGTGGGGCTATCTTTATTGAGCAGATTATCGAAGTTCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTTCGTAACGAAGCAAAAAGCCTCGATTTGACGGTGTGTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGa                                                                                                                          <  2:541023/251‑1 (MQ=255)
                                                                                                                               aTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGCAGTACCGGCCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGCGACCAAAGGGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCGGTGGAAggtacg                                                                                                                   >  1:651350/1‑247 (MQ=255)
                                                                                                                                 aGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggtacgat                                                                                                                 <  1:504640/251‑7 (MQ=255)
                                                                                                                                         ggAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagt                                                                                                         <  1:290244/250‑1 (MQ=255)
                                                                                                                                           ggCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtac                                                                                                       <  1:1879/250‑1 (MQ=255)
                                                                                                                                                ggggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagta                                                                                                   >  1:96062/1‑250 (MQ=255)
                                                                                                                                                 gggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAc                                                                                                 >  1:58327/1‑251 (MQ=255)
                                                                                                                                                  ggCTATCTTTATTGGGCAGATTAGCGAAGTACCGGCCGGCGCGATCCTGATTTTCCCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGGCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAcc                                                                                                <  2:43184/251‑1 (MQ=255)
                                                                                                                                                                 gCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACTGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGt                                                                                 <  1:51101/251‑1 (MQ=255)
                                                                                                                                                                    gATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATc                                                                              <  1:233382/251‑1 (MQ=255)
                                                                                                                                                                     tttAGAGAAGAACCGGACGGCGCGATCCTGCTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAAAGAAGCAAAAAGCCGCGATTTGACGGTGTTTGCTGCCACCTGTCCGCTGGTGACCAACGTGAATCTGGAAGTAGCACGCGCCAG‑CCGTCGTGGCGGAGAGTCGCTTCTCATAGGCCACGCCGGGCATCCGGAAGTGGAAGGTACAATGGGCCAGCACAGTAACCCGGCAggggggctgattct                                                                             >  2:239633/2‑242 (MQ=255)
                                                                                                                                                                          cGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCg                                                                        <  1:204726/251‑1 (MQ=255)
                                                                                                                                                                              gTACCGGACGGCGCGATCCTTATTTTCTCCGCACACGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATGTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATc                                                                    <  2:194581/251‑1 (MQ=255)
                                                                                                                                                                               gaCCGGACGGCGAGATTCTGTTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTTGCGTTGTGGCGGTGTTTGATGCCACCGGTCCGCTGGGGTCCAAAGGGCATATGGAAGCCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTCCAGTAACCCGGAAGGGGGAATGTATCTGGTAGaagcg                                                                   <  2:376735/250‑4 (MQ=255)
                                                                                                                                                                                aCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCg                                                                   <  1:17025/250‑1 (MQ=255)
                                                                                                                                                                                  cGGACAGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCg                                                                >  1:489675/1‑251 (MQ=255)
                                                                                                                                                                                    gACGGCGCGATCCGTATTTCTTCCGCACCCGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGAGTCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCTCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCAGGCCGGGCATCCGGAAGTGGAAGGTACGATGTGCCAGTACAGTAGCCCGGAAGGGGGAATTTATCTGGTCGAATTGCCGGa                                                              <  2:640303/251‑1 (MQ=255)
                                                                                                                                                                                     aCGGCGCGATCCTGATTTTCTCCGCTCACGTTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                             <  1:524642/251‑1 (MQ=255)
                                                                                                                                                                                     aCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                             >  1:511134/1‑251 (MQ=255)
                                                                                                                                                                                         cgcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtg                                                         <  1:618366/251‑1 (MQ=255)
                                                                                                                                                                                          gcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtgt                                                        >  1:194208/1‑251 (MQ=255)
                                                                                                                                                                                            gATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGCAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGGATGTATCTGGTCGAATCGCCGGACGAtgtgt                                                      >  1:567684/1‑251 (MQ=255)
                                                                                                                                                                                                   atttctgccGCACACGGTGTTTCTCAGGCGGTCCGTAACGAAGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAG‑CCGTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGGATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAAca                                               <  2:33752/244‑2 (MQ=255)
                                                                                                                                                                                                    atttatgccgactaggttgttCTTAGGCGGGACGTAAAGAAGCAAAAAGCAGCGATTTGACGGGGTTTTATGCAACCGGCCTGCTGGTGCCCAAAGTGCATATGGAAGGCGCCCGCGCCAGCTTGCAG‑GGCGAAGAGTCGATTCCCAACGGACACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCTGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTGGAAACTg                                              <  2:126036/232‑1 (MQ=255)
                                                                                                                                                                                                     tttCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCTCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTg                                              <  1:16671/250‑1 (MQ=255)
                                                                                                                                                                                                           cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTc                                        >  1:589189/1‑250 (MQ=255)
                                                                                                                                                                                                           cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCa                                       <  1:324194/251‑1 (MQ=255)
                                                                                                                                                                                                             cacaCGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGTAGTGGTAGGTACGATGGGCCAGTACAGTAACCCGGGAGGGGGAATGTATCTGGTCGAATCGCCGGACGAGTTGGGGGAACTAACGGTCaa                                      >  2:253061/1‑250 (MQ=255)
                                                                                                                                                                                                                acGGTGGTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGAGTTGACGGTGTGTGAGGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAc                                  <  2:58327/251‑1 (MQ=255)
                                                                                                                                                                                                                    tgtTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGa                                >  1:78523/1‑249 (MQ=255)
                                                                                                                                                                                                                            aGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGACGAGAAGctc                      >  1:483803/1‑251 (MQ=255)
                                                                                                                                                                                                                             ggCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGctct                     <  1:382809/251‑1 (MQ=255)
                                                                                                                                                                                                                               gtgtatgTAACCAAGCAAAAAGCCGTGAGTTGAAGTTGTTTGAGGCACCCTGTCCGCTGGTGCCCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGTAGTTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGGGTGGAAACTGACGGTCAAAATCGAAGAGAAGctctac                   <  2:479171/245‑3 (MQ=255)
                                                                                                                                                                                                                                   accgaaCGAAGCAAAAAGCCGTGATTTGACGGTTTTTGATGCCGCCTTGCCGCTGGTGTCCAAATTGCATTTGGAAGTCGCCCGCGCCAG‑CCGTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                <  2:430311/246‑1 (MQ=255)
                                                                                                                                                                                                                                   aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTa               <  1:146492/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGAc           >  1:397084/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCat                                                                                                                                     >  1:392346/1‑125 (MQ=255)
                                                                                                                                                                                                                                          gAAGCAAAAAGCCGCGATTTGACGGGGTTTGATGCCACCTTTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCTAAGAGTCGATTCTCATCGGCCACGCCGGGCat                                                                                                                                     <  2:392346/126‑2 (MQ=255)
                                                                                                                                                                                                                                                aaaaGCCGCGATTGGACTGTGTTTGTTGCTACCTGGCCGCTGGTGACCAAAGTGCTTATGGAAGTCGCCCGCGCCAG‑CCGGCGTGGCGAAGAGGCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTATCCCGGAAGGGGGAATGTATCTGGTAGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAc  <  2:356409/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aaaGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CCGTCGTGGCGAAGAGTCGATTCTCATCGGCAACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGAcaa   >  1:389067/1‑247 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
TTGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAC  >  NC_000913/26289‑26779

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