Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,663,209 G→T 55.0% intergenic (+11/‑11) U712_08550 → / → U712_08555 FlaA locus uncharacterized protein ylxG/Flagellar basal‑body rod protein flgG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,663,2090GT55.0% 2.0 / 74.4 58intergenic (+11/‑11)U712_08550/U712_08555FlaA locus uncharacterized protein ylxG/Flagellar basal‑body rod protein flgG
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base T (14/18);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 7.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

CAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCGTGTGACCTTTAAAGA  >  CP006881/1663074‑1663347
                                                                                                                                       |                                                                                                                                          
cAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCtggtgg                                                                                                                                          <  8:18759/138‑1 (MQ=255)
cAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTggggg                                                                                                                                          >  8:314366/1‑138 (MQ=255)
cAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTggggg                                                                                                                                          <  5:288271/138‑1 (MQ=255)
cAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         <  7:146825/139‑1 (MQ=255)
cAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGa                                                                                                                                         <  7:179965/139‑1 (MQ=255)
cAAAAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         <  5:174242/139‑1 (MQ=255)
  aCAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAat                                                                                                                                       <  8:275155/139‑1 (MQ=255)
     gATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAatata                                                                                                                                    <  3:185955/139‑1 (MQ=255)
      aTAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         >  7:200132/1‑133 (MQ=255)
      aTAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         <  8:200132/133‑1 (MQ=255)
               gCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTagg                                                                                                                          >  7:277992/1‑137 (MQ=255)
                cacaGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGt                                                                                                                         <  6:229947/139‑1 (MQ=255)
                cacaGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGt                                                                                                                         <  4:151550/139‑1 (MQ=255)
                      aaGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACt                                                                                                                    <  2:26272/138‑1 (MQ=255)
                         ctctGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTgggg                                                                                                                                           >  3:361194/1‑112 (MQ=255)
                         ctctGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTgggg                                                                                                                                           <  4:361194/112‑1 (MQ=255)
                         ctctGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         <  5:318907/114‑1 (MQ=255)
                         ctctGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         >  6:318907/1‑114 (MQ=255)
                         ctctGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTAt                                                                                                                >  2:312243/1‑139 (MQ=255)
                            tGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCt                                                                                                             >  8:115399/1‑139 (MQ=255)
                             gTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTg                                                                                                            >  3:190641/1‑139 (MQ=255)
                             gTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTg                                                                                                            >  5:364077/1‑139 (MQ=255)
                             gTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTg                                                                                                            >  2:110159/1‑139 (MQ=255)
                             gTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTg                                                                                                            >  7:301987/1‑139 (MQ=255)
                                      tttAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATTTATTATGTTACGTTCACTTTATTCTGGATGCAGCg                                                                                                   >  2:373932/1‑139 (MQ=255)
                                      tttAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCg                                                                                                   >  5:399332/1‑139 (MQ=255)
                                      tttAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCg                                                                                                   >  3:209038/1‑139 (MQ=255)
                                           aGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTAt                                                                                               >  5:184269/1‑138 (MQ=255)
                                              aaaTTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGaaa                                                                                           <  1:56715/139‑1 (MQ=255)
                                                attattATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         <  2:403225/91‑1 (MQ=255)
                                                attattATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGa                                                                                                                                         >  1:403225/1‑91 (MQ=255)
                                                attattATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGaaaa                                                                                          >  5:34478/1‑138 (MQ=255)
                                                           gTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTcaaaca                                                                               <  7:314366/138‑1 (MQ=255)
                                                            tccTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTcaaaca                                                                               >  7:185519/3‑137 (MQ=255)
                                                                      tgGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTa                                                                   <  6:364077/139‑1 (MQ=255)
                                                                         gACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCTTCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTGACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAAt                                                                 <  7:2198/138‑1 (MQ=255)
                                                                            cGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCgg                                                              >  8:105333/1‑138 (MQ=255)
                                                                                   aGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAacaaca                                                      >  2:291782/1‑139 (MQ=255)
                                                                                          cGAATGTCATGTCTGTTGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGcc                                                 <  8:301987/137‑1 (MQ=255)
                                                                                          cGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGcc                                                 <  7:105333/137‑1 (MQ=255)
                                                                                          cGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGcc                                                 <  8:185519/137‑1 (MQ=255)
                                                                                          cGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGcc                                                 <  7:115399/137‑1 (MQ=255)
                                                                                          cGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCaa                                               >  8:59050/1‑139 (MQ=255)
                                                                                          cGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCaa                                               >  8:244385/1‑139 (MQ=255)
                                                                                             atgtcatgtcTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGt                                            <  4:209038/139‑1 (MQ=255)
                                                                                               gtcatgtcTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTa                                           >  1:278732/1‑138 (MQ=255)
                                                                                                tcatgtcTGTGGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTaaa                                         <  2:278732/139‑1 (MQ=255)
                                                                                                  atgtcTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATa                                       <  6:34478/139‑1 (MQ=255)
                                                                                                    gtcTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACg                                                                     <  1:152133/107‑1 (MQ=255)
                                                                                                    gtcTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACg                                                                     >  2:152133/1‑107 (MQ=255)
                                                                                                         tgGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGt                                            <  6:372634/127‑1 (MQ=255)
                                                                                                         tgGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGt                                            >  5:372634/1‑127 (MQ=255)
                                                                                                         tgGGACAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGt                                   <  1:110159/136‑1 (MQ=255)
                                                                                                         tgGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGt                                   <  7:315594/136‑1 (MQ=255)
                                                                                                              cAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGcc                                                 <  5:269085/117‑1 (MQ=255)
                                                                                                              cAATCATCTAAATAAAAACATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGcc                                                 >  6:269085/1‑117 (MQ=255)
                                                                                                              cAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTACTCGGTAACAACATTGCCAACGTAAATACAGTAGGATTc                            >  4:13909/1‑138 (MQ=255)
                                                                                                                 tcatcTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGAt                              <  2:186892/133‑1 (MQ=255)
                                                                                                                 tcatcTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGAt                              >  1:186892/1‑133 (MQ=255)
                                                                                                                 tcatcTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCaaa                         >  4:123706/1‑138 (MQ=255)
                                                                                                                        aataaaAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAagag                    <  1:312243/136‑1 (MQ=255)
                                                                                                                            aaaaCATCTGGTGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCgtgt              <  6:184269/138‑1 (MQ=255)
                                                                                                                                 aTCTTGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCGTGTGACCtt        <  2:343506/139‑1 (MQ=255)
                                                                                                                                      ggggAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCGTGTGACCttt       <  7:59050/135‑1 (MQ=255)
                                                                                                                                      ggggAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCGTGTGACCttt       <  7:244385/135‑1 (MQ=255)
                                                                                                                                       gggAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCGTGTGACCTTTAAAGa  <  8:277992/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CAACAGATAAAACAGGCACAGTAAGCTCTGTTAAACATTTTAAAGGAAATTATTATCTCGTTCTTGATGATGGGACCGAGATCAGTCCTGCGAATGTCATGTCTGTGGGACAATCATCTAAATAAAAACATCTGGGGGAATATATTATGTTACGTTCACTTTATTCTGGAATCAGCGGTATGAAAAACTTTCAAACAAAGCTTGACGTAATCGGTAACAACATTGCCAACGTAAATACAGTAGGATTCAAAAAGAGCCGTGTGACCTTTAAAGA  >  CP006881/1663074‑1663347

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: