Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,750,060 C→T 7.1% intergenic (‑81/+48) typA ← / ← thiI ribosome associated GTPase/tRNA sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,750,0600CT 7.1% 163.2 / 2.2 56intergenic (‑81/+48)typA/thiIribosome associated GTPase/tRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base T (2/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

TCGACCAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAGT  >  NC_002947/5749912‑5750205
                                                                                                                                                    |                                                                                                                                                  
tCGACCAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCcggc                                                                                                                                                    >  5:332231/1‑149 (MQ=255)
  gACCAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCcggcgg                                                                                                                                                  >  2:161807/1‑149 (MQ=255)
  gACCAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCcggcgg                                                                                                                                                  >  3:118833/1‑149 (MQ=255)
   aCCAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATAAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGt                                                                                                                                                 >  6:235961/1‑149 (MQ=255)
     cAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCg                                                                                                                                               <  8:340539/149‑1 (MQ=255)
            gTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGc                                                                                                                                        >  2:136705/1‑149 (MQ=255)
                 aCCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATgg                                                                                                                                   >  8:245077/1‑149 (MQ=255)
                  ccATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATgg                                                                                                                                   <  5:260053/148‑1 (MQ=255)
                     tGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGaa                                                                                                                               <  5:67171/149‑1 (MQ=255)
                       gTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATa                                                                                                                             <  3:78521/149‑1 (MQ=255)
                       gTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATa                                                                                                                             <  4:118833/149‑1 (MQ=255)
                       gTCAACGTGGGCGATGATGGCGAAGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATa                                                                                                                             <  8:311253/149‑1 (MQ=255)
                        tCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATAc                                                                                                                            >  7:345908/1‑149 (MQ=255)
                                 gCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCAt                                                                                                                   <  7:245077/149‑1 (MQ=255)
                                  cGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGTCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCAGCCAt                                                                                                                   >  5:215727/1‑148 (MQ=255)
                                           cGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCgg                                                                                                         <  8:345908/149‑1 (MQ=255)
                                           cGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCCTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCgg                                                                                                         >  7:208358/1‑149 (MQ=255)
                                                    gCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAggg                                                                                                                                         >  8:184763/1‑108 (MQ=255)
                                                    gCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAggg                                                                                                                                         <  7:184763/108‑1 (MQ=255)
                                                    gCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATg                                                                                                >  1:91167/1‑149 (MQ=255)
                                                    gCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGCAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGaa                                                                                                                               >  8:214086/1‑118 (MQ=255)
                                                        aTTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGc                                                                                                        <  6:63790/137‑1 (MQ=255)
                                                        aTTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGc                                                                                                        >  5:63790/1‑137 (MQ=255)
                                                              gATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATa                                                                                      <  6:110545/149‑1 (MQ=255)
                                                                            ctcGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCCGCGGGCCCGGCGGTCGGGAGGGCGATGGCGATTACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTgg                                                                       >  6:292362/1‑149 (MQ=255)
                                                                              cGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGc                                                                      <  1:136705/148‑1 (MQ=255)
                                                                               gATCAGAGGTTTCGGGTGCCGGCCAGTGTATGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATTGCGAATACTGCCGCCATACATCCCCGGCGTCTTATTTCGGAAGATAAACACGCACATTg                                                                        <  6:215727/145‑1 (MQ=255)
                                                                               gATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCAt                                                                    >  5:134797/1‑149 (MQ=255)
                                                                                   agagGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTcc                                                                >  7:353701/1‑149 (MQ=255)
                                                                                    gagGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTccc                                                               >  2:288392/1‑149 (MQ=255)
                                                                                       gATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCa                                                            >  1:195744/1‑149 (MQ=255)
                                                                                        aTTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCAc                                                           >  8:325691/1‑149 (MQ=255)
                                                                                         ttCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACt                                                          >  1:9640/1‑149 (MQ=255)
                                                                                            ggTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGa                                                        <  8:208358/148‑1 (MQ=255)
                                                                                             gTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCTGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGc                                                      <  5:235961/149‑1 (MQ=255)
                                                                                                 gctcccAGTGTATGCGGCCAATAGGTACCGTTTAATGCTCAGCGGGCCCGGCGTTCGGGAGGGCGATTGCTAATACTGCAGCCATACATCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTACCTCACTGAGCAGAt                                                  <  8:110425/146‑1 (MQ=255)
                                                                                                  cgccCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATg                                                 <  7:100358/149‑1 (MQ=255)
                                                                                                      cAGTGTAGGCGGCCAATATGTACCGTTTAACGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGgtgt                                             <  4:338466/149‑1 (MQ=255)
                                                                                                             ggcggcCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTgtgcgtgc                                       <  1:161807/148‑1 (MQ=255)
                                                                                                                 gcCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACg                                  >  7:94102/1‑149 (MQ=255)
                                                                                                                    aaTATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCt                               >  4:299268/1‑149 (MQ=255)
                                                                                                                        tGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACCACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCtca                             <  1:288392/148‑1 (MQ=255)
                                                                                                                        tGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCtcatc                           <  8:94102/149‑1 (MQ=255)
                                                                                                                          tACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCtcatca                          >  6:115030/1‑148 (MQ=255)
                                                                                                                          tACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCtcatca                          >  8:29868/1‑148 (MQ=255)
                                                                                                                            acGTTTAAGGCTCAGCGGGCCCCGCGGTCGGGAGGGCGAGGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATcaca                        <  6:231144/147‑1 (MQ=255)
                                                                                                                             cGTTTAAGGCTCAGCGGGCCCGGTGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGgtgt                                             <  5:177814/126‑1 (MQ=255)
                                                                                                                             cGTTTAAGGCTCAGCGGGCCCGGTGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGgtgt                                             >  6:65551/1‑126 (MQ=255)
                                                                                                                             cGTTTAAGGCTCAGCGGGCCCGGTGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGgtgt                                             >  6:177814/1‑126 (MQ=255)
                                                                                                                             cGTTTAAGGCTCAGCGGGCCCGGTGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGgtgt                                             <  5:65551/126‑1 (MQ=255)
                                                                                                                                      cTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACAATGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGc              >  5:228357/1‑148 (MQ=255)
                                                                                                                                       tCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCgg                                                                                            >  7:160895/1‑70 (MQ=255)
                                                                                                                                       tCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCgg                                                                                            <  8:160895/70‑1 (MQ=255)
                                                                                                                                       tCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCt                                                              >  6:34224/1‑99 (MQ=255)
                                                                                                                                       tCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCt                                                              <  5:34224/99‑1 (MQ=255)
                                                                                                                                         aGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAAt           <  7:29868/148‑1 (MQ=255)
                                                                                                                                            gggCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAATACAg       >  4:234967/1‑149 (MQ=255)
                                                                                                                                                ccGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAg   >  7:102054/1‑149 (MQ=255)
                                                                                                                                                ccGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAg   >  2:96590/1‑149 (MQ=255)
                                                                                                                                                  ggcggTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAGt  >  2:368725/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGACCAGAGTGGTTTTACCATGGTCAACGTGGGCGATGATGGCGATGTTACGCAGATTTTCGATCACAAGTGTATCTCGATCAGAGGATTCGGTTGCCGCCCAGTGTAGGCGGCCAATATGTACCGTTTAAGGCTCAGCGGGCCCGGCGGTCGGGAGGGCGATGGCGAATACTGCCGCCATACAGCCCCGGCGTCTTATGTCGGACGATAAACACGCACATTGGCATGTCCCTCACTGAGCAGATGGTGTGCGTGCAAACGGCTCATCACACCCTTGTCGCAATACAGCAAGT  >  NC_002947/5749912‑5750205

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: