Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,075,454 G→A *262* (TAG→TAA yihW → putative transcriptional regulator for sulphoquinovose utilization

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,075,4540GA100.0% 246.0 / NA 75*262* (TAG→TAAyihWputative transcriptional regulator for sulphoquinovose utilization
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (39/36);  total (39/36)

ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTG  >  NC_000913/4075318‑4075586
                                                                                                                                        |                                                                                                                                    
aTGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAt                                                                                                                                    >  1:385134/1‑139 (MQ=255)
    ttCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAAttttt                                                                                                                                >  2:1037296/1‑139 (MQ=255)
            cggTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGt                                                                                                                        >  2:1072976/1‑139 (MQ=255)
                  aaaTTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  <  1:1320487/139‑1 (MQ=255)
                                    aGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                <  2:467283/139‑1 (MQ=255)
                                     ggTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcaca                                                                                               <  1:1288732/139‑1 (MQ=255)
                                         aCTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCca                                                                                                 >  1:812280/1‑133 (MQ=255)
                                         aCTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCca                                                                                                 <  2:812280/133‑1 (MQ=255)
                                          cTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGcc                                                                                          >  2:359974/1‑139 (MQ=255)
                                                  tttACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  >  1:1154559/1‑107 (MQ=255)
                                                  tttACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  <  2:1154559/107‑1 (MQ=255)
                                                    tACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  1:260772/1‑93 (MQ=255)
                                                    tACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  2:260772/93‑1 (MQ=255)
                                                        gacgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           >  2:798396/1‑124 (MQ=255)
                                                        gacgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           <  1:798396/124‑1 (MQ=255)
                                                         acgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTcc                                                                           >  1:1176929/1‑139 (MQ=255)
                                                         acgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTcc                                                                           >  1:1455973/1‑139 (MQ=255)
                                                          cgaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTccc                                                                          >  2:1440358/1‑139 (MQ=255)
                                                           gaGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCg                                                                         <  2:1484817/139‑1 (MQ=255)
                                                            aGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGc                                                                        >  1:1473786/1‑139 (MQ=255)
                                                               tGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTAc                                                                     <  1:1072976/139‑1 (MQ=255)
                                                                 ccgcccgcTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAg                                                                   <  2:1473786/139‑1 (MQ=255)
                                                                    cccgcTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTa                                                                <  2:1176929/139‑1 (MQ=255)
                                                                      cgcTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCCGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAAt                                                              <  1:1037296/139‑1 (MQ=255)
                                                                                     aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaa                                               >  2:188032/1‑139 (MQ=255)
                                                                                     aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaa                                               >  2:282503/1‑139 (MQ=255)
                                                                                     aCGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaa                                               >  2:876382/1‑139 (MQ=255)
                                                                                          tACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAaggag                                          >  1:746676/1‑139 (MQ=255)
                                                                                            cAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAagga                                           <  1:279433/136‑1 (MQ=255)
                                                                                            cAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAagga                                           >  2:279433/1‑136 (MQ=255)
                                                                                                 cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaaa                                              >  2:1045490/1‑128 (MQ=255)
                                                                                                 cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGaaa                                              <  1:1045490/128‑1 (MQ=255)
                                                                                                 cAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCaaaaa                                   >  2:984611/1‑139 (MQ=255)
                                                                                                   gCCAAATTGTAATCATTCTTCCCCAGGAAGACGCTTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                                 <  1:984611/139‑1 (MQ=255)
                                                                                                   gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTaa                                                                                                                     <  1:1428525/55‑1 (MQ=255)
                                                                                                   gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTaa                                                                                                                     >  2:1428525/1‑55 (MQ=255)
                                                                                                   gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                   >  1:1177482/1‑57 (MQ=255)
                                                                                                   gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                   <  2:1177482/57‑1 (MQ=255)
                                                                                                   gCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATg                                 <  2:793586/139‑1 (MQ=255)
                                                                                                    ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                   >  1:1150541/1‑56 (MQ=255)
                                                                                                    ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAAcc                                                                                                                   <  2:1150541/56‑1 (MQ=255)
                                                                                                    ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           >  2:2124/1‑80 (MQ=255)
                                                                                                    ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           <  1:2124/80‑1 (MQ=255)
                                                                                                    ccAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGc                                >  1:43122/1‑139 (MQ=255)
                                                                                                     cAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              <  2:1252352/44‑1 (MQ=255)
                                                                                                     cAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTtgt                                                                                                                              >  1:1252352/1‑44 (MQ=255)
                                                                                                       aaTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGa                                                <  1:844825/120‑1 (MQ=255)
                                                                                                       aaTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGa                                                >  2:844825/1‑120 (MQ=255)
                                                                                                        aTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCta                            <  1:876382/139‑1 (MQ=255)
                                                                                                        aTTGAAATCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCta                            <  1:282503/139‑1 (MQ=255)
                                                                                                              aTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  >  1:440623/1‑47 (MQ=255)
                                                                                                              aTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  <  2:1124302/47‑1 (MQ=255)
                                                                                                              aTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  <  2:440623/47‑1 (MQ=255)
                                                                                                              aTCATTCTTCCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCt                                                                                                                  >  1:1124302/1‑47 (MQ=255)
                                                                                                                  ttcttcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           >  1:10732/1‑66 (MQ=255)
                                                                                                                  ttcttcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGc                                                                                           <  2:10732/66‑1 (MQ=255)
                                                                                                                  ttcttcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGAt                  >  1:904458/1‑139 (MQ=255)
                                                                                                                    cttcCCCAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGAttt                <  1:359974/139‑1 (MQ=255)
                                                                                                                       ccccAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCaaa                                                                                       >  1:107620/1‑65 (MQ=255)
                                                                                                                       ccccAGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCaaa                                                                                       <  2:107620/65‑1 (MQ=255)
                                                                                                                           aGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGt                                                                  >  1:561942/1‑82 (MQ=255)
                                                                                                                           aGGAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGt                                                                  <  2:561942/82‑1 (MQ=255)
                                                                                                                             gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                <  1:535415/50‑1 (MQ=255)
                                                                                                                             gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                >  2:535415/1‑50 (MQ=255)
                                                                                                                             gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGct                               <  1:777342/115‑1 (MQ=255)
                                                                                                                             gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGct                               >  2:777342/1‑115 (MQ=255)
                                                                                                                             gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAg              <  2:230858/132‑1 (MQ=255)
                                                                                                                             gAAGACGCGTAAATTTTTGTGACCGTTAACCTGGCTCCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAg              >  1:230858/1‑132 (MQ=255)
                                                                                                                               aGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                <  2:658702/48‑1 (MQ=255)
                                                                                                                               aGACGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCcac                                                                                                >  1:658702/1‑48 (MQ=255)
                                                                                                                                 aCGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGAt                  <  2:1076924/124‑1 (MQ=255)
                                                                                                                                 aCGCGTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGAt                  >  1:1076924/1‑124 (MQ=255)
                                                                                                                                  cgcgTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTg  >  2:357258/1‑139 (MQ=255)
                                                                                                                                     gTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTtgt                                                       <  2:39884/83‑1 (MQ=255)
                                                                                                                                     gTAAATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTtgt                                                       >  1:39884/1‑83 (MQ=255)
                                                                                                                                        |                                                                                                                                    
ATGCTTCGGCGGCGGTTGAAATTGGTAACGTGGCACAGGTCACTGCGCTCTTTACCGACGAGCTGCCGCCCGCTGCGCTAAAATCACGCTTACAAGACAGCCAAATTGAAATCATTCTTCCCCAGGAAGACGCGTAGATTTTTGTGACCGTTAACCTGGCTTCATACCTTGCCACATAGCCAAACCCATCCTTTCCCGCTACAGTTAATTTCTTGTGGCGCGAAAGGAGGCAAAAATGCTCTATATCTTTGATTTAGGTAATGTGATTG  >  NC_000913/4075318‑4075586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: