Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,078,189 A→G 10.3% G223G (GGT→GGC cmtA ← putative mannitol‑specific PTS IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,078,1890AG10.3% 136.4 / 6.6 58G223G (GGT→GGCcmtAputative mannitol‑specific PTS IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (3/3);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

ATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGC  >  NC_000913/3078046‑3078333
                                                                                                                                               |                                                                                                                                                
atGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                             <  1:368167/149‑1 (MQ=255)
atGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTCCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTgg                                                                                                                                             <  2:323213/149‑1 (MQ=255)
  ggTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTTGTGGCTCCCACCGAGGAAGTGAATAATCCTACCCCCGTGCGCAGCAGGTCGTCTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTgggg                                                                                                                                           <  1:189958/149‑2 (MQ=255)
   gTCAGCGGCTTCATCAGCATATACGGGAAACACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggc                                                                                                                                          <  2:663084/149‑3 (MQ=255)
    tCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGTCCGGGCCTggggct                                                                                                                                         <  1:345007/149‑4 (MQ=255)
    tCAGCGGCTCCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACTGAGGAAGCAAATAATCATTCCCCCGGGCACACTACGCGTAGCGCTCCCTTTACCAACGAAGGTAGACGCCAGCCGCAGCCCCAGGCCCGGCCCGGGGCTa                                                                                                                                         >  2:326072/1‑149 (MQ=255)
      aGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTaga                                                                                                                                       <  2:788957/149‑1 (MQ=255)
            ttCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccac                                                                                                                                 <  1:742586/149‑1 (MQ=255)
             tcatcaGCACATACGGAAAATGCAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                                                                <  1:491408/149‑1 (MQ=255)
             tcagccgCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGccacc                                                                                                                                <  1:319867/143‑1 (MQ=255)
              catcaGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCACCAGGCCCGGACCTGGGTTAGAGGccacca                                                                                                                               <  1:418572/149‑1 (MQ=255)
                 caGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATaa                                                                                                                            <  1:739801/149‑1 (MQ=255)
                   gCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAaag                                                                                                                          <  1:691570/149‑1 (MQ=255)
                      caTACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGAAAACGCCAGCAGCAGCCCCAGGCCCGGGCCTGGGTTAGAGGCCCCCATAaagaag                                                                                                                       >  1:344135/1‑149 (MQ=255)
                            ggAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCAACGCCCCGGGACCAGAACGGTTACTCATCCCCTCACCAAACAAGGGAAACGCCAGCAGCAACCCCAGGCCCGCACCTGGGTTAAAGGCCACCAAAAACAACATGGAt                                                                                                                 >  2:349862/1‑149 (MQ=255)
                             gAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGAtt                                                                                                                >  1:128493/1‑149 (MQ=255)
                               aaTACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTa                                                                                                              <  1:534036/149‑1 (MQ=255)
                                aTACAGTTCGGGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGTCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTAc                                                                                                             <  1:389154/149‑1 (MQ=255)
                                    aGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGtt                                                                                                         <  1:316108/149‑1 (MQ=255)
                                     gTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTa                                                                                                        >  2:677194/1‑149 (MQ=255)
                                     gTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAAAGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCACCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTa                                                                                                        >  2:177445/1‑149 (MQ=255)
                                      ttCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTaa                                                                                                       <  1:230797/149‑1 (MQ=255)
                                         gTGGATCCCACCGAGGAAGTGAATAATCATAGCCCCGAGCGCATCCAGTTTACTCACACGTGTACCAAACAAGGTAAACGCAAGCATCACTCTCACGCTCGATCTTGGCTTAGAGGCCATCTTTAAGAGGTTGCGTTTACCGTTCAcca                                                                                                    >  2:790387/1‑148 (MQ=255)
                                          tGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGa                                                                                                   <  1:70598/149‑1 (MQ=255)
                                           ggATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGaa                                                                                                  <  1:7680/149‑1 (MQ=255)
                                            gATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGACAGCAGCAGCCCCAGGCCCGGACGTGGGTGAGAGGCCACCATACAGAAGACGGATTTACCGATAACCGAAg                                                                                                 >  2:645243/1‑149 (MQ=255)
                                                  aCCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCCGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGtt                                                                                           <  1:614853/149‑1 (MQ=255)
                                                   ccGAGGAAGTGCATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGGGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTg                                                                                          <  2:610978/149‑1 (MQ=255)
                                                          aGTGACCAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTccc                                                                                   >  2:531675/1‑149 (MQ=255)
                                                            tGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTGAGACTCATGGCACACTGCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGGCTGTTGCATTCCCAg                                                                                 >  1:790229/1‑149 (MQ=255)
                                                            tGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCCGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAg                                                                                 >  1:464863/1‑149 (MQ=255)
                                                             gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGTGGGATTTCCCTTTAACCGAAGCCTGTTGCATTCCAAGc                                                                                >  1:150120/1‑149 (MQ=255)
                                                             gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGc                                                                                >  1:197756/1‑149 (MQ=255)
                                                             gAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGc                                                                                >  2:394028/1‑149 (MQ=255)
                                                                     atcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGataata                                                                        <  1:121177/149‑1 (MQ=255)
                                                                      tcGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCCTTCCCAGCGGATAATAg                                                                       >  2:79828/1‑149 (MQ=255)
                                                                                     gaCGTTTCCTCAGCCCTTTACCAAACAAGGTAAACGCCAGCAGCACCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTcatccctc                                                        <  2:450301/148‑5 (MQ=255)
                                                                                       cGTTTACTCATCCCTTTACCAAACAACGTAAACAAAAATTGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGc                                                      >  1:792953/1‑149 (MQ=255)
                                                                                               cctCCCTTTACCAAACAAGGTCAACGCCAGCAGCAGCCCCAGGCCCCGCCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                              <  2:153157/147‑1 (MQ=255)
                                                                                               cATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGa                                              <  1:73120/149‑1 (MQ=255)
                                                                                                aTCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAg                                             >  1:464359/1‑149 (MQ=255)
                                                                                                    cTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGGTCGATCGCCTTATTGAGGaaa                                         >  2:683015/1‑149 (MQ=255)
                                                                                                        aCCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATGCAGGAAAAGCa                                     >  2:538607/1‑149 (MQ=255)
                                                                                                          cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGGTAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCAcg                                   >  1:713652/1‑148 (MQ=255)
                                                                                                          cAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACt                                   >  1:799745/1‑149 (MQ=255)
                                                                                                           aaaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACtt                                  >  2:735821/1‑149 (MQ=255)
                                                                                                            aaCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACttt                                 >  2:42141/1‑149 (MQ=255)
                                                                                                             aCAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAGGCACTTTc                                <  1:88853/149‑1 (MQ=255)
                                                                                                              cAAGGTAAACGCTATCAGCATCCCCAGGCCCGGACGTGGGTTAGATGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCg                               <  1:219634/149‑1 (MQ=255)
                                                                                                                          cagcagcagCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATtga                   >  1:72451/1‑149 (MQ=255)
                                                                                                                              ctctaacccAGGTCCGGGCCTGGGGTAGAGGCCAGCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATGGCATGCTAGAGGAAAAGCACTTTCGCCGGCTCATGTAtgac               >  1:403272/7‑149 (MQ=255)
                                                                                                                                 agcCCCAGGACCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGAtgactga            <  1:499815/149‑1 (MQ=255)
                                                                                                                                   acccAGGTCCGGGCCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGaca          >  1:739735/2‑149 (MQ=255)
                                                                                                                                     ccAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacacc        >  1:273569/1‑147 (MQ=255)
                                                                                                                                      cAGGACCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGacaaca       <  2:581729/149‑1 (MQ=255)
                                                                                                                                          cccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATgg   <  1:685745/149‑1 (MQ=255)
                                                                                                                                          cccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTATCCCAGCGTTATTGGTGACGGacaaccaggg   >  1:742957/1‑144 (MQ=255)
                                                                                                                                           ccGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGc  >  1:656018/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ATGGTCAGCGGCTTCATCAGCACATACGGGAAATACAGTTCGTGGATCCCACCGAGGAAGTGAATAATCATCGCCCCGGGCGCAGAACGTTTACTCATCCCTTTACCAAACAAGGTAAACGCCAGCAGCAGCCCCAGGCCCGGACCTGGGTTAGAGGCCACCATAAAGAAGATGGATTTACCGTTAACCGAAGCCTGTTGCATTCCCAGCGGATAATAGACGCCCTGATCGATCGCATTATTGAGGAAAAGCACTTTCGCCGGTTCATTGATGACTGACAACAATGGC  >  NC_000913/3078046‑3078333

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