Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 909,319 G→A 7.0% intergenic (‑111/‑174) cra ← / → fruB catabolite DNA‑binding transcriptional regulator/PTS fructose transporter subunit EI/HPr/IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947909,3190GA 7.0% 164.8 / 2.7 57intergenic (‑111/‑174)cra/fruBcatabolite DNA‑binding transcriptional regulator/PTS fructose transporter subunit EI/HPr/IIA
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base A (2/2);  total (28/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 = 8.44e-01

ACACACCGGCCAGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTCACAACAATACTCCA  >  NC_002947/909172‑909459
                                                                                                                                                    |                                                                                                                                            
acacacCGGCCAGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGCCCATGGTCGCAATTGTCCGACAGGATgggg                                                                                                                                             >  6:261171/1‑149 (MQ=255)
    acCGGCCAGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCTGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCttt                                                                                                                                         >  8:78765/1‑149 (MQ=255)
    acCGGCCAGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCttt                                                                                                                                         >  3:148144/1‑149 (MQ=255)
    acCGGCCAGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCttt                                                                                                                                         >  7:226622/1‑149 (MQ=255)
           aGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGACGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTgg                                                                                                                                  >  2:84259/1‑149 (MQ=255)
             aCGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACt                                                                                                                                    >  2:65597/1‑145 (MQ=255)
             aCGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACt                                                                                                                                    <  1:65597/145‑1 (MQ=255)
               gggCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAg                                                                                                                              <  1:49224/149‑1 (MQ=255)
                       tCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGCTCCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATc                                                                                                                        <  1:84259/148‑1 (MQ=255)
                         gCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAg                                                                                                                     <  5:85682/148‑1 (MQ=255)
                                    ccGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGt                                                                                                         >  3:46788/1‑149 (MQ=255)
                                     cGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTc                                                                                                        <  7:78765/149‑1 (MQ=255)
                                     cGAATTTTCCTTGTTATTACCGGAGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCAGTc                                                                                                        >  8:278765/1‑149 (MQ=255)
                                        aTTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTg                                                                                                      <  4:284187/148‑1 (MQ=255)
                                         ttttCCTTGTTATTACCGGGGCTGCCTGAAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGt                                                                                                    >  1:194343/1‑149 (MQ=255)
                                              cTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGa                                                                                                <  2:194343/148‑1 (MQ=255)
                                                  ttattaCCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTaaa                                                                                           >  8:157212/1‑149 (MQ=255)
                                                   tattaCCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTaaa                                                                                           <  3:237762/148‑1 (MQ=255)
                                                   tattaCCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGAGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTaaa                                                                                           <  5:161476/148‑1 (MQ=255)
                                                   tattaCCGGAGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTaaa                                                                                           <  7:278765/148‑1 (MQ=255)
                                                    attaCCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTa                                                                                                                                      >  7:97320/1‑104 (MQ=255)
                                                    attaCCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTa                                                                                                                                      <  8:97320/104‑1 (MQ=255)
                                                    attaCCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAAc                                                                                          <  8:42500/148‑1 (MQ=255)
                                                              cTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGc                                                                               <  7:157212/149‑1 (MQ=255)
                                                                    gTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTCGTCAGATTAAACGTTTCAGCTGCCGAAtt                                                                         >  6:13713/1‑149 (MQ=255)
                                                                     tAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAAttt                                                                        >  7:72590/1‑149 (MQ=255)
                                                                        ggCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCt                                                                     >  6:72599/1‑149 (MQ=255)
                                                                            ggACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCt                                                                 <  8:226622/149‑1 (MQ=255)
                                                                                  ggTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTg                                                           <  4:148144/149‑1 (MQ=255)
                                                                                  ggTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTg                                                           <  8:72590/149‑1 (MQ=255)
                                                                                            cccGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCcgcg                                                  >  7:235570/1‑148 (MQ=255)
                                                                                             caggTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCgcgc                                                 >  6:298298/3‑148 (MQ=255)
                                                                                               ggTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAg                                              >  3:260705/1‑149 (MQ=255)
                                                                                                gTACAGGCATCACGCGCAGACGACCATTGTCGCAATAGTCCGACAGGATGGGGCTTTTTACTGGGCAGAATATCGAGTAACGTGGCCTTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTATCTGCCTCGGCTGCCACCCGCGCAAgg                                             >  5:39698/1‑149 (MQ=255)
                                                                                                   cAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGtt                                           <  4:260705/148‑1 (MQ=255)
                                                                                                    aGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGAt                                                                                               >  6:185789/1‑95 (MQ=255)
                                                                                                      gCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGc                                         <  5:52714/147‑1 (MQ=255)
                                                                                                      gCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGc                                         >  6:52714/1‑147 (MQ=255)
                                                                                                         tCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAAc                                    >  3:151832/1‑149 (MQ=255)
                                                                                                             gcgcAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGc                                >  4:117769/1‑149 (MQ=255)
                                                                                                               gcAGACGACCATTGTCGCAATTGTCCGACAGGATGGAGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGccc                              >  8:19166/1‑149 (MQ=255)
                                                                                                               gcAGACGACCATTGTCGCAATTGTCCGACAGGATGGAGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGccc                              >  8:14195/1‑149 (MQ=255)
                                                                                                                 aGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCtt                            >  3:316885/1‑149 (MQ=255)
                                                                                                                  gacgacCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCttt                           <  5:341065/149‑1 (MQ=255)
                                                                                                                  gacgacCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCttt                           <  5:240999/149‑1 (MQ=255)
                                                                                                                  gacgacCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCtt                            >  6:341065/1‑148 (MQ=255)
                                                                                                                   acgacCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTa                          <  5:298298/149‑1 (MQ=255)
                                                                                                                    cgacCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTAt                         <  3:117769/149‑1 (MQ=255)
                                                                                                                       acCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCTTCTTGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTAtgt                       <  5:13713/149‑1 (MQ=255)
                                                                                                                          aTTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGaa                    >  3:181128/1‑149 (MQ=255)
                                                                                                                           ttGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAAt                   >  5:246667/1‑149 (MQ=255)
                                                                                                                           ttGTCGCAATTGTCCGACAGCATTGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGGCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAAt                   <  4:283755/149‑1 (MQ=255)
                                                                                                                               cGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTca               >  3:155379/1‑149 (MQ=255)
                                                                                                                                gCAATTGTCCGACAGGATGGAGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTca               <  7:19166/148‑1 (MQ=255)
                                                                                                                                gCAATTGTCCGACAGGATGGAGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTca               <  7:14195/148‑1 (MQ=255)
                                                                                                                                      gTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTCAcaacaa         <  4:41753/148‑1 (MQ=255)
                                                                                                                                            cAGGATGGGGGTTTTTACTGGGCAGATTTTCGAGTTACTTGGCCGTCTGGTCTGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGTGTGCCACCGGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTCACAACAATACTCCa  <  6:39698/149‑1 (MQ=255)
                                                                                                                                            cAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTCACAACAATACTCCa  <  4:46788/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
ACACACCGGCCAGACGGGCGATATCGCTGAGTTTCACCGAATTTTCCTTGTTATTACCGGGGCTGCCTGTAAGGCCGGACCGGGTAGTTTGACCCGGTACAGGCATCACGCGCAGACGACCATTGTCGCAATTGTCCGACAGGATGGGGCTTTTTACTGGGCAGATTATCGAGTAACGTGGCCGTCTGGTCAGATTAAACGTTTCAGCTGCCGAATTTTCTGCCTCGGCTGCCACCCGCGCAAGGTTGCTGAACTGGCCCTTTATGTGAATTTCACAACAATACTCCA  >  NC_002947/909172‑909459

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