Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,217,099 T→C 7.2% intergenic (+103/‑229) ymgC → / → ycgG Blue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,217,0990TC 7.2% 223.9 / 6.5 69intergenic (+103/‑229)ymgC/ycgGBlue light, low temperature and stress induced protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (27/37);  new base C (2/3);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

ACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGA  >  NC_000913/1216951‑1217209
                                                                                                                                                    |                                                                                                              
acTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgt                                                                                                                <  3:236991/149‑1 (MQ=255)
        atatCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGaa                                                                                                        <  2:126642/149‑1 (MQ=255)
         tatCATTATTAATGACAGGGAATCCGTGCTCTCTTAGCACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAAt                                                                                                       <  3:1230653/149‑1 (MQ=255)
         tatCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAAt                                                                                                       <  1:820137/149‑1 (MQ=255)
         tatCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAAt                                                                                                       <  1:762389/149‑1 (MQ=255)
            cATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    <  2:1801813/149‑1 (MQ=255)
              ttattaATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  >  3:1815191/1‑149 (MQ=255)
              ttattaATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  <  2:78756/149‑1 (MQ=255)
              ttattaATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  <  4:487452/149‑1 (MQ=255)
              ttattaATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  <  4:1815191/149‑1 (MQ=255)
                attaATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagagc                                                                                                >  1:604854/1‑148 (MQ=255)
                  taATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtc                                                                                              <  3:1664194/149‑1 (MQ=255)
                   aaTGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtct                                                                                             <  4:1196181/149‑1 (MQ=255)
                    aTGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtctc                                                                                            >  2:493719/1‑149 (MQ=255)
                     tGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagagcgagaa                                                                                           <  2:604854/149‑7 (MQ=255)
                         aGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGctct                                                                                       <  4:492743/149‑1 (MQ=255)
                          gggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGaa                                                                                                        <  4:1243251/131‑1 (MQ=255)
                          gggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGaa                                                                                                        >  3:1243251/1‑131 (MQ=255)
                           ggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATg                                                                                                      >  2:1549129/1‑132 (MQ=255)
                           ggAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATg                                                                                                      <  1:1549129/132‑1 (MQ=255)
                              aTCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCAtt                                                                                  <  3:245582/149‑1 (MQ=255)
                                 cGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                               <  1:493719/149‑1 (MQ=255)
                                   tGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  <  4:1490351/128‑1 (MQ=255)
                                   tGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  >  3:1490351/1‑128 (MQ=255)
                                      tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGTGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCa                                                                                    >  2:47513/1‑139 (MQ=255)
                                      tctctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGTGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCa                                                                                    <  1:47513/139‑1 (MQ=255)
                                        tctTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCt                                                                        >  2:1094819/1‑149 (MQ=255)
                                            aGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGagagcgagaa                                                                                           >  2:180959/1‑120 (MQ=255)
                                            aGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGagagcgagaa                                                                                           >  2:1204913/1‑120 (MQ=255)
                                            aGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGagagcgagaa                                                                                           <  1:180959/126‑7 (MQ=255)
                                            aGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGCGTCAAGAATGGagagcgagaa                                                                                           <  1:1204913/126‑7 (MQ=255)
                                                      tAAGGTTTTAGAATAATCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagag                                                                                                 <  1:540358/110‑1 (MQ=255)
                                                      tAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATg                                                                                                      >  2:442244/1‑105 (MQ=255)
                                                      tAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATg                                                                                                      <  1:442244/105‑1 (MQ=255)
                                                      tAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagag                                                                                                 >  2:540358/1‑110 (MQ=255)
                                                             ttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    <  3:1837829/100‑1 (MQ=255)
                                                             ttAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    >  4:1837829/1‑100 (MQ=255)
                                                                gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    <  3:1670471/97‑1 (MQ=255)
                                                                gAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    >  4:1670471/1‑97 (MQ=255)
                                                                   tAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    <  4:38220/94‑1 (MQ=255)
                                                                   tAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    >  3:38220/1‑94 (MQ=255)
                                                                     aaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATg                                                                                                      <  4:867687/90‑1 (MQ=255)
                                                                     aaCATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATg                                                                                                      >  3:867687/1‑90 (MQ=255)
                                                                        ttCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  >  2:1406868/2‑91 (MQ=255)
                                                                        ttCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga                                                                                                  <  1:1406868/90‑1 (MQ=255)
                                                                        aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                               >  1:1350717/1‑110 (MQ=255)
                                                                        aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                               <  2:999653/110‑1 (MQ=255)
                                                                        aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                               <  2:1350717/110‑1 (MQ=255)
                                                                        aTCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCtt                                                                               >  1:999653/1‑110 (MQ=255)
                                                                                tATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    <  3:772323/81‑1 (MQ=255)
                                                                                tATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    >  4:772323/1‑81 (MQ=255)
                                                                                  tCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagag                                                                                                 <  4:390297/82‑1 (MQ=255)
                                                                                  tCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagag                                                                                                 >  3:390297/1‑82 (MQ=255)
                                                                                     ccATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtctc                                                                                            <  3:170426/84‑1 (MQ=255)
                                                                                     ccATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTtctc                                                                                            >  4:170426/1‑84 (MQ=255)
                                                                                            gATGCAAATTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaattaa                     >  3:1743397/1‑148 (MQ=255)
                                                                                            gATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAAtaattaa                     <  4:1743397/148‑1 (MQ=255)
                                                                                               gCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTc                              >  3:1550056/1‑136 (MQ=255)
                                                                                               gCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTc                              <  4:1550056/136‑1 (MQ=255)
                                                                                                     tGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGAt           >  4:53894/1‑149 (MQ=255)
                                                                                                          aTTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGa      <  1:1094819/149‑1 (MQ=255)
                                                                                                           tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat                                                                     <  4:1151957/85‑1 (MQ=255)
                                                                                                           tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat                                                                     >  3:1151957/1‑85 (MQ=255)
                                                                                                           tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat                                                                     >  2:1550992/1‑85 (MQ=255)
                                                                                                           tttCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGat                                                                     <  1:1550992/85‑1 (MQ=255)
                                                                                                              cAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGga  <  3:53894/149‑1 (MQ=255)
                                                                                                                      aTGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    >  2:1092933/1‑43 (MQ=255)
                                                                                                                      aTGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga                                                                                                    <  1:1092933/43‑1 (MQ=255)
                                                                                                                                             tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                >  2:298750/1‑40 (MQ=255)
                                                                                                                                             tatCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCt                                                                                <  1:298750/40‑1 (MQ=255)
                                                                                                                                                    |                                                                                                              
ACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTTCTCGCTCTCCATTCTTGACGCCTGATATCCCGCCTAACTTATTTTGTACTGCCTACATTATGCTCAATAATTAATGGGTTAGATCAAGATGGA  >  NC_000913/1216951‑1217209

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