Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 619,414 T→G 7.3% G262G (GGA→GGC) ‡ fepC ← ferrienterobactin ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913619,4140TG 7.3% 173.5 / 7.0 55G262G (GGA→GGC) ‡fepCferrienterobactin ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base G (2/2);  total (26/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.57e-01

GGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAT  >  NC_000913/619267‑619562
                                                                                                                                                   |                                                                                                                                                    
ggCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTc                                                                                                                                                     <  1:417903/149‑1 (MQ=255)
 gCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTcc                                                                                                                                                    >  1:527411/1‑149 (MQ=255)
           tgATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGcacc                                                                                                                                          <  4:267049/149‑1 (MQ=255)
                 tttCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAgc                                                                                                                                    <  4:331448/149‑1 (MQ=255)
                      gcgTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGt                                                                                                                               >  2:511125/1‑149 (MQ=255)
                          tGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCg                                                                                                                           <  4:193684/149‑1 (MQ=255)
                          tGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCg                                                                                                                           <  3:533209/149‑1 (MQ=255)
                               ggcggGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAAc                                                                                                                      <  2:52134/149‑1 (MQ=255)
                                    gATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGAt                                                                                                                 <  3:97531/149‑1 (MQ=255)
                                          ggCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGAAGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCaa                                                                                                           >  1:486908/1‑149 (MQ=255)
                                               gtggtgGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATc                                                                                                      <  3:100777/149‑1 (MQ=255)
                                                       gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgcg                                                                                              >  1:363073/1‑149 (MQ=255)
                                                       gTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGCTCGGACGAACGGc                                                                                                                                             <  1:365003/102‑1 (MQ=255)
                                                        tGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgcgc                                                                                             >  2:66125/1‑149 (MQ=255)
                                                           tGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGa                                                                                          <  2:47853/149‑1 (MQ=255)
                                                            gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGAc                                                                                         <  1:537963/149‑1 (MQ=255)
                                                            gTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGAc                                                                                         <  1:395903/149‑1 (MQ=255)
                                                                   tATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAAt                                                                                  >  4:158425/1‑149 (MQ=255)
                                                                       cGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATgcgc                                                                              <  4:344622/149‑1 (MQ=255)
                                                                       cGATGGAATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCAGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATgcgc                                                                              <  4:354027/149‑1 (MQ=255)
                                                                              aTCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCa                                                                       >  2:559594/1‑149 (MQ=255)
                                                                               tCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCATCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgc                                                                                               <  3:562340/124‑1 (MQ=255)
                                                                               tCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAgc                                                                                               >  4:562340/1‑124 (MQ=255)
                                                                                  aGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTc                                                                   >  4:202298/1‑149 (MQ=255)
                                                                                      aaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATaa                                                                                    <  2:210187/128‑1 (MQ=255)
                                                                                      aaCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATaa                                                                                    >  1:210187/1‑128 (MQ=255)
                                                                                       aCAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTcagcag                                                              >  1:233228/1‑149 (MQ=255)
                                                                                        cAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGt                                                             <  4:146248/149‑1 (MQ=255)
                                                                                              gCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAAt                                                       <  1:559594/149‑1 (MQ=255)
                                                                                                      ggCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCTGCACCACAAGCGGCGGTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCTCTCAATCAGTTCAGCAGTGACAATCTCCTTCg                                               <  4:367945/149‑1 (MQ=255)
                                                                                                      ggCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCg                                               <  3:388102/149‑1 (MQ=255)
                                                                                                       gCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCgg                                              >  1:173745/1‑149 (MQ=255)
                                                                                                         aGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGgcg                                            >  1:14593/1‑149 (MQ=255)
                                                                                                             cACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGa                                                                                          <  4:43286/99‑1 (MQ=255)
                                                                                                             cACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGa                                                                                          >  3:43286/1‑99 (MQ=255)
                                                                                                              aCTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCTATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTccc                                       <  2:233228/149‑1 (MQ=255)
                                                                                                                 tAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCAGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                    <  1:227835/149‑1 (MQ=255)
                                                                                                                    tttAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                    >  2:185435/1‑146 (MQ=255)
                                                                                                                    tttAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGa                                    <  1:185435/146‑1 (MQ=255)
                                                                                                                     ttAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCAAACGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCACAGACCATAAAGGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCGCTGAGcaa                                >  4:360119/1‑149 (MQ=255)
                                                                                                                      tAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGcaac                               >  1:549275/1‑149 (MQ=255)
                                                                                                                          tATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtt                           >  1:54590/1‑149 (MQ=255)
                                                                                                                            tttGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAAtttt                         <  2:363073/149‑1 (MQ=255)
                                                                                                                              tGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcc                       >  3:221240/1‑149 (MQ=255)
                                                                                                                                  gTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTttcccttc                   >  1:7794/1‑149 (MQ=255)
                                                                                                                                      gtggTGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgc               >  1:246076/3‑149 (MQ=255)
                                                                                                                                        ggTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaa             >  3:474037/1‑149 (MQ=255)
                                                                                                                                        ggTGCCGCTCGGACGAACGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaa             >  4:453723/1‑149 (MQ=255)
                                                                                                                                         gTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATgcgc                                                                              <  2:530273/83‑1 (MQ=255)
                                                                                                                                         gTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATgcgc                                                                              >  1:530273/1‑83 (MQ=255)
                                                                                                                                         gTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaat            <  1:183297/149‑1 (MQ=255)
                                                                                                                                         gTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaat            <  1:511125/149‑1 (MQ=255)
                                                                                                                                             cGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgcaa        <  2:54590/149‑1 (MQ=255)
                                                                                                                                              gTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCgcaatgcaat       >  2:328024/1‑149 (MQ=255)
                                                                                                                                               ttCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATc      >  1:8094/1‑149 (MQ=255)
                                                                                                                                                   tCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGTCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAt  >  4:512870/1‑149 (MQ=255)
                                                                                                                                                   tCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAt  >  1:346226/1‑149 (MQ=255)
                                                                                                                                                   tCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAt  >  4:548017/1‑149 (MQ=255)
                                                                                                                                                   tCCGAGCGGCAACACAAGCGGCGTTCCGGCAACTGGAACGTCAAAGATCATGCAGCGCAGACAATAAATGCGCGCAATCAATTCAGCAGTGACAATCTCCTTCGGCGCCCCCAGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAt  <  3:233191/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGCGATTATTGTGATCCTTTCCGCGTTGATCGGCGGGATGGTGGCTTGTGGTGGCGTGCTGTTGCGCTATGCGATGGCATCCAGAAAACAGGATGCCATGATGGCAGACCACTTAGTTTAACTATTTGCGGTTGAAGGTGCCGTTCGTCCGAGCGGCACCACAAGCGGCGTTCCGGCAACTGGATCGTCAATGATCATGCAGCGCAGACCATAAATGCGCTCAATCAGTTCAGCAGTGACAATCTCCTTCGGCGCTCCCTGAGCAACAATTTTCCCTTCCCGCAATGCAATCAAAT  >  NC_000913/619267‑619562

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