Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 557,870 (T)5→6 intergenic (+421/‑3) yddT → / → ydzN putative exported phage protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome557,8651.T93.5% 208.8 / 6.6 62intergenic (+416/‑8)yddT/ydzNputative exported phage protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base T (30/28);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTT  >  BS168_reference_genome/557729‑557997
                                                                                                                                          |                                                                                                                                      
tgAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTATTTAATTAAAAAGGAGtt                                                                                                                                      >  8:42474/1‑139 (MQ=255)
tgAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGtt                                                                                                                                      >  2:17319/1‑139 (MQ=255)
tgAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGtt                                                                                                                                      >  6:50759/1‑139 (MQ=255)
tgAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGtt                                                                                                                                      >  4:223963/1‑139 (MQ=255)
      tAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGttttttaa                                                                                                                               <  1:324788/139‑4 (MQ=255)
              tGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATATTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGtt                                                                                                                        <  2:323183/138‑1 (MQ=255)
                tAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTg                                                                                                                       >  8:114125/1‑137 (MQ=255)
                    aTCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTAATTGCGTTGGAGTaaa                                                                                                                <  6:94682/139‑1 (MQ=255)
                     tCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTaaa                                                                                                                <  2:47586/139‑1 (MQ=255)
                            aGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCa                                                                                                         <  4:299746/139‑1 (MQ=255)
                                 ggTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtt                                                                                                    <  4:234678/139‑1 (MQ=255)
                                   tttAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttt                                                                                                  >  1:204003/1‑139 (MQ=255)
                                   tttAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttt                                                                                                  <  7:114125/139‑1 (MQ=255)
                                     tAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGtttttt                                                                                                >  4:198555/1‑139 (MQ=255)
                                        gATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGta                                                                                             <  6:145698/139‑1 (MQ=255)
                                            aGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                          >  4:48580/1‑138 (MQ=255)
                                            aGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                          >  2:1087/1‑138 (MQ=255)
                                            aGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                          <  5:201282/138‑1 (MQ=255)
                                             gATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTAtag                                                                                          >  8:102636/1‑137 (MQ=255)
                                                  tGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCAtt                                                                                   >  5:282067/1‑139 (MQ=255)
                                                      aaaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGa                                                                               <  4:27177/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              <  5:200504/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              <  2:92585/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              <  2:87635/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              <  2:341399/139‑1 (MQ=255)
                                                       aaTTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              <  6:282067/139‑1 (MQ=255)
                                                              tatatTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatata                                                                       <  8:316452/139‑1 (MQ=255)
                                                                  ttttCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATgg                                                                    <  6:292621/138‑1 (MQ=255)
                                                                   tttCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGtt                                                                  >  5:192293/1‑139 (MQ=255)
                                                                                 tGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAata                                                                         >  1:341367/1‑118 (MQ=255)
                                                                                 tGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAata                                                                         <  2:341367/118‑1 (MQ=255)
                                                                                      cTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAAcc                                                >  1:148142/1‑138 (MQ=255)
                                                                                       tataATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCaa                                              <  7:271070/139‑1 (MQ=255)
                                                                                       tataATAAATGTATACAGATTGAATGGAGATGTAATTAATTAAAAATTATTTTTTTAATTGAGTTTCAGTAAATGGTTCAATTTTTTTTGTATAGTATAATTGGTATCAATATATGTTTATAAATTATTTACTAAcaa                                               >  2:315636/1‑136 (MQ=255)
                                                                                         taataaATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCaa                                              <  7:102636/137‑1 (MQ=255)
                                                                                                 gTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTCTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGta                                    >  8:192690/1‑139 (MQ=255)
                                                                                                    taCAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGtacta                                 <  2:204003/138‑1 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatata                                                                       <  4:82100/96‑1 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatata                                                                       >  3:82100/1‑96 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGt                            >  8:173219/1‑139 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGt                            >  7:71046/1‑139 (MQ=255)
                                                                                                         aTTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGAAGCATTGGGATCAATATATGGTTATAAATTATTTACTACCCAAGAAGTCAGTACTAGTCGt                            >  6:184415/1‑139 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTTAAAAGGAGTTTTTGAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTAAAAATTATTTACTAAACAAGAAGTCAGTACTAGTCGtt                           >  3:330274/1‑138 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAGCCAAGAAGTCAGAACTACTCGtt                           >  3:228767/1‑138 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                         >  8:248468/1‑124 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                         <  7:248468/124‑1 (MQ=255)
                                                                                                           tGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGtt                           >  1:147278/1‑138 (MQ=255)
                                                                                                            gTATGGAGATGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGttt                          >  6:170901/1‑138 (MQ=255)
                                                                                                                  agaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTtaataa                     >  2:240846/1‑137 (MQ=255)
                                                                                                                  agaTGTAATTAATTAAAAAGGAGTTATTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAACTTATTTACTCACCAAGATGTCAGTACGACTCGTTtagtaa                     >  4:345988/1‑137 (MQ=255)
                                                                                                                   gaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                         <  4:265101/116‑1 (MQ=255)
                                                                                                                   gaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGt                                         >  3:265101/1‑116 (MQ=255)
                                                                                                                   gaTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTtaataa                     >  6:104816/1‑136 (MQ=255)
                                                                                                                    aTGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGc                 <  1:1087/139‑1 (MQ=255)
                                                                                                                     tGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatata                                                                       >  4:128776/1‑84 (MQ=255)
                                                                                                                     tGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAatata                                                                       <  3:128776/84‑1 (MQ=255)
                                                                                                                     tGTAATTAATTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGc                 <  5:104816/138‑1 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              >  7:183087/1‑75 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAt                                                                              <  8:183087/75‑1 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc              >  5:68868/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc              >  6:3757/1‑139 (MQ=255)
                                                                                                                       taattaatTAAAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGAGCATTATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTc              >  1:276023/1‑139 (MQ=255)
                                                                                                                            aattaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTaa         >  5:138073/1‑139 (MQ=255)
                                                                                                                             attaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCaa                                              >  6:252891/1‑101 (MQ=255)
                                                                                                                             attaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCaa                                              <  5:252891/101‑1 (MQ=255)
                                                                                                                              ttaaAAAGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAAtt       >  1:153127/1‑139 (MQ=255)
                                                                                                                                   aaGGAGTTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTCTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAAttgtgtt  <  7:192690/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                      
TGAAAGTAATGATTTGTAGGATCAGTTAAGAAAGGTTTAAGATTAGATATTGAAAAATTTTATATATTTTCAATAAAATACTGGTACTATAATAAATGTATACAGATTGTATGGAGATGTAATTAATTAAAAAGGAGTTTTTAATTGCGTTGGAGTAAATGGTTCAATGTTTTTTGTATAGTAGCATTGGGATCAATATATGGTTATAAATTATTTACTAACCAAGAAGTCAGTACTACTCGTTTAATAATAGCTTCTGTAATTGTGTT  >  BS168_reference_genome/557729‑557997

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