Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,723,150 G→A 7.2% R136H (CGC→CAC)  PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,723,1500GA 7.2% 157.4 / 2.2 55R136H (CGC→CAC) PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (2/2);  total (27/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.60e-01

AACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGACATTG  >  NC_002947/3723004‑3723298
                                                                                                                                                  |                                                                                                                                                    
aaCCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAACCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTgcgc                                                                                                                                                     >  2:212613/1‑148 (MQ=255)
         ggtgTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCACGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGaa                                                                                                                                           <  4:6906/149‑1 (MQ=255)
          gtgtGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGGTCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGaa                                                                                                                                           <  2:30672/148‑1 (MQ=255)
               ggTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCt                                                                                                                                     >  3:320555/1‑149 (MQ=255)
                     gcagcCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGATGCGCGAAACGCGAATACCCTTGCAcc                                                                                                                               >  4:249977/1‑149 (MQ=255)
                      cagcCCTTGAATGAGATCCTGATACACTTGCGTGAGAAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCt                                                                                                                              >  1:6001/1‑149 (MQ=255)
                       agcCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTg                                                                                                                             <  7:278488/149‑1 (MQ=255)
                        gcCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTgg                                                                                                                            <  6:254200/149‑1 (MQ=255)
                        gcCCTTGAATGAGATCCTGATACACTTGCGTGAACAGCCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTgg                                                                                                                            <  6:317409/149‑1 (MQ=255)
                         cccTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATa                                                                                                                                         <  5:333496/135‑1 (MQ=255)
                         cccTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATa                                                                                                                                         >  6:333496/1‑135 (MQ=255)
                           cTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCACGAAACGCGAATACCCTTGCACCTGGTa                                                                                                                          >  8:193145/1‑148 (MQ=255)
                           cTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCACGAAACGCGAATACCCTTGCACCTGGTa                                                                                                                          >  6:150046/1‑148 (MQ=255)
                              gAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACATGGTacac                                                                                                                       >  3:302951/1‑148 (MQ=255)
                                aTGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACaa                                                                                                                     <  3:164136/148‑1 (MQ=255)
                                aTGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAg                                                                                                                    <  6:122332/149‑1 (MQ=255)
                                 tGAGATCCTGATACACTTGCGTGAGGAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAgg                                                                                                                   <  2:6001/149‑1 (MQ=255)
                                   agaTCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGgc                                                                                                                  <  6:31391/148‑1 (MQ=255)
                                      tCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGgccagc                                                                                                              >  1:233689/1‑149 (MQ=255)
                                       ccTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGgccagc                                                                                                              >  5:44405/1‑148 (MQ=255)
                                          gATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATg                                                                                                          >  5:166146/1‑149 (MQ=255)
                                           aTACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGc                                                                                                         >  5:295440/1‑149 (MQ=255)
                                            tACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCg                                                                                                        <  5:8269/149‑1 (MQ=255)
                                             acacTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTac                                                                                                                         >  2:75339/1‑131 (MQ=255)
                                             acacTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTac                                                                                                                         <  1:75339/131‑1 (MQ=255)
                                                 ttGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGc                                                                                                    >  7:14723/1‑148 (MQ=255)
                                                 ttGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGc                                                                                                    >  1:295021/1‑148 (MQ=255)
                                                            aCCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCATTGTTCACTTCGCTCGATATCCAGTTGCTGCGCGAACCGCGAATACCATAGCACCTGCTACACAAGGCCAGCCATCCGTTGCCACGCAGGGTTc                                                                                        >  5:89449/1‑149 (MQ=255)
                                                                   gCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTggcggc                                                                                 >  2:333450/1‑149 (MQ=255)
                                                                    ccATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGgcggcg                                                                                >  8:15603/1‑149 (MQ=255)
                                                                      aTGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGt                                                                              <  6:44405/149‑1 (MQ=255)
                                                                       tGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTg                                                                             <  2:233689/149‑1 (MQ=255)
                                                                         cTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGa                                                                           <  2:295021/149‑1 (MQ=255)
                                                                         cTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGa                                                                           <  1:42260/149‑1 (MQ=255)
                                                                         cTGATCAAAGAGCTCGAGCATGAACCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGa                                                                           <  1:212613/149‑1 (MQ=255)
                                                                          tGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGAc                                                                          >  5:317957/1‑149 (MQ=255)
                                                                                 agagCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCa                                                                   <  4:302951/149‑1 (MQ=255)
                                                                                      tCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGAACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCg                                                              >  1:222284/1‑149 (MQ=255)
                                                                                         aGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCggag                                                           <  8:14723/149‑1 (MQ=255)
                                                                                          gCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCggag                                                           >  5:188437/1‑148 (MQ=255)
                                                                                              gAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGAt                                                        >  3:148527/1‑147 (MQ=255)
                                                                                                      gtggATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGaa                                               <  7:15603/148‑1 (MQ=255)
                                                                                                      gtggATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAg                                              >  8:192835/1‑149 (MQ=255)
                                                                                                      gtggATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCACGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGaa                                               <  5:150046/148‑1 (MQ=255)
                                                                                                      gtggATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCACGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGaa                                               <  7:193145/148‑1 (MQ=255)
                                                                                                      gtggATGCGCGAATTTTTCACTTCGGTGGATATCCAGTTGCTGCGCTAACCGCGAATACCCTTGCACGTGGTACACATGGCCAGCCATGCGTTGCCAGGCAGTTTTCTGGCGGCGGTGTACCTGTCCAGGCCGGAGGATCAGTTTGAAg                                              <  6:295440/149‑1 (MQ=255)
                                                                                                             cgcgcAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCaa                                       <  7:192835/149‑1 (MQ=255)
                                                                                                              gcgcAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCggag                                                           >  4:308860/1‑128 (MQ=255)
                                                                                                              gcgcAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCggag                                                           <  3:308860/128‑1 (MQ=255)
                                                                                                                         aCTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGatca                            <  6:166146/148‑1 (MQ=255)
                                                                                                                                  gATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAgc                  >  8:31339/1‑149 (MQ=255)
                                                                                                                                    tatCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAgcgc                <  1:333450/149‑1 (MQ=255)
                                                                                                                                       ccAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTg              >  5:306110/1‑148 (MQ=255)
                                                                                                                                              ctgcGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGa       <  6:188437/148‑1 (MQ=255)
                                                                                                                                               tgcGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGACa     >  1:204773/1‑149 (MQ=255)
                                                                                                                                               tgcGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTCCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGACa     >  8:104328/1‑149 (MQ=255)
                                                                                                                                                 cgcgAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGACAtt   <  7:31339/149‑1 (MQ=255)
                                                                                                                                                  gcgAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGACATTg  >  8:118536/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGATACACTTGCGTGAGCAACCGTTTGCCATGCTGATCAAAGAGCTCGAGCATGAGCCTTGGTGGATGCGCGCAATGTTCACTTCGCTGGATATCCAGTTGCTGCGCGAAACGCGAATACCCTTGCACCTGGTACACAAGGCCAGCCATGCGTTGCCACGCAGGGTTCTGGCGGCGGTGGACCTGTCCAGGCCGGAGGATCAGTTTGAAGGCTTCAACGACCAGATCATCAGCGAAGCGCTGAAGCTGACATTG  >  NC_002947/3723004‑3723298

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