Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,868,469 C→T 6.9% G32D (GGT→GAT)  lgt ← prolipoprotein diacylglyceryl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,868,4690CT 6.9% 171.0 / 2.0 58G32D (GGT→GAT) lgtprolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base T (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

GTAGGCGTGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGAACAA  >  NC_002947/5868321‑5868612
                                                                                                                                                    |                                                                                                                                               
gTAGGCGTGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCAc                                                                                                                                                 >  6:153428/1‑149 (MQ=255)
       tGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCgccgat                                                                                                                                          <  5:89393/149‑1 (MQ=255)
         cAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCgccgatgc                                                                                                                                        <  5:255456/149‑1 (MQ=255)
              cGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATc                                                                                                                                   <  7:185932/149‑1 (MQ=255)
               gTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTCCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCgccgatgccgat                                                                                                                                    >  3:100153/1‑147 (MQ=255)
                tAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAg                                                                                                                                 >  5:240001/1‑149 (MQ=255)
                  gAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACCCCCGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGt                                                                                                                               <  5:126613/149‑1 (MQ=255)
                  gAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAgg                                                                                                                                >  2:240861/1‑148 (MQ=255)
                        cacaTAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACAt                                                                                                                           >  1:163800/1‑147 (MQ=255)
                        cacaTAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACAt                                                                                                                           >  1:169488/1‑147 (MQ=255)
                           aTAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGc                                                                                                                      <  2:13816/149‑1 (MQ=255)
                                        ccgccgACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGAt                                                                                                         >  7:97537/1‑149 (MQ=255)
                                           ccgACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGAtttt                                                                                                      <  5:68027/149‑1 (MQ=255)
                                                attCACCTCGATCGCCAGCCAGACCCGGAGGTCGGACAGCTTCGCTCGGTTGCAGGTGGGGTCGAAGGGTTTGAGCCGGCGTGAGGCCAGCTGCCAGGGACCGCCGATGCGGATCAGGTACATCATGCCGTACCAATGGATTTTCAggg                                                                                                 <  4:219668/147‑3 (MQ=255)
                                                  tCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGc                                                                                               >  7:47023/1‑149 (MQ=255)
                                                   cACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTCCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGcc                                                                                              <  4:100153/149‑1 (MQ=255)
                                                     cccccATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCg                                                                                            <  1:122351/149‑1 (MQ=255)
                                                        ccATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATg                                                                                         <  1:240861/149‑1 (MQ=255)
                                                              aCAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   >  4:15608/1‑149 (MQ=255)
                                                                aGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   >  5:251103/1‑147 (MQ=255)
                                                                aGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   <  6:251103/147‑1 (MQ=255)
                                                                        cacCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGt                                                                         <  8:97537/149‑1 (MQ=255)
                                                                        cacCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGt                                                                         <  3:15608/149‑1 (MQ=255)
                                                                         acCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTc                                                                        <  2:36151/149‑1 (MQ=255)
                                                                          ccAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCTTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCt                                                                       <  5:153428/149‑1 (MQ=255)
                                                                           cAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACg                                                                             >  1:17296/1‑142 (MQ=255)
                                                                           cAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACg                                                                             <  2:17296/142‑1 (MQ=255)
                                                                           cAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTa                                                                      >  5:332297/1‑149 (MQ=255)
                                                                            aGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTAt                                                                     <  6:240001/149‑1 (MQ=255)
                                                                                 ggACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGc                                                                 <  8:47023/148‑1 (MQ=255)
                                                                                 ggACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCg                                                                >  4:293140/1‑149 (MQ=255)
                                                                                    cAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCggg                                                              >  1:114178/1‑148 (MQ=255)
                                                                                        ttCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAgg                                                         >  2:262505/1‑149 (MQ=255)
                                                                                          ctcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAgg                                                         <  3:293140/147‑1 (MQ=255)
                                                                                          ctcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGgc                                                        <  3:339466/148‑1 (MQ=255)
                                                                                           tcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATg                                                                                                            >  8:320872/1‑95 (MQ=255)
                                                                                           tcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATg                                                                                                            <  7:320872/95‑1 (MQ=255)
                                                                                                   ccAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAAc                                              >  3:258044/1‑149 (MQ=255)
                                                                                                               gAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTc                                           >  2:250087/1‑140 (MQ=255)
                                                                                                               gAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTc                                           <  1:250087/140‑1 (MQ=255)
                                                                                                                aaGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   <  3:151886/99‑1 (MQ=255)
                                                                                                                aaGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   >  4:151886/1‑99 (MQ=255)
                                                                                                                         aGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATtcaaatcaaa                         >  5:127425/1‑148 (MQ=255)
                                                                                                                          gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATg                                                                                         <  5:110069/83‑1 (MQ=255)
                                                                                                                          gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATg                                                                                         >  6:110069/1‑83 (MQ=255)
                                                                                                                          gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCaaaaa                       <  2:163800/149‑1 (MQ=255)
                                                                                                                          gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCaaaaa                       <  2:169488/149‑1 (MQ=255)
                                                                                                                           ccGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAg                      >  4:282901/1‑149 (MQ=255)
                                                                                                                                  gAGGCCAGCAGCCAGGCACCGCCGATGCCGTTCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGc               >  1:33265/1‑149 (MQ=255)
                                                                                                                                    ggCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCa             <  4:258044/149‑1 (MQ=255)
                                                                                                                                     gCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGa                  <  3:295098/143‑1 (MQ=255)
                                                                                                                                     gCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGa                  >  4:295098/1‑143 (MQ=255)
                                                                                                                                          cagcCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAg       <  2:114178/149‑1 (MQ=255)
                                                                                                                                            gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                >  6:67610/5‑106 (MQ=255)
                                                                                                                                            gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                <  7:116694/102‑1 (MQ=255)
                                                                                                                                            gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                <  5:67610/102‑1 (MQ=255)
                                                                                                                                            gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                >  8:116694/5‑106 (MQ=255)
                                                                                                                                             ccAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGaa     >  3:220025/1‑148 (MQ=255)
                                                                                                                                               aGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGACCTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGaacaa  >  8:328047/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTAGGCGTGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGAACAA  >  NC_002947/5868321‑5868612

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