Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,718,092 C→A 9.2% L167L (CTG→CTT) ‡ paaG ← 2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,718,0920CA 9.2% 163.5 / 4.6 53L167L (CTG→CTT) ‡paaG2‑(1,2‑epoxy‑1,2‑dihydrophenyl)acetyl‑CoA isomerase
Reads supporting (aligned to +/- strand):  ref base C (19/29);  new base A (3/2);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

GCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGG  >  NC_002947/3717945‑3718238
                                                                                                                                                   |                                                                                                                                                  
gCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGcccg                                                                                                                                                   <  8:26461/149‑2 (MQ=255)
 cTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAg                                                                                                                                                  >  2:108210/1‑149 (MQ=255)
  tGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAg                                                                                                                                                  >  3:157282/1‑148 (MQ=255)
     gcTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATg                                                                                                                                              <  1:188721/149‑1 (MQ=255)
              gtgtCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGcc                                                                                                                                     <  3:226343/149‑1 (MQ=255)
                 tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCGAAGGCGCTCGCCCAGCATGGCCAGCGCCtt                                                                                                                                   <  5:348430/148‑1 (MQ=255)
                 tCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTg                                                                                                                                  <  4:157282/149‑1 (MQ=255)
                 tCAGGCCGTAGGTGGGAGGGGCGGCGAGGGGGCGGGCGAGGGTGATGGCTTCATCACGCAGCGCGGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGATCGCCCAGCATGGCCAGCGCCTTg                                                                                                                                  <  3:250598/149‑1 (MQ=255)
                     gCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTGGCCATGCTGGGCGAGCGCCTTggtgc                                                                                                                              <  1:166739/149‑4 (MQ=39)
                        gTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCggg                                                                                                                           <  1:210298/149‑1 (MQ=255)
                           ggTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGcc                                                                                                                         <  1:358348/148‑1 (MQ=255)
                           ggTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGcc                                                                                                                         <  5:378062/148‑1 (MQ=255)
                              gggcggggcgGCGAGGTGGCGGGCGAGGGCGCGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGc                                                                                                                     <  6:101730/149‑1 (MQ=255)
                                           aGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCgg                                                                                                        >  2:36817/1‑149 (MQ=255)
                                             gTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgca                                                                                                      <  6:6549/149‑1 (MQ=255)
                                                  gggcgAGGGTGCGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCagcagc                                                                                                 <  7:21008/149‑1 (MQ=255)
                                                  gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                                   <  4:148180/147‑1 (MQ=255)
                                                  gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGgcagca                                                                                                   >  3:148180/1‑147 (MQ=255)
                                                  gggcgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCagcagc                                                                                                 >  1:63170/1‑149 (MQ=255)
                                                     cgAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCaa                                                                                              <  5:222286/149‑1 (MQ=255)
                                                      gAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAg                                                                                             <  7:100675/149‑1 (MQ=255)
                                                          gTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAcc                                                                                         <  1:125405/149‑1 (MQ=255)
                                                            gAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgc                                                                                       >  1:160988/1‑149 (MQ=255)
                                                              gggCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgccg                                                                                     <  2:115426/149‑1 (MQ=255)
                                                              gggCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgccg                                                                                     <  2:160988/149‑1 (MQ=255)
                                                                gCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgcc                                                                                      <  8:88551/146‑1 (MQ=255)
                                                                gCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTAccgcc                                                                                      >  7:88551/1‑146 (MQ=255)
                                                                     atcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTcc                                                                               <  1:108210/148‑1 (MQ=255)
                                                                      tcaCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCGATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCgg                                                                             <  4:19164/149‑1 (MQ=255)
                                                                        aCGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCa                                                                           >  6:216030/1‑149 (MQ=255)
                                                                         cGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCAc                                                                          >  5:169346/1‑149 (MQ=255)
                                                                                cTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGcc                                                                    >  7:163138/1‑148 (MQ=255)
                                                                                          caccacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAg                                                          >  3:207702/1‑148 (MQ=255)
                                                                                           accacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAg                                                          >  6:106451/1‑147 (MQ=255)
                                                                                           accacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGACAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGTCCCGGGCCATGACGACCAAGCGCGGCAGCAGCCAAGTACCGCCTGAGTCCGGCACCAGGCCGATCATGCAg                                                          >  8:348289/1‑147 (MQ=255)
                                                                                            ccacTCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAg                                                       <  1:330261/149‑1 (MQ=255)
                                                                                            ccacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                             <  3:191136/143‑1 (MQ=255)
                                                                                            ccacGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTg                                                             >  4:191136/1‑143 (MQ=255)
                                                                                                            cccccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTgg                                       >  1:147118/1‑149 (MQ=255)
                                                                                                             ccccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGc                                      <  4:207702/149‑1 (MQ=255)
                                                                                                              cccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGcc                                     <  8:163138/149‑1 (MQ=255)
                                                                                                               ccATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGAAGAAGGCCTGGATGAAACTGGcc                                     >  8:277810/1‑148 (MQ=255)
                                                                                                                      tCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgggc                             >  5:47149/1‑149 (MQ=255)
                                                                                                                        ggCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAgcgggcgg                           >  1:262342/1‑149 (MQ=255)
                                                                                                                         gCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGggcggc                          <  2:63170/149‑1 (MQ=255)
                                                                                                                              ttCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAg                       <  5:106451/147‑1 (MQ=255)
                                                                                                                                ccgggcCAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGcc                         >  1:292114/5‑143 (MQ=38)
                                                                                                                                ccgggcCAAGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGcc                         <  2:292114/139‑1 (MQ=38)
                                                                                                                                 ggCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCAcc                   <  6:169346/148‑1 (MQ=255)
                                                                                                                                    gcCAATGCGCTCGCACAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAATCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGTCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGt               <  7:348289/149‑1 (MQ=255)
                                                                                                                                       aaGGCGCTGGCCATGCTGGGCGAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCg             >  8:136761/1‑148 (MQ=255)
                                                                                                                                         ggCGCTGGCCATGCTGGGCGAGCGCCTTGGCGCCGAACAGGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAg          >  6:359728/1‑149 (MQ=37)
                                                                                                                                         ggCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAg          <  2:262342/149‑1 (MQ=255)
                                                                                                                                            gcTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGc        >  8:14266/1‑148 (MQ=255)
                                                                                                                                                  ccAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAggg  >  3:281151/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTGCGCTTGATCAGTGTCAGGCCGTAGGTGGGCTGGGCGGCGAGGTGGCGGGCGAGGGTGAGGGCTTCATCACGCAGCGCTGCATCATCCACCACGCGGTAGATCAGCCCCCATTGCTCGGCCTGTTCGGCGCCAAGGCGCTCGCCCAGCATGGCCAGCGCCTTGGCCCGGGCCATGCCGACCAAGCGCGGCAGCAGCCAAGTACCGCCGGAGTCCGGCACCAGGCCGATCTTGCAGAAGGCCTGGATGAAACTGGCCGAGCGGGCGGCCAGCACCAGGTCGCAGGCCAAGGG  >  NC_002947/3717945‑3718238

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