Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,667,359 A→T 6.9% L163* (TTG→TAG)  hcaT ← putative 3‑phenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,667,3590AT 6.9% 157.8 / 3.5 58L163* (TTG→TAG) hcaTputative 3‑phenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base A (24/30);  new base T (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

ACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGCCAGTGCAT  >  NC_000913/2667211‑2667503
                                                                                                                                                    |                                                                                                                                                
acaaacGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCa                                                                                                                                                  >  2:1001208/1‑149 (MQ=255)
  aaacGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAAc                                                                                                                                                >  2:31281/1‑149 (MQ=255)
   aacGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGAGGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACa                                                                                                                                               <  1:1016003/149‑1 (MQ=255)
            cAGAAAGCGCCAGTTCTGGCGAACCAGCGCCCACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGa                                                                                                                                      >  1:238649/1‑149 (MQ=255)
                aaGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAc                                                                                                                                  <  1:26479/149‑1 (MQ=255)
                aaGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAc                                                                                                                                  <  1:1506297/149‑1 (MQ=255)
                 cgcgcCACGTCTGGCTAACCAGCGCCTACCAGGCAGACCACCCGGTGCTCTCCTTCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAcc                                                                                                                                 <  2:999783/148‑1 (MQ=255)
                 aGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCATCGTCAACAGCGCGAGGATCAcc                                                                                                                                 <  1:61182/149‑1 (MQ=255)
                      cAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATa                                                                                                                            >  2:1066325/1‑149 (MQ=255)
                      cAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATa                                                                                                                            <  2:515596/149‑1 (MQ=255)
                               cGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACAt                                                                                                                   >  1:694291/1‑149 (MQ=255)
                                    cAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGa                                                                                                              >  1:931255/1‑149 (MQ=255)
                                          cAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCTACAGCGCGAGGATCAccc                                                                                                                                <  2:1047/125‑1 (MQ=255)
                                          cAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCTACAGCGCGAGGATCAccc                                                                                                                                >  1:846912/1‑125 (MQ=255)
                                          cAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCTACAGCGCGAGGATCAccc                                                                                                                                >  1:1047/1‑125 (MQ=255)
                                          cAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCTACAGCGCGAGGATCAccc                                                                                                                                <  2:846912/125‑1 (MQ=255)
                                          cAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGtt                                                                                                        <  1:55217/149‑1 (MQ=255)
                                             ccACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCCCCCGATAATCAAACATAGTGACCAGTTTGc                                                                                                     >  1:1385771/1‑149 (MQ=255)
                                             ccACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGc                                                                                                     <  1:1001208/149‑1 (MQ=255)
                                             ccACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGc                                                                                                     <  1:1066325/149‑1 (MQ=255)
                                                    gACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAg                                                                                              >  1:851551/1‑149 (MQ=255)
                                                             gTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGcc                                                                                     >  1:1570375/1‑149 (MQ=255)
                                                                ctctCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAAt                                                                                  <  2:238649/149‑1 (MQ=255)
                                                                  ctcCTGCTGGCGACTTGCCCCTTGTGGCGGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGa                                                                                <  2:371067/149‑1 (MQ=255)
                                                                      tgctgGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACaaa                                                                            >  2:875870/1‑149 (MQ=255)
                                                                            cGACTTGCCCCTTGTGGCTGATTCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCACGTCAGCGCCGAGCCAATGACAAACGCCAc                                                                      <  2:1116438/149‑1 (MQ=255)
                                                                                 tGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGc                                                                 <  1:1263731/149‑1 (MQ=255)
                                                                                    cccTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGcccc                                                              >  2:953562/1‑149 (MQ=255)
                                                                                      cTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCcac                                                            >  1:512161/1‑149 (MQ=255)
                                                                                           ggCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGCGCCAATGACAAACGCCACCGAGCCCCACAGTCg                                                       >  1:795814/1‑149 (MQ=255)
                                                                                            gCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCCCCGAGCCCCACAGTCGc                                                      >  1:1499711/1‑149 (MQ=255)
                                                                                             cTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCa                                                     >  2:1278274/1‑149 (MQ=255)
                                                                                              tGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCAc                                                    <  2:1528096/149‑1 (MQ=255)
                                                                                               gAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACt                                                   <  2:339290/149‑1 (MQ=255)
                                                                                                aaTCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCCATAATCAAACAAAGTGACCACTTTCCCCGTCAGCTCCGACCCAATGACAAACGCCACCGAGCCCCACAGTCGCCCtt                                                  >  1:1430230/1‑149 (MQ=255)
                                                                                                 aTCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCa                                                                                                           >  2:1056953/1‑91 (MQ=255)
                                                                                                 aTCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCa                                                                                                           <  1:1056953/91‑1 (MQ=255)
                                                                                                  tcgtcgGACGGATGAGACAGCCGAGCAGCATGGATGCCACGCACAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTg                                                <  2:1018386/149‑1 (MQ=255)
                                                                                                            gaTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGTATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCaa                                      <  2:927170/149‑1 (MQ=255)
                                                                                                             aTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATTATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCATTGACAAACGCCACCGAGCCCCACAGTCGCACTATGCCGTAATCAAg                                     >  2:281298/1‑149 (MQ=255)
                                                                                                                 gaAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCggg                                 >  1:664021/1‑149 (MQ=255)
                                                                                                                  aaaGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGa                                >  2:1467775/1‑149 (MQ=255)
                                                                                                                           cagcaTGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCtttt                       <  2:1499007/149‑1 (MQ=255)
                                                                                                                                tGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCAc                  <  1:1403056/149‑1 (MQ=255)
                                                                                                                                tGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCAc                  <  1:875870/149‑1 (MQ=255)
                                                                                                                                tGGATGCCACGCCCAACGTCAAAAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCAc                  <  1:1063345/149‑1 (MQ=255)
                                                                                                                                   aTGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTa               >  1:525826/1‑149 (MQ=255)
                                                                                                                                   aTGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTa               >  1:419638/1‑149 (MQ=255)
                                                                                                                                      ccACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTc            >  2:439454/1‑149 (MQ=255)
                                                                                                                                      ccACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTc            <  2:694291/149‑1 (MQ=255)
                                                                                                                                        aCGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGc          >  1:809047/1‑149 (MQ=255)
                                                                                                                                          gCCCAACGTGAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACATACGCCACCGAGCCACTCAGTCGCTCTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTTTTCGCCa        <  1:578615/149‑1 (MQ=255)
                                                                                                                                           cccAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGCCAg       <  2:525826/149‑1 (MQ=255)
                                                                                                                                           cccAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGCCAg       <  2:664021/149‑1 (MQ=255)
                                                                                                                                            ccATCGGCAACATCGCGAGGATCACACGATAATCAAACATAGTGACCAGTTTGCCCGTAATCTCTGAGACAATGACAAACGCCACCGAGCCGCACAGTCGCACTTTGCCGTAATCAAGCGGGAATTGCTTTTGCCACGTATTCGCCAGt      <  2:259602/149‑1 (MQ=255)
                                                                                                                                            ccAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATATTGACCAGTTTGCCCGTCAGCGCCTAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTTATCAAGCGGGAATTGCTTTTGCCACGTATTCGCCAGt      <  2:358364/149‑1 (MQ=255)
                                                                                                                                            ccAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACACACGCCCCCGAGCCCCACAGTCGCACTTTGCGGTAATCAAGCGGGAACGGCTTTTGCCACTTTGTCGCCAGt      >  2:314896/1‑149 (MQ=255)
                                                                                                                                             cAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGGCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGCCAGTg     <  2:851551/149‑1 (MQ=255)
                                                                                                                                                cGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGCCAGTGCAt  <  1:31281/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCACCCGATAATCAAACATAGTGACCAGTTTGCCCGTCAGCGCCGAGCCAATGACAAACGCCACCGAGCCCCACAGTCGCACTTTGCCGTAATCAAGCGGGAACTGCTTTTGCCACGTATTCGCCAGTGCAT  >  NC_000913/2667211‑2667503

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