Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 10,718 G→A S213S (AGC→AGT yaaW ← putative enzyme‑specific chaperone YaaW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091310,7180GA100.0% 232.3 / NA 74S213S (AGC→AGTyaaWputative enzyme‑specific chaperone YaaW
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (42/32);  total (42/32)

ATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAT  >  NC_000913/10471‑10959
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
aTTAGCCAACTTAGTGTTGCCCATAATTCATCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTg                                                                                                                                                                                                                                                <  1:70220/251‑1 (MQ=255)
    gCCAACTTAGTGTTGCCCATAATTCCTCAAACATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGc                                                                                                                                                                                                                                            >  1:162223/1‑251 (MQ=255)
       aaCTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTg                                                                                                                                                                                                                                          >  2:8855/1‑250 (MQ=255)
        aCTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGAc                                                                                                                                                                                                                                        <  1:375541/251‑1 (MQ=255)
        aCTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGAc                                                                                                                                                                                                                                        >  1:262805/1‑251 (MQ=255)
         cTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGAcc                                                                                                                                                                                                                                       >  2:399379/1‑251 (MQ=255)
          ttAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCg                                                                                                                                                                                                                                      >  1:235199/1‑251 (MQ=255)
           tAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTGTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCg                                                                                                                                                                                                                                      <  1:558649/250‑1 (MQ=255)
             gtgtTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCtt                                                                                                                                                                                                                                   <  1:576581/251‑1 (MQ=255)
                 tGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAAc                                                                                                                                                                                                                               >  1:444425/1‑251 (MQ=255)
                        aaTTCCTCAAACATCATCATCGAATGAATGGTGAACTAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCAttt                                                                                                                                                                                                                        >  2:265862/1‑251 (MQ=255)
                         aTTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTa                                                                                                                                                                                                                       <  1:590593/251‑1 (MQ=255)
                          ttCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTa                                                                                                                                                                                                                       >  1:206460/1‑250 (MQ=255)
                              tCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGct                                                                                                                                                                                                                  >  2:331567/1‑250 (MQ=255)
                                 aaaaTCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAAATGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTaatgatgcg                                                                                                                                                                                                               >  2:445946/1‑243 (MQ=255)
                                             gaatgaatGGTGAAATAATTGCCCTGAATAACAGTAGTGTTGTCAGTGCGCGGCAGAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCTATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCt                                                                                                                                                                                                   <  2:218394/251‑1 (MQ=255)
                                                  aatgGTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCGGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAg                                                                                                                                                                                              <  1:265815/251‑1 (MQ=255)
                                                     ggTGAAATAATTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACAGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAgacc                                                                                                                                                                                           >  2:39454/1‑248 (MQ=255)
                                                           ataagTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACccgcgcc                                                                                                                                                                                     <  1:404159/251‑1 (MQ=255)
                                                             aatTTCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACccgcgccgc                                                                                                                                                                                   >  2:81798/1‑251 (MQ=255)
                                                                 tCCCTGAATAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCgca                                                                                                                                                                                <  1:412074/250‑1 (MQ=255)
                                                                    ctgaataacagcagcgtactCAGCGCGTGGCATAATAATCAGCCAGTGGGGCAGTGTATACGATCTTTTAAGGAGAAAATGAAAATTTTCCCCGGTTTCCGGAGTCAGACCTGAGCGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAAACCCATTTAATGCTGCGCCCACCGGGCCTCCCAGACCCGCGCCGCgcagca                                                                                                                                                                             >  2:567459/20‑250 (MQ=255)
                                                                       aaTAACAGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAAcc                                                                                                                                                                         >  2:192709/1‑251 (MQ=255)
                                                                          aacagTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAACATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACAGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATg                                                                                                                                                                      >  2:617166/5‑251 (MQ=255)
                                                                             agTAGAGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGccc                                                                                                                                                                   >  2:215432/2‑251 (MQ=255)
                                                                                       tCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTc                                                                                                                                                         >  1:213127/1‑251 (MQ=255)
                                                                                          gggCGCGGCATAATAATCAGGCAGTTGGGCAGTGTCTACGATCTGTTGAGGGGAAAATGAAAATTTGCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCagt                                                                                                                                                      <  2:470107/251‑3 (MQ=255)
                                                                                          gggCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCAtt                                                                                                                                                      <  2:543496/251‑1 (MQ=255)
                                                                                          aggCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGCCCCGCCCCGCGCAGCAAACCAAGCCCAAGCACGCTCAtt                                                                                                                                                      >  2:342499/2‑251 (MQ=255)
                                                                                               cgGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGGCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATtgctgc                                                                                                                                                 <  2:98841/251‑1 (MQ=255)
                                                                                                    taataaTCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTggg                                                                                                                                            >  1:216038/1‑251 (MQ=255)
                                                                                                     aataatCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATGTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTggg                                                                                                                                            <  1:625303/250‑1 (MQ=255)
                                                                                                        aatCAGCCAGTGGGGCAGTGTCTCCGATCTTTTGAGAGGAAAATGAAAATTTTCACCGGTTTACGGTGTCAGACCTGAGTGGAGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCAGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCg                                                                                                                                        >  2:75065/1‑251 (MQ=255)
                                                                                                         atCAGCCAGTGGGGCAGTGTCTGCGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCAGTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCCCGCTCATTGCTGCGTGAGTGCGt                                                                                                                                       <  2:422526/251‑1 (MQ=255)
                                                                                                               cAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATg                                                                                                                                 >  2:17919/1‑251 (MQ=255)
                                                                                                                   ggggCAGTGTCTACGATCTTTTGAGGGGAACATGAAAATTTTCCCCGGGTTACGGTGGCAGACCTGAGGGGCGCTAAACATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCAAGCGATAGGCACTGCCGCTGCCCGCTTTCACCCCATTTAATGCTGCGCCCACAGAGACTCCCAGCCCCGCTCCGCGCCGCAAACCATGCCCCAACACGCACATTGCTGCGTGGGTGCGGAAAATGAggg                                                                                                                             >  2:335582/1‑251 (MQ=255)
                                                                                                                       cAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCAGCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCGCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCAGTGCTGCGTGGGTGCGTAAAATGCGGGTCaa                                                                                                                         <  2:125052/251‑2 (MQ=255)
                                                                                                                       cAGTGTCTACGATCTTTTGAGAGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGCTAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCCCAGCAAACCATGCCCCAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCac                                                                                                                         >  2:130943/1‑250 (MQ=255)
                                                                                                                           gtCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTg                                                                                                                      <  1:200133/250‑1 (MQ=255)
                                                                                                                               aCGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTGTCGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCAGTGAATGCTGCGCCCACAGGGCCTCCCAGACGCGCGCCGCGCAGCAAACCATGCCCAAGCACGTTCAGTGCTGCGTGGGTGCGTAAAATGCTGGTCAATTGGCt                                                                                                                   <  2:53067/249‑1 (MQ=255)
                                                                                                                               aCGATCTTTTGAGGGGAAAATGAAAATTGTCCCCGGTTTACGGTGTCATACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAATCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTgg                                                                                                                 <  1:574900/251‑1 (MQ=255)
                                                                                                                                cGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTgg                                                                                                                 >  1:221403/1‑250 (MQ=255)
                                                                                                                                 gATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGGGGGTGCGTAAAATGCGGGTCAATTGGCTGGa                                                                                                                >  2:442755/1‑250 (MQ=255)
                                                                                                                                  aTCTTTTGAGGGGAAAATGAACATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGACGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGCTGCGTAAAATGCGGGTCAATTGGCTGGaa                                                                                                               >  2:428285/1‑250 (MQ=255)
                                                                                                                                   tCTTTTGCGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGaa                                                                                                               >  2:97259/1‑249 (MQ=255)
                                                                                                                                      tttGAGGGGAAAATGAACATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCa                                                                                                           >  2:90784/1‑250 (MQ=255)
                                                                                                                                        tGAGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAAt                                                                                                        >  1:194678/1‑251 (MQ=255)
                                                                                                                                          aggtgaaaattaaaatTTTCCCCGGTTTCCGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATACGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGGCGAGCAGCAAACCATGCCCAAGGACGCTCATGGCGGCGTTGGTGCGTAAAATGCTGTTCAATGGGCTGTAAAGCAAATgt                                                                                                      <  2:630034/247‑2 (MQ=255)
                                                                                                                                          aGGGGAAAATGAAAATTTTCCCCGGTTTACGGTGGCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGc                                                                                                      <  1:304951/251‑1 (MQ=255)
                                                                                                                                             ggAAAATGAAAATTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGac                                                                                                   <  1:39454/251‑1 (MQ=255)
                                                                                                                                                      aaatTTTCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCtttt                                                                                          >  1:164920/1‑251 (MQ=255)
                                                                                                                                                         ttttCCACGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGcc                                                                                       >  1:157498/1‑251 (MQ=255)
                                                                                                                                                          tttCCCCGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGcc                                                                                       >  2:544372/1‑250 (MQ=255)
                                                                                                                                                             ccccGGTTTACGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTACTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCATTTGACaat                                                                                    >  2:6411/1‑250 (MQ=255)
                                                                                                                                                                 ggtgtacGGTGTCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAAttt                                                                               <  2:427543/245‑1 (MQ=255)
                                                                                                                                                                         gtgtCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTtca                                                                       <  2:129254/248‑1 (MQ=255)
                                                                                                                                                                         ggggcAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTtca                                                                       <  2:206460/247‑1 (MQ=255)
                                                                                                                                                                          tgtCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTCCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATACTTTGTCTTtcat                                                                      >  1:231692/3‑251 (MQ=255)
                                                                                                                                                                           gtCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATGTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTtcat                                                                      <  2:444425/249‑1 (MQ=255)
                                                                                                                                                                                    gAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGgc                                                            <  2:216038/251‑1 (MQ=255)
                                                                                                                                                                                      gTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGgcag                                                          <  1:498833/251‑1 (MQ=255)
                                                                                                                                                                                         gcgcTAACAATCCGGCGCAGGCAGGCGATTTGAAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTCACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGgcagc                                                         >  2:328003/1‑249 (MQ=255)
                                                                                                                                                                                          cgcTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCagcagc                                                      >  1:528753/1‑251 (MQ=255)
                                                                                                                                                                                                 cATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTccc                                               >  1:255042/1‑250 (MQ=255)
                                                                                                                                                                                                        cgcAGGCAGGCGCTTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGACGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCAAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCAtt                                        >  1:555653/1‑251 (MQ=255)
                                                                                                                                                                                                         ccaggcaggCGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTc                                       <  1:6411/250‑1 (MQ=255)
                                                                                                                                                                                                            ggcaggcGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCCGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCAcc                                    >  2:307067/1‑251 (MQ=255)
                                                                                                                                                                                                             gcaggcGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCCTTCAccc                                   >  2:538095/1‑251 (MQ=255)
                                                                                                                                                                                                              caggcGATTTGCCGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCt                                  >  1:143180/1‑251 (MQ=255)
                                                                                                                                                                                                               aggcGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTg                                 <  2:213127/251‑1 (MQ=255)
                                                                                                                                                                                                                 gcGATTTGCAGTACTGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGTCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGc                               <  2:450129/251‑1 (MQ=255)
                                                                                                                                                                                                                  cGATTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCa                              <  1:235973/251‑1 (MQ=255)
                                                                                                                                                                                                                    aTTTGCAGTACGGCTGGAATTGTCACGCGATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGAGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATc                            >  2:274130/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gcgATAGGCACTGCCGCTGACCGCTTTAACCCCATTTAATGCTGCGCCCACAGGGCCTCCCAGACCCGCGCCGCGCAGCAAACCATGCCCAAGCACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAATTGGCTGGAAAGCAAATGCGACACGCCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAt  <  2:142701/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
ATTAGCCAACTTAGTGTTGCCCATAATTCCTCAAAAATCATCATCGAATGAATGGTGAAATAATTTCCCTGAATAACTGTAGTGTTTTCAGGGCGCGGCATAATAATCAGCCAGTGGGGCAGTGTCTACGATCTTTTGAGGGGAAAATGAAAATTTTCCCCGGTTTCCGGTATCAGACCTGAGTGGCGCTAACCATCCGGCGCAGGCAGGCGATTTGCAGTACGGCTGGAATCGTCACGCGATAGGCGCTGCCGCTGACCGCTTTAACCCCATTTAGTGCCGCACCTACAGGGCCTCCCAGCCCCGCGCCGCGCAGCAAACCATGCCCAAGTACGCTCATTGCTGCGTGGGTGCGTAAAATGCGGGTCAGTTGGCTGGAAAGCAAATGCGACACACCTTTTGCCAATAATTTGTCTTTCATCAGCAGCGGCAGCAGCTCTTCCAGCTCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTAT  >  NC_000913/10471‑10959

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