Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 690,948 G→A V279V (GTC→GTT)  ybeX ← CorC‑HlyC family protein YbeX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913690,9480GA94.1% 104.4 / ‑5.1 35V279V (GTC→GTT) ybeXCorC‑HlyC family protein YbeX
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (20/12);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCTTCTTC  >  NC_000913/690702‑691193
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
gCCAATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACAt                                                                                                                                                                                                                                                   >  1:91179/1‑251 (MQ=255)
 caaaaagaggCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCCGACGCCCCATCAGCGGAATAATCGCCGCCGGCCAGACGACGTAAGGAGAGAAAGCCAGCGCTCCGCAGGCACAGAATAATAACCCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGGAGGTATCTATCCAGTTTCGGTTTTATTCATCCAGCGTAGGCTGGGGTGAGTCATCCGGGATTTTAACATg                                                                                                                                                                                                                                                  >  2:256681/9‑251 (MQ=255)
            aGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGaata                                                                                                                                                                                                                                        >  1:237014/1‑250 (MQ=255)
                  gAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGcc                                                                                                                                                                                                                                 <  2:110387/251‑1 (MQ=255)
                                            ccTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAAc                                                                                                                                                                                                       <  2:73780/251‑1 (MQ=255)
                                              tGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTg                                                                                                                                                                                                     >  1:132396/1‑251 (MQ=255)
                                                             gAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTGTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACGTTGAACTGGTAACCggcgatgtc                                                                                                                                                                                      <  1:367558/251‑1 (MQ=255)
                                                                  aatCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTGTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATa                                                                                                                                                                                  <  1:258995/250‑1 (MQ=255)
                                                                                aGACGTAGTAAGGAGAGCAGGACAGCAATCCGCCGGCCCCGCATAATAACGACCGCAGCAGGCGAACGCGCTGGCGTTCAATTCATGAGGCCAAACCCATGTAGTTATCTATCCAGTTTAGGTTTTATTCATCCAACTTCGGCTGGGGTGACTCATCCGGGCTTTTCACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCCTCGATGTCGCTAGTTTCGCAACGCGcc                                                                                                                                                                   >  2:183224/1‑251 (MQ=255)
                                                                                 gACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCg                                                                                                                                                                  >  1:15789/1‑251 (MQ=255)
                                                                                     tagtAAGGAGAGCAGGCCAGCGTTCAGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGACTAATACGCCGACTGTCGGCCATCGCCACTTTGCACTGGTAACCGACGATGTCGATAGTTTCGCCACGCGCCGGa                                                                                                                                                                >  2:6580/3‑249 (MQ=255)
                                                                                       ggaaGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGAGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACGGTCGTTCATCGCCTCTTTGAACTGGTATCCGTCGATGTCGATAGGTTCGCCACGCGCCGGAAGAt                                                                                                                                                            <  2:247882/249‑1 (MQ=255)
                                                                                                             ccGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCAcaag                                                                                                                                      <  1:76216/251‑4 (MQ=255)
                                                                                                                                   ccAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGActtct                                                                                                                <  1:134566/251‑1 (MQ=255)
                                                                                                                                                      cTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGc                                                                                             >  1:40142/1‑251 (MQ=255)
                                                                                                                                                        tgCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATTCCTGCATCACCAGACCACCGATAGTGTCGACCTCTGCATCGCTATAGTGGGTgacg                                                                                           >  2:25800/2‑248 (MQ=255)
                                                                                                                                                              cAATTAATGAGGCAAAAGACATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCACGCTGCGGCTGGGGTGAGTCCTCCAGGATTTTAACATGATCCTGCATAATACGCCGACTGTCAGACATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCCCGAGCCGAAAGATGCACAAATGCACGCATCACCAGCCCACAGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGACGACCGCt                                                                                     >  2:376129/1‑251 (MQ=255)
                                                                                                                                                                aTTAATAAGGCAAATGCCATGTAGTTATCTAGCCAGTTGCGGTTTTATTCATCCAGCGTCGGCTGGGGGGAGGCATCCGGGATTTTAACATGAACCTCAATAATACGCCGACTGGCGGCCATCGCCACTTTGAACTGGGAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGTATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGTACGCTTc                                                                                   <  2:433073/251‑1 (MQ=255)
                                                                                                                                                                  tAATGAGGCAAAAGCCATGTAGTTATCTATCAAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAAACTGACTAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATCGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                                                                 >  2:31869/1‑251 (MQ=255)
                                                                                                                                                                     tGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCAAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGCACGCTTCGTTgc                                                                              >  2:108508/1‑250 (MQ=255)
                                                                                                                                                                           aaaaGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCtt                                                                        >  1:335401/1‑251 (MQ=255)
                                                                                                                                                                              agccatggaggaatctatcCAGTTTCGGTGTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCTCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATAGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCa                                                                      >  2:96873/13‑250 (MQ=255)
                                                                                                                                                                                 cATGTAGTTATATATCCAGATTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACAACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATgg                                                                  >  2:208407/1‑251 (MQ=255)
                                                                                                                                                                                      aGTTATCTATCCAGTTTCGGTATTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGATCCTGAATAATACGCCGACTGTCGGGCAGCGCCACTTTGAAGTGGTAACCGTCGATGGCGATAGTTTCGCGACGCGTCGGAAGATTCCCAGATGCCGGCAGCACCAGATCACCGATAGTGTCGACTTCTTCATCGCTAAAGGGGGTTCCGAACGCTTCGTTGAAGTCTTCAATGGAAGc                                                              <  1:193886/250‑1 (MQ=255)
                                                                                                                                                                                            ctatcCAGGTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGCACATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAATTGGGTGCCGAACGCTTCGTTGAAGACTTCAATGTGATCCAGTgcg                                                       >  2:274439/1‑251 (MQ=255)
                                                                                                                                                                                               tcCAGTTTCGGTTTTATTCATCCCGCTTCGACTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGACTAATACGCCGCCTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCtt                                                                        >  1:381806/1‑231 (MQ=255)
                                                                                                                                                                                               tcCAGTTTCGGTTGTATTCAGCCAGCGTCGGCTGGGGTGAGTCATCCGGGATTTGGACATGAGCCTGAATAATACGCCGACTGTCGGCCATCGCCTCTGTGAACTGGTAACCGGCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCAGCACCAGACCACCGATAGTGTCGACTTCTGCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCtg                                                                        <  2:381806/231‑2 (MQ=255)
                                                                                                                                                                                                ccAGTGTCGGTTTGATTCCTCCAGCTTCGGCTGAGGTGAGTCAGCCGGGATTTTAACAGGAACATGAAAAATACGCCGACTGTCGGCCAGCGCCACTTTGAACAGGTAACCGTCGAGGTCGATAGTTTCGCCACGCGCCGGAAGAGGCCCAAATGCCTGCATCACCAGACCACCGATATTGTCGACTTCTTCATCGCTAAAGTGGGTGCAGAACGCTGCGTTAAAGTCTTCAATGGAAGCCAGTGCGCGTa                                                   <  2:13371/251‑1 (MQ=255)
                                                                                                                                                                                                              aTTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATg                                      >  1:24917/1‑250 (MQ=255)
                                                                                                                                                                                                                             gCTGGGGTGAGTCATCCGGGATTTTAACATGAACCTGCATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCTTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCACTGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACgg                      >  1:50640/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gTCATCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGatat            >  1:54261/1‑251 (MQ=255)
                                                                                                                                                                                                                                          aTCCGGGATTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCCTCACCAGACCACCGATAGTGTCGACTTCTTCCTCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTCCGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATatcat         >  2:422637/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aTTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCCTCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttattc  >  2:149993/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aTTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttcttc  <  1:31869/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 aTTTTAACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCttcttc  <  1:229315/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     
GCCAATAGCGGCAGACTGGAGTGGACGGCGGTTAAAGGTCAACGCCTGAAGCCCCATCAGCGAAATAATCGCCGCAGGCCAGACGTCGTAAGGAGAGAAGGCCAGCGTTCCGCAGGCACCGAATAATAACGCCAGCAGCAGGCGAATGCGCTGGCGTTCAATTAATGAGGCAAAAGCCATGTAGTTATCTATCCAGTTTCGGTTTTATTCATCCAGCTTCGGCTGGGGTGAGTCATCCGGGATTTTGACATGAACCTGAATAATACGCCGACTGTCGGCCATCGCCACTTTGAACTGGTAACCGTCGATGTCGATAGTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGGTCCAGGTATGACGACTCAGCTGACGGAAGTCGATATCATCTTCTTC  >  NC_000913/690702‑691193

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: ‑