Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 924,303 T→G 27.0% D4790A (GAC→GCC)  CJ019_00881 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122924,3030TG27.0% 75.8 / 5.9 63D4790A (GAC→GCC) CJ019_00881hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base G (6/13);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTGCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGT  >  ACB122/924173‑924431
                                                                                                                                  |                                                                                                                                
ttGTGCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCCGGCCGCTGACCACGCTGCCGTCGGc                                                                                                                              >  4:188457/1‑135 (MQ=255)
ttGTGCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGc                                                                                                                              >  5:428749/1‑135 (MQ=255)
 tgtgCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGc                                                                                                                              <  1:220107/134‑1 (MQ=255)
    gCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGtt                                                                                                                           <  4:44415/134‑1 (MQ=255)
    gCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGtt                                                                                                                           <  7:133757/134‑1 (MQ=255)
    gCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGtt                                                                                                                           <  2:65795/134‑1 (MQ=255)
     cGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGa                                                                                                                         <  5:221564/135‑1 (MQ=255)
         tGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATcac                                                                                                                     <  1:276768/135‑1 (MQ=255)
                           gCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTc                                                                                                   <  1:116969/135‑1 (MQ=255)
                                cATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGGCGGCGTTGATCACCACTTGGGTGAGTGGGTcgggcg                                                                                              <  8:13468/135‑1 (MQ=255)
                                 aTCTGCAAGTACGAAACCCGTGCCCCGTACGGTACCGGCCGCGTTACTGACGATTCCCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCGGTCGGCGTTGATCACCACTTGGGTGAGTGGGTcgggcgg                                                                                             <  7:190986/135‑1 (MQ=255)
                                 aTCAGCCAGTACGATACCCGGGCCCAGTACGGTACCGGCCGCGTTACTGCCGATTACCGTGGCGCCGGCCGCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTcgggcgg                                                                                             <  3:319482/135‑1 (MQ=255)
                                           aCGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGAGGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTcc                                                                                   <  6:397391/135‑1 (MQ=255)
                                           aCGATACCCGGGCCCCGGACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCGGCGCCCAGGCCGCTGACCAGGCTGCCGGCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTcc                                                                                   <  8:175113/135‑1 (MQ=255)
                                            cGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCCTGGCGCCGGCCTCGCCCCCGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCg                                                                                  >  2:90756/1‑135 (MQ=255)
                                                  ccGTGCCCAGTACGGTACCGGCCGCGTTTATGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACCCTGCCGTCGGCGTTGATGACCACGTGGGTGAGTGGGTGGGGGGGGGTGAAGGACGGGGttt                                                                            >  8:436923/1‑135 (MQ=255)
                                                   cGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGCCCACGCTGCCGGCGGCGTTGATAACCACGTGGGTGTGTGGGTCGGGCGGGGTGAAGGCCGGGGGTTg                                                                           >  8:197986/1‑135 (MQ=255)
                                                       cccAGTACGGTACCGGCCGCGTTACTGACGATTTCCCTGGCGCCGGCCACGCCCAGGCCCCCGCCCACCCTGGCGGCGGGGTGTGTGCCCCCATGGGGGTGTTGGTGGGGCGGGGGTGATTCCGGGGGGTGCAgg                                                                       >  6:53098/1‑135 (MQ=255)
                                                       cccAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCCGGCCGCGGCCCACGCTGCCGGCGGCGGTGATCCCCACATGGGTGAGTGGGTCGGGGGGGGTGGAGTCCGGGGGTTGGAgg                                                                       >  7:390949/1‑135 (MQ=255)
                                                       cccAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTcgggcggg                                                                                            <  3:220134/114‑1 (MQ=255)
                                                       cccAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTcgggcggg                                                                                            >  4:220134/1‑114 (MQ=255)
                                                       cccAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGGTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAgg                                                                       >  3:103758/1‑135 (MQ=255)
                                                           gTACGGTACCGGCGGCGTTACTGACGATTACCGTGGCGCCGGCCGCGCGCTGGCCGCGGACCAGGGTGCCGGCGGCGTTTATCACCACTTGGGTGAGTGGGGCGGGCGGGGGGAAGTCCGGGGTTTGCAGGGTAg                                                                   <  8:82696/135‑1 (MQ=255)
                                                               ggTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCCCGCTGCCGTCGGCGTTGATCACCACTTGGGGGAGTGGGGCGGGCGGGGGGAAGTGCGGGGTTTGCAGGGTAGCCAc                                                               >  8:246766/1‑135 (MQ=255)
                                                                  aCCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCGCGCCCAGGCCGCTGACCCCCCTGCCGGCCGCGTTGATCACCCCGTGGGTGAGTGGGTGGGGCGGGGGGAAGTCCGGGGGTTGGAGGGTAGCCACAGc                                                            >  4:14260/1‑135 (MQ=255)
                                                                    cGGCCGGGTTACTCACGATAACCGTGGCGCCGGCCTCGACCAGGCCGGTGACCACGCTGCCGGCGGCGTTAATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCg                                                          <  4:389745/135‑1 (MQ=255)
                                                                    cGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCCCCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGGTTGCAGGGTAGCCACTGCCg                                                          >  8:61334/1‑135 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCCGGCCGCTGACCCCGCTGCCGTCGGCGTTGATAACCACTTTGGTGAGTGGGTCGGGCGGGGTGAATTCCGGGGTTTGCAGGGTAGCCCCAGCCGACACg                                                     >  4:51059/1‑135 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGCCCACCCTGCCGGCGGCGTTGATCACCCCGTGGGGGAGTGGGGGGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACg                                                     >  7:221653/1‑135 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCCCGCTGCCGTCGGCGTTGATAACCACGTGGGGGAGGGGGGCGGGCGGGGGGGAGTCCGGGGGTTGGAGGGTTGCCCCGGCCgccccg                                                     >  1:141019/1‑130 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCCCGCTGCCGGCGGCGTTGATGCCCACTTGGGTGATTGGGTCGGGCGGGGTGTAGTCCGGGGGTTTCAGGGTAGCCCCAGCCGccacc                                                     >  4:289796/1‑134 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGtt                                                                                                                           <  5:216314/65‑1 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGtt                                                                                                                           >  6:216314/1‑65 (MQ=255)
                                                                         gcgTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGGTGATCACCACGTGGGTGAGTGGGTCGGGGGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCAATGCCGCCACg                                                     >  8:53039/1‑135 (MQ=255)
                                                                             tACTGACGATTCCCGTGGCGCCGGCCTCGCCCAGGCCGCGGACCACGCTGCCGGCGGCGTTGAGCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTgc                                                 <  8:390949/135‑1 (MQ=255)
                                                                                  aCGATTAGCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCAGGCTGCCGTCGGGGTCGATCACCACTTGGGTGAGTGGGTCGGGCGTGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTgccggcc                                            <  7:322975/135‑1 (MQ=255)
                                                                                   cGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCGGACCACGCTGCCGGCGGCGTTGATCACCACTTGGGGGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTgccggccg                                           <  4:118280/135‑1 (MQ=255)
                                                                                       tACCGGGGCGCCGGCCGCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGGGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATc                                       <  3:51059/135‑1 (MQ=255)
                                                                                        aCCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCg                                      >  1:91593/1‑135 (MQ=255)
                                                                                               cgcCGGCCTCGCCCAGGCCGCTGCCCACGCTGCCGTCGGCGGTGATCACCCCTTGGGTGAGTGGGGCGGGCGGGGTGCAGTCCGGGGGTTGCAGGGTAGCCCCTGCCGACACGTTGCCGGCCGCATCGGCCgggg                               >  5:351982/1‑135 (MQ=255)
                                                                                               cgcCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGa                                                                                                                         <  8:285784/45‑1 (MQ=255)
                                                                                               cgcCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGa                                                                                                                         >  7:285784/1‑45 (MQ=255)
                                                                                               cgcCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCCTCGGCCTggg                               >  8:114274/1‑135 (MQ=255)
                                                                                               cgcCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGGTTGCAGGGTAGCCACTGCCGACCCGTTGCCGGCCGCCTCGGCCgggg                               >  3:36043/1‑135 (MQ=255)
                                                                                                 ccGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGt                              >  2:330743/1‑134 (MQ=255)
                                                                                                           cccggCCGCTGCCCACGGTGCGGTCGGGGTTGATCACCACTTGGGTGTGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGagca                   <  3:449188/132‑1 (MQ=255)
                                                                                                           cccggCCGCTCCCCACGGGGCGGGCGGCGTGTATCACCCCTTGGGTGTGGGGGGCGGGCGGGGGGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGagca                   <  5:53098/132‑1 (MQ=255)
                                                                                                           cccggCCGCGGCCCACGCTGCCGGCGGCGTTGAGCACCCCTTGGGTGCGTGGGGCGGGCGGGGGGAAGGCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGagca                   <  3:291849/132‑1 (MQ=255)
                                                                                                           ccAGGCCGCTGACCACGCTGCCGTCGGCTTTGAGCACCACTGGGGGGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTGGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGagca                   <  1:354805/135‑1 (MQ=255)
                                                                                                                ccGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGGCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCg              <  1:90756/135‑1 (MQ=255)
                                                                                                                 cGCTGACCAGGCGGCCGGCGGCGTTGATCACCCCTTGGGTGAGGGGGGCGGGCGGGGGGAAGTCCGGGGGTTGCAGGGGAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCgg             <  5:82050/135‑1 (MQ=255)
                                                                                                                 cGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGCCACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCgg             >  4:104944/1‑135 (MQ=255)
                                                                                                                  gCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGCGCAGCTCGGc            >  6:6243/1‑135 (MQ=255)
                                                                                                                  gCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGc            >  5:264411/1‑135 (MQ=255)
                                                                                                                     gcccccgcgGCCGTCGGCGTTGCTCCCCACTTGGGTGAGTGGTGCGGGCGGGGTGAAGGCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGGGAAGCTGAGCAGCTCGGCGtt         <  5:153432/130‑1 (MQ=255)
                                                                                                                     gACCAGGGTGCCGGCGGCTTGGATCACCGCTTGGGGGAGGGGGGCGGGCGGGGTGAAGGCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGtt         <  5:381332/135‑1 (MQ=255)
                                                                                                                     gACCACGCGGCCGTGGGCGTTGATCACCTCTTGGGGGAGTGGGGCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGtt         <  3:188457/135‑1 (MQ=255)
                                                                                                                      aCCACGCTGCCGTCGGCGTTGATCACCCCTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGGCTGGGTGAAGCTGCGCAGCTCGGCGTTg        >  6:350272/1‑135 (MQ=255)
                                                                                                                      aCCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGCAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGGGAAGCTGCGCAGCTCGGCGTTg        >  8:335020/1‑135 (MQ=255)
                                                                                                                       ccACGCTGCCGGCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGc       <  8:384626/135‑1 (MQ=255)
                                                                                                                           gcggCCGTCGGCGTTCAGAACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAg   <  7:424680/132‑1 (MQ=255)
                                                                                                                           gCTGCCGGGGGCGTTCATCCCCACTGGGGTGAGTGGGGCGGGCGGGGGGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAg   <  1:438277/135‑1 (MQ=255)
                                                                                                                            cTGCCGGCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGGGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGt  <  3:296597/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                
TTGTGCGGGTGACAGTTCGACGCTGAAGCTGCCATCAGCCAGTACGATACCCGTGCCCAGTACGGTACCGGCCGCGTTACTGACGATTACCGTGGCGCCGGCCTCGCCCAGGCCGCTGACCACGCTGCCGTCGGCGTTGATCACCACTTGGGTGAGTGGGTCGGGCGGGGTGAAGTCCGGGGTTTGCAGGGTAGCCACTGCCGACACGTTGCCGGCCGCATCGGCCTGGGTGAAGCTGAGCAGCTCGGCGTTGCCCAGT  >  ACB122/924173‑924431

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