Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,551,327 A→C 17.0% intergenic (+44/+82) groL → / ← pykA chaperonin GroEL/pyruvate kinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,551,3270AC17.0% 106.6 / 22.5 59intergenic (+44/+82)groL/pykAchaperonin GroEL/pyruvate kinase II
Reads supporting (aligned to +/- strand):  ref base A (26/23);  new base C (5/5);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

ATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGT  >  NC_002947/1551185‑1551468
                                                                                                                                              |                                                                                                                                             
atCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgg                                                                                                                                        <  2:290904/150‑1 (MQ=255)
 tCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCggg                                                                                                                                       <  2:1215939/150‑1 (MQ=255)
  cagcgccGAAGCCATGGTTTGCGACCTGCCGGAAGCCGAGCCAGCTGCTGGCATGCCTGGTATTGGTGGCATGGGTGGCCTGCCCGGCATGATGTAAGCCAGCCTTGCACTCCTGCATCATCAAAAAGCCCCGGTTCTCGCGCACCgggg                                                                                                                                      <  2:1241830/145‑1 (MQ=255)
  caccacCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCgggg                                                                                                                                      <  1:416957/150‑1 (MQ=255)
  caccacCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATAAAAAAGCCCCGGTTCTCGCGAACCgggg                                                                                                                                      <  2:567322/150‑1 (MQ=255)
          aaGCCATGGTTGCCGGCCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTa                                                                                                                              <  2:306339/150‑1 (MQ=255)
            gCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTAtt                                                                                                                            <  2:691934/150‑1 (MQ=255)
                 ggTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCga                                                                                                                                               >  1:1105747/1‑126 (MQ=255)
                 ggTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCga                                                                                                                                               <  2:1105747/126‑1 (MQ=255)
                      ccGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCgggg                                                                                                                  <  2:608114/150‑1 (MQ=255)
                                ggAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCgca                                                                                                        >  1:182667/1‑150 (MQ=255)
                                            aGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCGACAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTgc                                                                                            <  1:487421/150‑1 (MQ=255)
                                               tgctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCt                                                                                         <  2:182667/150‑1 (MQ=255)
                                                 ctgGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTg                                                                                       >  1:712186/1‑150 (MQ=255)
                                                         ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtgt                                                                               <  1:309387/150‑1 (MQ=255)
                                                         ccTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTtgt                                                                               >  1:165228/1‑150 (MQ=255)
                                                             atatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCGGGGc                                                                                                                                     >  2:1181291/1‑92 (MQ=255)
                                                             atatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCTCGCGAACCGGGGc                                                                                                                                     <  1:1181291/92‑1 (MQ=255)
                                                              tatGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCtgt                                                                          <  2:712186/150‑1 (MQ=255)
                                                                 gggtggcatgggtggcatGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGaa                                                                       >  1:558192/1‑150 (MQ=255)
                                                                      gcatgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAg                                                                  >  2:42954/1‑150 (MQ=255)
                                                                         tgggtggcatgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTgg                                                               >  1:552909/1‑150 (MQ=255)
                                                                                 atgggAGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCACCCGGTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTTAGCCTTGTGCTGGGAAATGAGTGGATCAACCg                                                       >  2:500547/1‑150 (MQ=255)
                                                                                  tgggCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGa                                                      >  1:3116/1‑150 (MQ=255)
                                                                                            gatgTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATc                                            >  1:78771/1‑150 (MQ=255)
                                                                                               gTAAGCCAGCCTTCCCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCAcc                                         <  2:442005/150‑1 (MQ=255)
                                                                                                tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGTCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAAATAGTGGTGCAACCGACCAGCGGATCACCg                                        >  1:914053/1‑150 (MQ=255)
                                                                                                tAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAACCTTTTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCg                                        >  1:914143/1‑150 (MQ=255)
                                                                                                 aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                             >  1:464384/1‑144 (MQ=255)
                                                                                                 aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                             >  1:463874/1‑144 (MQ=255)
                                                                                                 aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                             <  2:464384/144‑1 (MQ=255)
                                                                                                 aaGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGAt                                             <  2:463874/144‑1 (MQ=255)
                                                                                                  atccagccTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGAc                                      <  1:500547/148‑1 (MQ=255)
                                                                                                    ccagccTTACCCCCCTGCACCATAAAGACATGCTGGCCAGCGATGACCGTGGCGTCTTTATTGCTGTCGGGGCCTTTCCGAATTCCCGGTCGGCGCCGGAGCCTTGTGCTGTGAAACGAGTGGATCAACCCGCCAACGGAGCACCGACGt                                    >  1:1254847/1‑150 (MQ=255)
                                                                                                          ttACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGa                              >  1:821904/1‑150 (MQ=255)
                                                                                                           tACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTgcgc                                                                                          >  2:953183/1‑89 (MQ=255)
                                                                                                           tACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTgcgc                                                                                          <  1:953183/89‑1 (MQ=255)
                                                                                                            aCCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGc                                                                                                                 <  2:864906/65‑1 (MQ=255)
                                                                                                            aCCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGc                                                                                                                 >  1:864906/1‑65 (MQ=255)
                                                                                                              cccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCtt                          >  2:165237/1‑150 (MQ=255)
                                                                                                               ccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGc                                                                                                                                     >  2:1277300/1‑42 (MQ=255)
                                                                                                               ccccTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGc                                                                                                                                     <  1:1277300/42‑1 (MQ=255)
                                                                                                                   tGCACCATCAAAAAGCCCCGGTTCGCGAGAACGGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGc                     <  2:558192/150‑1 (MQ=255)
                                                                                                                        cATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTg                <  2:1254847/150‑1 (MQ=255)
                                                                                                                         aTCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTgg               >  2:548318/1‑150 (MQ=255)
                                                                                                                              aaaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATGGCTGTCGGGGCCTCTTCGCAGCACAAGGCGGCGCTCGAGCCTTGGGCAGTGAAACGAGTGGATTAACCGACCATCGGATCACCGACGTGCATGAACTTCATACCATTGGTGCCa          >  2:902775/1‑150 (MQ=255)
                                                                                                                               aaGCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc         >  2:862993/1‑150 (MQ=255)
                                                                                                                               aaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc         <  2:78771/150‑1 (MQ=255)
                                                                                                                               aaGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCAc         <  1:42954/150‑1 (MQ=255)
                                                                                                                                 gCCCCGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATGAACCGACCAGCGGATCAGCGACGTGCCGGATCTTCATGCCCTTGGTGCCACCg       >  1:58685/1‑150 (MQ=255)
                                                                                                                                  ccccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                              >  1:1226837/1‑94 (MQ=255)
                                                                                                                                  ccccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGa                                                              <  2:1226837/94‑1 (MQ=255)
                                                                                                                                  ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGc                                  <  1:889013/122‑1 (MQ=255)
                                                                                                                                  ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGc                                  >  2:889013/1‑122 (MQ=255)
                                                                                                                                  ccccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGa      >  1:236327/1‑150 (MQ=255)
                                                                                                                                   cccGGTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGAt     >  2:321531/1‑150 (MQ=255)
                                                                                                                                   cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCaa                                                          <  1:155338/97‑1 (MQ=255)
                                                                                                                                   cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCaa                                                          >  2:155338/1‑97 (MQ=255)
                                                                                                                                   cccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGAt     >  1:132134/1‑150 (MQ=255)
                                                                                                                                    ccGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATg    >  2:563065/1‑150 (MQ=255)
                                                                                                                                      ggTTCTCGCGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGgt  <  1:321531/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
ATCACCACCGAAGCCATGGTTGCCGACCTGCCGGAAGACAAGCCAGCTGCTGGCATGCCTGATATGGGTGGCATGGGTGGCATGGGCGGCATGATGTAAGCCAGCCTTACCCCCCTGCACCATCAAAAAGCCCCGGTTCGCGAGAACCGGGGCTTTTTTATTGCTGTCGGGGCCTCTTCGCAGCACAAGGCTGCGCTTGAGCCTTGTGCTGTGAAACGAGTGGATCAACCGACCAGCGGATCACCGACGTGCAGGATCTTCATGCCATTGGTGCCACCGATGGT  >  NC_002947/1551185‑1551468

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: