Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,277,898 A→G T144A (ACC→GCC)  soxR → DNA‑binding transcriptional dual regulator SoxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,277,8980AG97.1% 98.9 / ‑4.2 34coding (430/465 nt)soxRDNA‑binding transcriptional dual regulator SoxR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/17);  minor base . (0/1);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTCTTCGCCAGTGTAAAAGTATACTTTTTAACCGCAATATTTTT  >  NC_000913/4277656‑4278134
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            
cTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGgcgc                                                                                                                                                                                                                                       >  1:145189/1‑250 (MQ=255)
  cAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATTCCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTATCCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCAc                                                                                                                                                                                                                                     <  1:17997/250‑1 (MQ=255)
     cGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCCTAACCCGGGCGACCGCTAAGGAGAAGAAGGGGCCGGCGCCCGCtt                                                                                                                                                                                                                                 >  2:180295/1‑251 (MQ=255)
        aTTGGCATTCCGCTGGCGACCATTGGTGAAGCGTGTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGc                                                                                                                                                                                                                               <  2:62928/250‑1 (MQ=255)
           ggCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGa                                                                                                                                                                                                                           <  2:216144/251‑1 (MQ=255)
            gCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGa                                                                                                                                                                                                                           >  1:59426/1‑250 (MQ=255)
             cATTCCGCTGACGACCATTGGTGAAGCGGTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGATGGATCGGCGCATTCATACCTTAGCGGCTCTGCGTGACGAACTGGCCGGATGAATTCGTTGTGGCTGCCTTTCGCGCCAAGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGCGAAGAAGGAGCCGGCGCACGCTTGCTGGa                                                                                                                                                                                                                           >  2:276351/1‑249 (MQ=255)
                                          ttGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCTCGCGCTTTCTGGAAGATGAACAAAACTAAAGCGCCACAAGGgc                                                                                                                                                                                             <  1:14737/250‑1 (MQ=255)
                                                    gCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGGCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCGTAGTGGCGCGGCGTTACGAACGGGACGGATGTAGGGGTTGGGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGGAACCCGGGCGACCGCTTAGAAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAATCGCCACAAGGGCGCTTTAggttg                                                                                                                                                                                  <  2:30592/251‑1 (MQ=255)
                                                     cccGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAgtttgt                                                                                                                                                                                 <  1:242438/251‑1 (MQ=255)
                                                         aaGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTc                                                                                                                                                                             <  2:244876/251‑1 (MQ=255)
                                                                  cGTTAAGTGCGAAATAGGGGAAACAGCTTTCGTCTCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGGGGCGCTGCGTGACGAACTGGACGGATGTATTGGTGGTGGCTGCCTTTCGCGCAGTGATTGTCCGTTGCGTAATCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTTCTGGAAGATGAACAAAACTAAAGCGCCACTAGGTCGCTTTAGTTTGTTTTCCGgtctttg                                                                                                                                                                    <  2:81413/251‑1 (MQ=255)
                                                                   ggtAAGTGCGCAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCTGCGCACGCTTGCTGGAATATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgt                                                                                                                                                                   <  2:237411/249‑1 (MQ=255)
                                                                    ttAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGgtctttgtc                                                                                                                                                                  >  1:241223/1‑251 (MQ=255)
                                                                          gCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTtct                                                                                                                                                             <  1:50658/250‑1 (MQ=255)
                                                                             aaaGAATGGAAACCGCTTTCGTCACTATGGAGAGAAGAGTTGGATACGCGAATTAATACCTTACTAGCGCTCCGTGAAGAAATGGACGGATGTACTGGATGTGGCTGCCTTTCGCGCAGTGACTGCCCGAAGCGTAAACCGGGCGACAGCTCTGGAGAAGAAGGTGACGGAGCCCGCTTGCTGGAAGATGAACCAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTACGGTCTTTTTCTTTcactct                                                                                                                                                         >  2:283629/1‑251 (MQ=255)
                                                                                      aaaCAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTg                                                                                                                                                <  2:72527/251‑1 (MQ=255)
                                                                                                 gTCCCAATGGCGAGCAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACACGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACagg                                                                                                                                       >  2:287775/1‑249 (MQ=255)
                                                                                                        tGGCGAGTAAAGATGGATCGCCGCATTCATACCTTAGTGGCGCTGCGTGACGACCTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTCGGAGGAGACGGTGCCGGCGCCCGCTTGCTGGAAGATGAACAAAACTAAAGCGCCCCAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTcc                                                                                                                              >  2:255249/1‑251 (MQ=255)
                                                                                                                      ggATCGGCGCATTCATTCCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacac                                                                                                                <  2:59426/251‑1 (MQ=255)
                                                                                                                          cGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCt                                                                                                            <  1:130956/251‑1 (MQ=255)
                                                                                                                                     tACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCAc                                                                                                 >  2:162609/1‑251 (MQ=255)
                                                                                                                                             tGGCGCTGCGTGCCGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCCTTCGAGCCCgtgggggg                                                                                         >  2:35813/1‑247 (MQ=255)
                                                                                                                                              ggCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCCTAACCCGGGCGACCGCTAAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGa                                                                                         <  1:180295/250‑1 (MQ=255)
                                                                                                                                                             acggtaCGGATGTATTGGGTGTGGCTGCCTTTCGCTCAGTGATTGCCCGTTGCTTATCCCTGGCGATCGCTTAGGAGAAGAAGTGTCCGGCGCTCGCTTGCTTGAAGTTGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGGCTTTGTCTTTCTCGCTATCCCGCTGTTGCACAGTAGGTTTTCCCCCGTCGTCATCACTCCTCATTCTAGCACgtggtgtattgtccggttggtttg                                                                         <  2:273946/246‑19 (MQ=255)
                                                                                                                                                              cTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTg                                                                        >  1:85288/1‑251 (MQ=255)
                                                                                                                                                              cTGGACGGAAGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTCACCCGGGCGACCGCTTAGGAGCAGAAGGTGCCGGCGCACGCTTGCTGGAAGCTGAACAAAACTAAAGCGCCACAAGGGAGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTAGAGCACGTGGAGGAGGTTCCGGTTGGTGTTg                                                                        >  2:254418/1‑251 (MQ=255)
                                                                                                                                                                  aCGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGc                                                                    <  1:243441/251‑1 (MQ=255)
                                                                                                                                                                                 gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacaca                                                                                                               <  2:71549/193‑1 (MQ=255)
                                                                                                                                                                                 gtgGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAacaca                                                                                                               >  1:71549/1‑193 (MQ=255)
                                                                                                                                                                                            tctcgcAGTGATTGACAGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGCAGGTGCCGGCGCACGCTTGCTGGAAGATGCACAAAACTAACGCGCTACAAGGGCGCTTCAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATTCCGCTGGTACACAGGCGGGTTTCCCCCGACGTCAACACAACTCATTCGAGCACGTGGTGGAGGTTCCTGTTGGTGTTGAGGCTTTAATTGTATGTCACCGACGCttct                                          >  2:2657/4‑251 (MQ=255)
                                                                                                                                                                                                gcAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCAGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTCTTCGc                                      >  1:210941/1‑251 (MQ=255)
                                                                                                                                                                                                    tGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACACGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTCTTCGCCAgt                                  >  1:175514/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aGGAGAAGAAGGTGCCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTCTTCGCCAGTGTAAAAGTATACTTTTTAACCGCAATAttttt  <  2:63988/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                            
CTCAGCGTATTGGCATTCCGCTGGCGACCATTGGTGAAGCGTTTGGCGTGTTGCCCGAAGGGCATACGTTAAGTGCGAAAGAGTGGAAACAGCTTTCGTCCCAATGGCGAGAAGAGTTGGATCGGCGCATTCATACCTTAGTGGCGCTGCGTGACGAACTGGACGGATGTATTGGTTGTGGCTGCCTTTCGCGCAGTGATTGCCCGTTGCGTAACCCGGGCGACCGCTTAGGAGAAGAAGGTACCGGCGCACGCTTGCTGGAAGATGAACAAAACTAAAGCGCCACAAGGGCGCTTTAGTTTGTTTTCCGGTCTTTGTCTTTCTCTCTATCCCGCTGGTACACAGGAGGGTTTCCCCCGACGTCAACACACCTCATTCGAGCACGTGGTGGAGGTTCCGGTTGGTGTTGATGCTTTAATTGTATGTCACCGACGTTTCTTCGCCAGTGTAAAAGTATACTTTTTAACCGCAATATTTTT  >  NC_000913/4277656‑4278134

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