Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 914,631 G→C 7.7% A448P (GCC→CCC) ‡ fruA → PTS fructose transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947914,6310GC 7.7% 155.7 / 3.9 52A448P (GCC→CCC) ‡fruAPTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base C (2/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-02

ATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAA  >  NC_002947/914484‑914779
                                                                                                                                                     |                                                                                                                                                      
aTGGGCACCACCAATGCCATCCTGCTGGGGGTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGGGGGCGGATCTAGAAGGCGGGGTATGCGTTGTCGGTGGTGCTGCTGGCGTAGTCGTGCTAAGCGCCAATGGAGGAGACGATGgc                                                                                                                                                       >  8:283308/1‑149 (MQ=255)
      accaccAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCgg                                                                                                                                                 >  6:143105/1‑149 (MQ=255)
          ccaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATg                                                                                                                                            <  5:194048/149‑1 (MQ=255)
          ccaATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATg                                                                                                                                            <  2:76434/149‑1 (MQ=255)
             aTGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTg                                                                                                                                        <  5:143105/149‑1 (MQ=255)
               gCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGCCGGCGGCCATGGTCGCCGCCAt                                                                                                                                             >  8:286728/1‑143 (MQ=255)
               gCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGCCGGCGGCCATGGTCGCCGCCAt                                                                                                                                             <  7:286728/143‑1 (MQ=255)
                   accTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTAGGCGTTCTCGGTGGGGCTGCTGGCCGCGTCGAGCTACGCGCCAAACGGCGGCACCATGCCGGCGGCCATGGTCgccgcc                                                                                                                                   <  1:234079/148‑1 (MQ=37)
                      tgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGAtt                                                                                                                               >  1:241373/1‑149 (MQ=255)
                        ctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTgg                                                                                                                             >  2:163969/1‑149 (MQ=255)
                        ctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCACCATGCCGGCGGCCATGGTCgccgccattg                                                                                                                               <  5:211059/147‑5 (MQ=38)
                           gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCtt                                                                                                                          >  8:92499/1‑149 (MQ=255)
                           gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCtt                                                                                                                          <  1:163969/149‑1 (MQ=255)
                           gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCtt                                                                                                                          <  6:195889/149‑1 (MQ=255)
                           gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCtt                                                                                                                          >  5:193922/1‑149 (MQ=255)
                             gCTGTTGCTGGGCGGCATGATGTGCGTAGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTAGGTTGGGCTGCAGGGCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCCGCAAGGGGCCGGCGAttcggatgg                                                                                                                        >  8:333047/1‑142 (MQ=255)
                                 ttGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTAt                                                                                                                    >  5:305315/1‑149 (MQ=255)
                                               gatgTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTgg                                                                                                      >  4:191459/1‑149 (MQ=255)
                                                     cGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGcc                                                                                                                >  6:227139/1‑133 (MQ=255)
                                                     cGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGcc                                                                                                                <  5:227139/133‑1 (MQ=255)
                                                          aCCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGTGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAg                                                                                           >  2:226178/1‑149 (MQ=255)
                                                            cTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGtt                                                                                         >  7:71543/1‑149 (MQ=255)
                                                                      cGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                               >  8:296304/1‑149 (MQ=255)
                                                                      cGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagc                                                                               >  4:151684/1‑149 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                         <  3:191459/149‑1 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                         <  7:92499/149‑1 (MQ=255)
                                                                            acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCCAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                         >  1:42111/1‑149 (MQ=255)
                                                                                   cgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                                  <  4:214909/149‑1 (MQ=255)
                                                                                      cTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                >  1:282831/1‑148 (MQ=255)
                                                                                                 gggTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                                  <  6:313021/134‑1 (MQ=255)
                                                                                                  cgTGGGGCGGCTGGCCGCGGCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATCGTGCCGCCGATTGGCGTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctgg                                                   <  5:345637/148‑1 (MQ=255)
                                                                                                        gctgctGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                             >  5:103108/1‑149 (MQ=255)
                                                                                                         ctgctgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                            >  7:193104/1‑149 (MQ=255)
                                                                                                         ctgctgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTggg                                            <  3:169858/149‑1 (MQ=255)
                                                                                                            ctgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTgggg                                           <  6:305315/147‑1 (MQ=255)
                                                                                                            ctgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCt                                         >  4:139962/1‑149 (MQ=255)
                                                                                                            ctgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCt                                         >  5:278485/1‑149 (MQ=255)
                                                                                                                ccTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                                     >  1:46850/1‑149 (MQ=255)
                                                                                                                ccTCGTCGAGCTACGCGCCAATGGCGGCGACCAGGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                                     <  2:282831/149‑1 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAATGGCTGAGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGTGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCa                                 >  8:162960/1‑149 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                <  2:24594/118‑1 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                >  1:24594/1‑118 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                  <  7:296304/148‑1 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                  <  1:241025/148‑1 (MQ=255)
                                                                                                                    gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                  >  2:241025/1‑148 (MQ=255)
                                                                                                                          cttcgcgcgATTGGCTGGGACCATGGCGGCGTGCATGGTGCCGCCGATTGGCTTGGGTATTGCCTCCTTCCTGGCGCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAtctc                             <  7:283308/144‑1 (MQ=255)
                                                                                                                            ccaagcCAATCGGCGGCACCATGCCGGCGGCCATGGTCGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGCGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTcc                            >  3:257347/5‑147 (MQ=25)
                                                                                                                                gcCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGTGTATTGCCACCTTCCTGGCCCGTCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAggggg                     <  2:241373/149‑1 (MQ=255)
                                                                                                                                  cAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGc                    >  2:319086/1‑148 (MQ=255)
                                                                                                                                           ggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGt            <  2:42111/149‑1 (MQ=255)
                                                                                                                                            ggcaCCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt           >  5:132805/4‑149 (MQ=38)
                                                                                                                                            gcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt           <  3:173678/149‑1 (MQ=255)
                                                                                                                                                  catgccggcGGCCATGGTCGCCGCCAATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCggg                                                                 <  6:191596/88‑1 (MQ=37)
                                                                                                                                                  aTGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCGCCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgc     >  1:248257/1‑149 (MQ=255)
                                                                                                                                                     gccgccGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCaa  <  3:139962/149‑1 (MQ=255)
                                                                                                                                                     gccgccGGCATGGTGCCGCCGATTGGCTGGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCaa  <  6:193922/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAA  >  NC_002947/914484‑914779

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: