Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,669,467 C→T 37.5% R204* (CGA→TGA)  U712_08600 → Flagellar biosynthetic protein fliR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,669,4670CT37.5% 56.8 / 64.3 64R204* (CGA→TGA) U712_08600Flagellar biosynthetic protein fliR
Reads supporting (aligned to +/- strand):  ref base C (13/27);  new base T (9/15);  total (22/42)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

AATTTCGGCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGTTTTCAGTTTAACGATTGAAACGAT  >  CP006881/1669335‑1669603
                                                                                                                                    |                                                                                                                                        
actttCGGCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTg                                                                                                                                    <  2:100320/137‑1 (MQ=255)
 aTTTCGGCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGt                                                                                                                                   >  1:142397/1‑139 (MQ=255)
      ggCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCg                                                                                                                               <  2:155290/138‑1 (MQ=255)
      ggCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCg                                                                                                                               <  7:64708/138‑1 (MQ=255)
      ggCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCg                                                                                                                               <  3:173586/138‑1 (MQ=255)
             aaaaGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTaa                                                                                                                        <  5:105224/138‑1 (MQ=255)
             aaaaGTTTGCCTATTTTATTGCGAAAAGCATTAATGCAAAGTTTATTATCGCTTATCAAATGTCAGCGCCGGTTGTGGCCAGTATGTTTCTACTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTaaa                                                                                                                       <  6:148570/139‑1 (MQ=255)
              aaaGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTaaa                                                                                                                       <  6:290947/138‑1 (MQ=255)
                               ttGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTggtcggtc                                                                                                      <  7:37014/138‑1 (MQ=255)
                                                 cAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGTCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGt                                                                                   <  7:319837/139‑1 (MQ=255)
                                                 cAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGt                                                                                   <  8:58533/139‑1 (MQ=255)
                                                   aTGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGACCTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGACTTTAAAGATCGCTGt                                                                                   >  7:76537/1‑137 (MQ=255)
                                                   aTGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGt                                                                                   >  3:67741/1‑137 (MQ=255)
                                                   aTGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGt                                                                                   >  3:207222/1‑137 (MQ=255)
                                                   aTGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGt                                                                                   >  6:126746/1‑137 (MQ=255)
                                                    tGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTggtcggtc                                                                                                      >  1:143051/1‑117 (MQ=255)
                                                    tGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTggtcggtc                                                                                                      <  2:143051/117‑1 (MQ=255)
                                                          ttctCGCTTTTCAAATGTCAGCGCAGGTTGTGGCCAGTTAGTTTCTAGGTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGGCGTTCTGCCTTTAAAGATCGCTGTCAGTTTCAt                                                                          <  8:44021/136‑1 (MQ=255)
                                                          ttatCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCAt                                                                          <  3:85154/139‑1 (MQ=255)
                                                               gCTTTTCAAATGTCAGCGCCGGTTGTTGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTTTCGTGGCTCGAACTGTTCCGCATTTAAATGTATTTGTGGTCGGTCTGCATTTAAAGATCTGTGTGAGTTTCATTATg                                                                      <  7:175265/138‑1 (MQ=255)
                                                               gCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      <  2:341884/138‑1 (MQ=255)
                                                               gCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      <  6:216500/138‑1 (MQ=255)
                                                               gCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      <  5:126746/138‑1 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    <  5:106652/139‑1 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    >  1:286498/1‑139 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCt                                                                                     <  6:248779/122‑1 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCt                                                                                     >  5:248779/1‑122 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    <  6:112725/139‑1 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    <  4:56546/139‑1 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    >  6:69195/1‑139 (MQ=255)
                                                                cTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    <  5:69195/139‑1 (MQ=255)
                                                                   ttCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGaa                                                                                                                                       >  7:215233/1‑69 (MQ=255)
                                                                   ttCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGaa                                                                                                                                       <  8:215233/69‑1 (MQ=255)
                                                                   ttCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  <  2:286498/138‑1 (MQ=255)
                                                                   ttCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  <  4:207222/138‑1 (MQ=255)
                                                                   ttCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  <  2:335334/138‑1 (MQ=255)
                                                                   ttCAAATGTCAGCGCCGGTTGTGGACAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  <  6:507/138‑1 (MQ=255)
                                                                            cAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTgtatgtat                                                        <  4:38814/139‑1 (MQ=255)
                                                                            cAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTgtatgtat                                                        <  3:354326/139‑1 (MQ=255)
                                                                                    gTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTa                                                >  8:215555/1‑139 (MQ=255)
                                                                                    gTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGATGTCAGTTTCATTATGCTCa                                                                  <  4:211558/121‑1 (MQ=255)
                                                                                    gTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGATGTCAGTTTCATTATGCTCa                                                                  >  3:211558/1‑121 (MQ=255)
                                                                                        tgGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCtt                                            <  6:87428/139‑1 (MQ=255)
                                                                                         ggCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTtc                                           >  7:178434/1‑139 (MQ=255)
                                                                                          gCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTtcg                                          >  1:246436/1‑139 (MQ=255)
                                                                                          gCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTtcg                                          >  8:248178/1‑139 (MQ=255)
                                                                                              gtttgtttCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATAGCTGTCAGTTTCATTATGCTCATTCTATGTATCTCCGTTATCTTcgtcgt                                      >  5:22980/1‑139 (MQ=255)
                                                                                                   tttCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGGTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTcc                                 <  6:197386/139‑1 (MQ=255)
                                                                                                    ttCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      >  4:199170/1‑101 (MQ=255)
                                                                                                    ttCTAGTTGATTTAGCATTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      <  3:199170/101‑1 (MQ=255)
                                                                                                         gTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGa                              <  5:171646/136‑1 (MQ=255)
                                                                                                             aTTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGttt                        <  7:215555/138‑1 (MQ=255)
                                                                                                             aTTTAGCATTAGGTAGCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGtttt                       <  2:132237/139‑1 (MQ=255)
                                                                                                                    aTTAGGTATCGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGGCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGTTTTCAGTTTa                <  5:320767/139‑1 (MQ=255)
                                                                                                                       aGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      <  1:112017/82‑1 (MQ=255)
                                                                                                                       aGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATg                                                                      >  2:112017/1‑82 (MQ=255)
                                                                                                                             cGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGTTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  >  4:81624/1‑80 (MQ=255)
                                                                                                                             cGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    <  2:310713/78‑1 (MQ=255)
                                                                                                                             cGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCt                                                                    >  1:310713/1‑78 (MQ=255)
                                                                                                                             cGTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  <  3:81624/80‑1 (MQ=255)
                                                                                                                             cGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  >  1:221259/1‑80 (MQ=255)
                                                                                                                             cGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCa                                                                  <  2:221259/80‑1 (MQ=255)
                                                                                                                              gTGGCTTGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGTTTTCAGTTTAACGATTGaaa      <  5:196389/139‑1 (MQ=255)
                                                                                                                                  cTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGTTTTCAGTTTAACGATTGAAACGAt  >  1:255149/1‑139 (MQ=255)
                                                                                                                                    |                                                                                                                                        
AATTTCGGCGATGAAAAGTTTGCCTATTTTATTGCGAAAAGCTTTAATGCAATGTTTATTATCGCTTTTCAAATGTCAGCGCCGGTTGTGGCCAGTTTGTTTCTAGTTGATTTAGCATTAGGTATCGTGGCTCGAACTGTTCCGCAGTTAAATGTATTTGTGGTCGGTCTGCCTTTAAAGATCGCTGTCAGTTTCATTATGCTCATTGTATGTATGTCCGTTATCTTCGTCGTTGTCCGGAATGTTTTCAGTTTAACGATTGAAACGAT  >  CP006881/1669335‑1669603

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