Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,098,442 A→C intergenic (‑141/+5) wzzB ← / ← ugd regulator of length of O‑antigen component of lipopolysaccharide chains/UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,098,4420AC100.0% 219.7 / NA 58intergenic (‑141/+5)wzzB/ugdregulator of length of O‑antigen component of lipopolysaccharide chains/UDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (37/21);  total (37/21)

CTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACAAGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGTTGAAGAATGA  >  NC_000913/2098295‑2098587
                                                                                                                                                   |                                                                                                                                                 
ctctCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGAt                                                                                                                                                >  1:801222/1‑151 (MQ=255)
  ctcATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGt                                                                                                                                              >  2:856597/1‑151 (MQ=255)
       catTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGt                                                                                                                                          <  2:705783/150‑1 (MQ=255)
       catTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTc                                                                                                                                         >  2:531549/1‑151 (MQ=255)
       catTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTc                                                                                                                                         >  1:913962/1‑151 (MQ=255)
       catTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTc                                                                                                                                         >  1:371984/1‑151 (MQ=255)
          tACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCt                                                                                                                                      >  1:529621/1‑151 (MQ=255)
          tACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGAGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCt                                                                                                                                      >  1:798295/1‑151 (MQ=255)
             ccTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGcc                                                                                                                                   >  2:647869/1‑151 (MQ=255)
                aaCTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCaaa                                                                                                                                <  2:429379/151‑1 (MQ=255)
                aaCTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCa                                                                                                                                  <  2:179972/149‑1 (MQ=255)
                 aCTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCaa                                                                                                                                 >  2:266741/1‑149 (MQ=255)
                          aaTTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCgg                                                                                                                      <  2:827480/151‑1 (MQ=255)
                                  aGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATAc                                                                                                              >  2:890755/1‑151 (MQ=255)
                                   gTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATAcc                                                                                                             >  2:461057/1‑151 (MQ=255)
                                         aTAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCt                                                                                                       >  1:759433/1‑151 (MQ=255)
                                               tAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCcac                                                                                                  <  2:877263/150‑1 (MQ=255)
                                                    aTAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCtt                                                                                            <  1:317463/151‑1 (MQ=255)
                                                      aaTCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTa                                                                                           >  2:502384/1‑150 (MQ=255)
                                                      aaTCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTa                                                                                           >  1:236633/1‑150 (MQ=255)
                                                       aTCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTa                                                                                           >  2:843770/1‑149 (MQ=255)
                                                       aTCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTa                                                                                           >  2:657946/1‑149 (MQ=255)
                                                       aTCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAg                                                                                         <  1:51212/151‑1 (MQ=255)
                                                        tCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGc                                                                                        >  1:375710/1‑151 (MQ=255)
                                                        tCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGc                                                                                        >  1:23423/1‑151 (MQ=255)
                                                           tCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGctct                                                                                     >  1:872460/1‑151 (MQ=255)
                                                               gAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTg                                                                                 <  2:592504/151‑1 (MQ=255)
                                                                   gCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCa                                                                              <  1:788221/150‑1 (MQ=255)
                                                                     tataGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTc                                                                           >  1:99790/1‑151 (MQ=255)
                                                                      ataGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCAtt                                                                            <  2:855854/149‑1 (MQ=255)
                                                                        aGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGt                                                                        >  2:862983/1‑151 (MQ=255)
                                                                          gATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGtt                                                                       >  2:274982/1‑150 (MQ=255)
                                                                           aTAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGtt                                                                       >  1:430711/1‑149 (MQ=255)
                                                                                    aaaaaTAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATg                                                             <  1:125748/150‑1 (MQ=255)
                                                                                    aaaaaTAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGa                                                            <  2:318062/151‑1 (MQ=255)
                                                                                       aataatGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGt                                                         <  1:484335/151‑1 (MQ=255)
                                                                                          aatGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCgg                                                      >  1:418748/1‑151 (MQ=255)
                                                                                          aatGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCg                                                       <  2:51519/150‑1 (MQ=255)
                                                                                              aaTTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCttg                                                  >  2:878506/1‑151 (MQ=255)
                                                                                                    tgtgTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGttt                                            >  1:18225/1‑151 (MQ=255)
                                                                                                     gtgtAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTg                                           <  2:846244/151‑1 (MQ=255)
                                                                                                        taagtaaCAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTg                                           >  2:577708/1‑148 (MQ=255)
                                                                                                          agtaaCAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGt                                      >  2:354128/1‑151 (MQ=255)
                                                                                                           gtaaCAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTg                                     <  1:128643/151‑1 (MQ=255)
                                                                                                           gtaaCAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTg                                     >  2:864014/1‑151 (MQ=255)
                                                                                                            taaCAAAAGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTgg                                    >  2:289731/1‑151 (MQ=255)
                                                                                                                aaaaGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCgag                                <  1:439825/151‑1 (MQ=255)
                                                                                                                 aaaGACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCg                                  >  2:297093/1‑148 (MQ=255)
                                                                                                                    gACAATCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCa                            <  2:11605/151‑1 (MQ=255)
                                                                                                                         tCAGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTc                       >  1:115117/1‑151 (MQ=255)
                                                                                                                           aGGGCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTcc                      >  1:371722/1‑150 (MQ=255)
                                                                                                                             ggCGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAgg                   <  2:903164/151‑1 (MQ=255)
                                                                                                                               cGTAAATAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTATTGTTTGAAGGTGGCGAGATCACGTTCCAGGCg                 >  1:783546/1‑151 (MQ=255)
                                                                                                                                     tAGCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGtt           >  2:832594/1‑151 (MQ=255)
                                                                                                                                       gCCCTGATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGTTg          >  2:172311/1‑150 (MQ=255)
                                                                                                                                            gATAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGTTGAAGAAt    <  2:719781/151‑1 (MQ=255)
                                                                                                                                              tAACACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGTTGAAGAATg   >  2:309773/1‑150 (MQ=255)
                                                                                                                                               aaCACGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGTTGAAGAATGa  <  1:271933/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCTCATCATTACCCTAACTGACGGAAATTTCCAAGTCCGCATAGTTTAGCGATAATCCTCGTGAATGCTATAGGATAAGTGATAAAAATAATGAATTGTTGTGTAAGTAACAAAAGACAATCAGGGCGTAAATAGCCCTGATAACAAGATGTTAGTCGCTGCCAAAGAGATCGCGGGTGTATACCTTATCTGCCACATCCTTAAGCTCTTCTGCCATTCGGTTAGAGATAATGACGTCGGCTTGTTGTTTGAAGGTGGCGAGATCACGTTCCAGGCGAGAGTTGAAGAATGA  >  NC_000913/2098295‑2098587

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