Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB143 120,983 (T)9→8 intergenic (+280/+92) polA → / ← engB DNA polymerase I/putative GTP‑binding protein EngB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143120,9750T.94.9% 253.3 / 6.0 59intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base . (42/14);  total (43/16)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  ACB143/120854‑121105
                                                                                                                         |                                                                                                                                  
aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCgg                                                                                                                      >  3:347855/1‑135 (MQ=255)
   ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAtt                                                                                                                  >  2:291137/1‑136 (MQ=255)
           aaGCCCTTCGCTGCTCCTATGTTGGGGTAAGTGGTGGCAGATATCTGGCCCCCGCCCTCGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCggg                                                                                                           <  6:74269/135‑1 (MQ=255)
           aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGACCTAGCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCAGGAATTTGACgg                                                                                                            >  5:296432/1‑134 (MQ=255)
                           cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                          <  7:133641/136‑1 (MQ=255)
                             aTGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  1:104394/136‑1 (MQ=255)
                              tgttgtGTTAAGTGTTTGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  8:76049/135‑1 (MQ=255)
                              tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGCAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  8:161681/135‑1 (MQ=255)
                              tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  2:359596/135‑1 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  1:285598/1‑135 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  6:183032/1‑135 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  5:338713/1‑135 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  4:342024/1‑135 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                >  2:336961/1‑135 (MQ=255)
                                      tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGCGTTTTTTTTT‑ACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAAGTc                                                                                >  4:199061/1‑135 (MQ=255)
                                        aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                             >  6:166253/1‑136 (MQ=255)
                                        aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                             >  3:59695/1‑136 (MQ=255)
                                        aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                             >  2:318589/1‑136 (MQ=255)
                                           gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATa                                                                           >  3:191032/1‑135 (MQ=255)
                                           gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATa                                                                           >  2:133709/1‑135 (MQ=255)
                                           gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATa                                                                           >  3:318661/1‑135 (MQ=255)
                                             tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATacac                                                                        >  8:159075/1‑136 (MQ=255)
                                             tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATacac                                                                        >  5:183892/1‑136 (MQ=255)
                                                cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgc                                                                     >  7:205631/1‑136 (MQ=255)
                                                cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgc                                                                     >  7:317759/1‑136 (MQ=255)
                                                cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgc                                                                     >  2:318426/1‑136 (MQ=255)
                                                 aGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                    >  6:28716/1‑136 (MQ=255)
                                                      tCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAtt                                                               >  4:268159/1‑136 (MQ=255)
                                                       cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTc                                                              >  5:295947/1‑136 (MQ=255)
                                                            ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCggg                                                         >  4:242630/1‑136 (MQ=255)
                                                                  ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                   >  5:320682/1‑136 (MQ=255)
                                                                  ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                   >  5:300768/1‑136 (MQ=255)
                                                                  ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                   >  5:152754/1‑136 (MQ=255)
                                                                  ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                   >  7:114015/1‑136 (MQ=255)
                                                                    tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                 >  8:108889/1‑136 (MQ=255)
                                                                      gCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAACTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                <  4:294250/135‑1 (MQ=255)
                                                                            cGGACTTGGACCCCGAACTTCCCCCTCCCCCTACCAACTCCGGGG‑TTTTTTTTGCCCGGAAATTGCCGGGAGCGGGTTTTCTCGCTAATTCTATGTCGAATAAACCTCAGAATTCGCGGGTAAACCCGGGTCAGgc                                         <  6:151150/136‑1 (MQ=255)
                                                                            cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         <  8:181727/136‑1 (MQ=255)
                                                                            cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         <  3:309169/136‑1 (MQ=255)
                                                                             ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                        <  5:187755/136‑1 (MQ=255)
                                                                               aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcgc                                       >  7:271220/1‑135 (MQ=255)
                                                                               aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcgc                                       >  3:296258/1‑135 (MQ=255)
                                                                                       cccGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACt                               >  6:147883/1‑135 (MQ=255)
                                                                                       cccGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACt                               >  1:49157/1‑135 (MQ=255)
                                                                                              ttCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCCGGCGCAAACCACTAAACGAcc                       >  5:348002/1‑136 (MQ=255)
                                                                                                cccccTCCCCATACGATGTCC‑GGGTTTTTTATGGCCCGTAATTGGCCGGGAGGGGGTTTTCTCGCTAATGCTATGTCGAATACCCCGCGGCATTGGCGGGAAAACCCGGGTCAGGCGCAAACCACAAAACGACCgg                     <  4:245918/136‑1 (MQ=255)
                                                                                                    cTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTg                 >  7:315337/1‑136 (MQ=255)
                                                                                                    cTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTg                 >  3:160945/1‑136 (MQ=255)
                                                                                                      ccccATCCGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCTAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCGGGGTCAGGCGAAAACCACTAAACGACCGGCTTGTc               <  6:341030/136‑1 (MQ=255)
                                                                                                      ccccATACTAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  2:118786/136‑1 (MQ=255)
                                                                                                      ccccATACGACGTCC‑GGGTTTTTTTTTGCCCGTAATTTGCCGGTAGCGGGTTTTCTCGCGAATGCTATGGCGAATACACCCCCTCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  4:267193/136‑1 (MQ=255)
                                                                                                      ccccATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGTGCGGGTTTTCTCGCTAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  4:16101/136‑1 (MQ=255)
                                                                                                          aTACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTc            >  1:56368/1‑135 (MQ=255)
                                                                                                          aTACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTc            >  4:176255/1‑135 (MQ=255)
                                                                                                          aTACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTc            >  3:139689/1‑135 (MQ=255)
                                                                                                           tACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTACACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTc            >  8:102860/1‑134 (MQ=255)
                                                                                                              gAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGAGGGTAAACCCGGGTCAGGCGCAAACCACAAAACGACAGGCTGCTCTTCCAGtt       >  6:128869/1‑136 (MQ=255)
                                                                                                                 gTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtcca    >  7:80473/1‑136 (MQ=255)
                                                                                                                  tCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccat   >  6:178525/1‑136 (MQ=255)
                                                                                                                     cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        >  2:323572/2‑49 (MQ=255)
                                                                                                                     cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  1:323572/48‑1 (MQ=255)
                                                                                                                     cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc  >  3:96953/2‑135 (MQ=255)
                                                                                                                         |                                                                                                                                  
AATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  ACB143/120854‑121105

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