Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,683,627 T→C intergenic (‑180/+3) dctA ← / ← pdeK C4 dicarboxylate/orotate:H(+) symporter/putative c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,683,6270TC100.0% 219.0 / NA 59intergenic (‑180/+3)dctA/pdeKC4 dicarboxylate/orotate:H(+) symporter/putative c‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/34);  total (25/34)

TACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAA  >  NC_000913/3683492‑3683763
                                                                                                                                       |                                                                                                                                        
tACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCt                                                                                                                                      <  2:193381/140‑1 (MQ=255)
tACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCt                                                                                                                                      <  2:1156685/140‑1 (MQ=255)
 aCCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTa                                                                                                                                     >  2:467022/1‑140 (MQ=255)
 aCCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTa                                                                                                                                     >  1:705803/1‑140 (MQ=255)
 aCCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTa                                                                                                                                     >  1:632571/1‑140 (MQ=255)
  ccTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTAc                                                                                                                                    >  1:316636/1‑140 (MQ=255)
    tGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACtt                                                                                                                                  <  2:203618/140‑1 (MQ=255)
    tGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACtt                                                                                                                                  <  1:467022/140‑1 (MQ=255)
     gCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGTATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACttt                                                                                                                                 >  2:1151755/1‑140 (MQ=255)
     gCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACttt                                                                                                                                 >  2:973070/1‑140 (MQ=255)
     gCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACttt                                                                                                                                 <  2:511481/140‑1 (MQ=255)
      cgtccTTAAAATAACACCCAGCCAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACtttt                                                                                                                                >  2:983152/4‑140 (MQ=255)
           ttAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTcc                                                                                                                           <  2:307062/140‑1 (MQ=255)
            tAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCt                                                                                                                                      <  2:859398/128‑1 (MQ=255)
            tAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCt                                                                                                                                      >  1:859398/1‑128 (MQ=255)
            tAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCa                                                                                                                          >  1:269333/1‑140 (MQ=255)
                 tAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTaa                                                                                                                     <  2:1069574/140‑1 (MQ=255)
                   acacCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAAct                                                                                                                   <  2:466222/140‑1 (MQ=255)
                         aGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTc                                                                                                             >  2:64973/1‑140 (MQ=255)
                          gACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            >  2:851120/1‑140 (MQ=255)
                          gACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCg                                                                                                            <  1:1151755/140‑1 (MQ=255)
                            caacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          <  2:953246/140‑1 (MQ=255)
                            caacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGaa                                                                                                          >  2:590784/1‑140 (MQ=255)
                              acaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        <  2:316636/140‑1 (MQ=255)
                              acaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        <  1:851120/140‑1 (MQ=255)
                              acaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGa                                                                                                        <  1:990841/140‑1 (MQ=255)
                                aTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt                                                                                                      <  2:269333/140‑1 (MQ=255)
                                 tCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGAttt                                                                                                     >  1:559473/1‑140 (MQ=255)
                                    caGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa                                                                                                  <  2:632571/140‑1 (MQ=255)
                                     aGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAAt                                                                                                 >  2:290836/1‑140 (MQ=255)
                                                tGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGc                                                                                      <  1:1067717/140‑1 (MQ=255)
                                                      aTAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCaa                                                                                >  2:7860/1‑140 (MQ=255)
                                                        aaTTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCaaaa                                                                              >  1:700585/1‑140 (MQ=255)
                                                           tGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                           <  1:590784/140‑1 (MQ=255)
                                                           tGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                           <  1:782172/140‑1 (MQ=255)
                                                           tGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                           <  2:700585/140‑1 (MQ=255)
                                                           tGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAgg                                                                           <  1:1104784/140‑1 (MQ=255)
                                                                  aaaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCt                                                                    >  2:874269/1‑140 (MQ=255)
                                                                   aaaGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTg                                                                   <  1:874269/140‑1 (MQ=255)
                                                                     aGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTggg                                                                 <  1:973070/140‑1 (MQ=255)
                                                                            aaaGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATAcc                                                          <  2:49305/140‑1 (MQ=255)
                                                                                tttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACg                                                      >  2:1039656/1‑140 (MQ=255)
                                                                                 ttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGc                                                     >  2:920450/1‑140 (MQ=255)
                                                                                 ttGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGc                                                     <  2:598141/140‑1 (MQ=255)
                                                                                               acaAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                       <  2:243843/140‑1 (MQ=255)
                                                                                               acaAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                       <  1:1039656/140‑1 (MQ=255)
                                                                                                  aaCCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGccagccag                                    <  2:559473/140‑1 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                       <  1:266796/126‑1 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTgccagc                                       >  2:266796/1‑126 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGTCAGTCGCGTTGTGc                         >  2:777219/1‑140 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtg                            <  1:955929/137‑1 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTtg                            >  2:955929/1‑137 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGc                         >  2:344857/1‑140 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGc                         >  2:189513/1‑140 (MQ=255)
                                                                                                             tttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGc                         >  1:1205602/1‑140 (MQ=255)
                                                                                                              ttAAAGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCt                        <  1:613434/140‑1 (MQ=255)
                                                                                                                 aaGTTTTGTAATCAGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCa                     <  1:718052/140‑1 (MQ=255)
                                                                                                                              aGTTTGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCtt        >  2:674393/1‑140 (MQ=255)
                                                                                                                                tttGGGGCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCg      <  2:1205602/140‑1 (MQ=255)
                                                                                                                                    gggCAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCaa  <  1:7860/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAGTGGGCGAGCAAAAAGGAAGCCCTGGGCAATACCAACGCCCGCTTTTGCCAGCCAGTCGCGTTGTGCTTCAGTCTCCACGCCTTCGGCAA  >  NC_000913/3683492‑3683763

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: