Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,650,018 T→A 9.1% intergenic (+1495/+196) PP_3214 → / ← PP_3216 hypothetical protein/cysteine dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,650,0180TA 9.1% 150.5 / 5.6 55intergenic (+1495/+196)PP_3214/PP_3216hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base A (2/3);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

GGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGG  >  NC_002947/3649880‑3650165
                                                                                                                                            |                                                                                                                                                   
ggCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGa                                                                                                                                            <  5:175172/148‑1 (MQ=255)
   aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGAggg                                                                                                                                         <  1:286024/148‑3 (MQ=255)
   aCGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCTACGATGACGGTGAggg                                                                                                                                         >  2:286024/1‑146 (MQ=255)
            tCACGAAACAACTCGCCGCCCGCAGCGCACCAGCCGGCAAGGCCGCTATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCccg                                                                                                                               <  6:210983/149‑1 (MQ=255)
                gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                           >  3:298449/1‑149 (MQ=255)
                gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                           >  2:28997/1‑149 (MQ=255)
                 aaaCAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCggg                                                                                                                          >  7:225120/1‑149 (MQ=255)
                     aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc                                                                                                                      >  3:179349/1‑149 (MQ=255)
                     aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc                                                                                                                      >  3:32290/1‑149 (MQ=255)
                       cTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgc                                                                                                                    >  4:232681/1‑149 (MQ=255)
                        tCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgcg                                                                                                                   >  3:260309/1‑149 (MQ=255)
                         cgccgcCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg                                                                                                                  >  1:279982/1‑149 (MQ=255)
                          gccgccCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAGACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGCGTGCGGGGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg                                                                                                                  >  1:235913/1‑149 (MQ=255)
                           ccgcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGc                                                                                                                <  1:28997/149‑1 (MQ=255)
                           ccgcccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCACCGGGGAACGCGGGc                                                                                                                >  3:89724/1‑149 (MQ=255)
                                   aGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGAttt                                                                                                        <  4:32290/149‑1 (MQ=255)
                                    gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAATGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTTCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTc                                                                                                       <  6:86345/149‑1 (MQ=255)
                                     cgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTcc                                                                                                      <  8:191895/149‑1 (MQ=255)
                                        gcCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGc                                                                                                   >  2:36064/1‑149 (MQ=255)
                                            gCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCa                                                                                               <  4:26848/149‑1 (MQ=255)
                                               ggCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGtt                                                                                            >  2:183131/1‑149 (MQ=255)
                                                gCAAGGCCGCCATCAAGCAAGGCA‑ACATTGAAGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGAGGTGTCCCGCCGCGGAACGCGGGCTAAGATTTCCAGCTCCAGttt                                                                                           >  4:30051/1‑149 (MQ=255)
                                                 cAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCACCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                          <  4:89724/149‑1 (MQ=255)
                                                 cAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGCCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTc                                                                                          <  5:278516/149‑1 (MQ=255)
                                                     gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACg                                                                                                                                                >  6:145931/1‑91 (MQ=255)
                                                     gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACg                                                                                                                                                <  5:145931/91‑1 (MQ=255)
                                                     gccgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACa                                                                                      <  3:232681/149‑1 (MQ=255)
                                                     cgcgccATCAAGCAAGGCT‑ACATTGACTTAGCCCCAGGCCAGTAAACGCTGAGCGGCGTCTGCGGCGA‑GCCCCTCACCGGCATCGTAGAGTGTGAGTGCGGTGTCCCGCGGGGGAACGCGGGCGAAGATTTCCAGCTGAAGTTTCgg                                                                                        >  8:216217/3‑146 (MQ=255)
                                                       cgccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAgg                                                                                    >  1:54475/1‑149 (MQ=255)
                                                        gccATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAggg                                                                                   <  6:86639/149‑1 (MQ=255)
                                                           aTCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTa                                                                                >  8:252365/1‑149 (MQ=255)
                                                                 caagGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTgg                                                                          <  1:183131/149‑1 (MQ=255)
                                                                    ggCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGGGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCcg                                                                       >  6:21361/1‑149 (MQ=255)
                                                                         aCATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGaa                                                                   <  1:36064/149‑1 (MQ=255)
                                                                               aCGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCg                                                              >  7:54203/1‑148 (MQ=255)
                                                                                    gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg                                                         <  3:74816/148‑1 (MQ=255)
                                                                                      ccAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGcc                                                      >  2:312897/1‑149 (MQ=255)
                                                                                        aaGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTg                                                    <  8:54203/149‑1 (MQ=255)
                                                                                         aGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCACCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTgg                                                   <  5:309010/149‑1 (MQ=255)
                                                                                           gCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc                                                 <  7:252365/149‑1 (MQ=255)
                                                                                           gCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc                                                 <  2:279982/149‑1 (MQ=255)
                                                                                               gTAAACGCTGCGCGGCGACTGCGGCGACGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGa                                               >  1:14704/1‑148 (MQ=255)
                                                                                                       tGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGTGAACAGCGGATGGGCCTGGGCGAAAGGGtcttc                                     >  4:187597/1‑149 (MQ=255)
                                                                                                         cgcgGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGtcttc                                     <  1:164535/147‑1 (MQ=255)
                                                                                                              cgACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGt                              <  2:54475/149‑1 (MQ=255)
                                                                                                                  ttcggcgA‑GCCCCTCACCGTCATCGTAGAAGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGt                          <  3:30051/147‑1 (MQ=255)
                                                                                                                     ggcgA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGTGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGa                       <  3:187597/149‑1 (MQ=255)
                                                                                                                      gcgA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGAt                      <  8:225120/149‑1 (MQ=255)
                                                                                                                        gA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGACGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAATGGTCTTCTTCGCGTGAGTCGAt                      >  2:304227/1‑147 (MQ=255)
                                                                                                                           gCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                  <  2:235913/149‑1 (MQ=255)
                                                                                                                              ccTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGcc               <  1:312897/149‑1 (MQ=255)
                                                                                                                              acTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa                             >  6:193638/2‑135 (MQ=255)
                                                                                                                              acTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa                             <  5:193638/134‑1 (MQ=255)
                                                                                                                                     gTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGctcct        >  7:101060/1‑149 (MQ=255)
                                                                                                                                           gTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGGCTTCTTCGCGTGCGTCGATCAGGGCCAGCTCCTCCTTgg  <  1:304227/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCA‑ACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGA‑GCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGG  >  NC_002947/3649880‑3650165

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