Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 170,165 G→A E130K (GAA→AAA)  fhuC → iron(III) hydroxamate ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913170,1650GA96.5% 191.5 / ‑3.4 58E130K (GAA→AAA) fhuCiron(III) hydroxamate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (33/23);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAC  >  NC_000913/170035‑170301
                                                                                                                                  |                                                                                                                                        
aGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTAtct                                                                                                                                 <  2:1222543/140‑1 (MQ=255)
        aTTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGtt                                                                                                                         >  1:285413/1‑140 (MQ=255)
              cGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTa                                                                                                                   >  2:79139/1‑140 (MQ=255)
                  gcagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAAc                                                                                                               <  1:79139/140‑1 (MQ=255)
                   cagcTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAAcc                                                                                                              >  2:632630/1‑140 (MQ=255)
                      cTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCt                                                                                                           <  2:597333/140‑1 (MQ=255)
                       ttCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTg                                                                                                          <  1:561487/140‑1 (MQ=255)
                        tcctccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTgg                                                                                                         >  1:225401/1‑140 (MQ=255)
                             cGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCAt                                                                                                    >  2:98206/1‑140 (MQ=255)
                              ggCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATc                                                                                                   <  2:225401/140‑1 (MQ=255)
                                cAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCgg                                                                                                 <  1:842665/140‑1 (MQ=255)
                                  gAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCt                                                                                               <  2:715337/140‑1 (MQ=255)
                                   aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                              >  1:627864/1‑140 (MQ=255)
                                   aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                                                                                              >  1:426115/1‑140 (MQ=255)
                                          gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                       >  1:682209/1‑140 (MQ=255)
                                           aCCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAg                                                                                      >  1:261785/1‑140 (MQ=255)
                                            ccGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGt                                                                                     <  1:1041412/140‑1 (MQ=255)
                                             cgtgcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtc                                                                                    >  1:283107/1‑140 (MQ=255)
                                               tgcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctc                                                                                  >  2:301671/1‑140 (MQ=255)
                                                 cgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTAGTCGATAGtctctc                                                                                >  2:813335/1‑140 (MQ=255)
                                                   tgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTg                                                                              >  1:135079/1‑140 (MQ=255)
                                                          gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                       >  1:653691/1‑124 (MQ=255)
                                                          gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                                                                                       <  2:653691/124‑1 (MQ=255)
                                                          gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGa                                                                       >  1:625372/1‑140 (MQ=255)
                                                          gtggAGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGa                                                                       >  1:485902/1‑140 (MQ=255)
                                                            ggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATGGTCTCTCTGGCGGCGAAc                                                                     <  2:1026477/140‑1 (MQ=255)
                                                             gCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACg                                                                    >  1:531516/1‑140 (MQ=255)
                                                               gATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTc                                                                  >  2:913171/1‑140 (MQ=255)
                                                                aTTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCa                                                                 >  2:82049/1‑140 (MQ=255)
                                                                 ttGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAg                                                                <  1:974803/140‑1 (MQ=255)
                                                                 ttGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAg                                                                <  2:426115/140‑1 (MQ=255)
                                                                 ttGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAg                                                                >  2:1036306/1‑140 (MQ=255)
                                                                 ttGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTAGTCGATAGTCTCTCTGGCGGCGAACGTCAg                                                                <  1:813335/140‑1 (MQ=255)
                                                                         tACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTg                                                        <  1:409937/140‑1 (MQ=255)
                                                                           cccGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGa                                                      <  2:942586/140‑1 (MQ=255)
                                                                               tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGc                                                  <  1:632630/140‑1 (MQ=255)
                                                                               tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGc                                                  >  1:853964/1‑140 (MQ=255)
                                                                               tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGc                                                  >  2:931848/1‑140 (MQ=255)
                                                                               tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGc                                                  >  2:206460/1‑140 (MQ=255)
                                                                               tGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGc                                                  >  1:736030/1‑140 (MQ=255)
                                                                                      gcgcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtg                                           >  1:1048324/1‑140 (MQ=255)
                                                                                       cgcgcTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtgg                                          <  1:527783/140‑1 (MQ=255)
                                                                                            tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                     <  2:135079/140‑1 (MQ=255)
                                                                                            tGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                     >  1:1091104/1‑140 (MQ=255)
                                                                                                 cgcTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGa                                <  2:682209/140‑1 (MQ=255)
                                                                                                   cTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                       >  1:1063572/1‑131 (MQ=255)
                                                                                                   cTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                       <  2:1063572/131‑1 (MQ=255)
                                                                                                   cTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                       >  1:425299/1‑131 (MQ=255)
                                                                                                   cTTTGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCtggtgg                                       <  2:425299/131‑1 (MQ=255)
                                                                                                      tGGGGCGGCAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGcc                           >  2:139837/1‑140 (MQ=255)
                                                                                                             gcAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                     <  2:660657/123‑1 (MQ=255)
                                                                                                             gcAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGgcg                                     >  1:660657/1‑123 (MQ=255)
                                                                                                             gcAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCt                    >  2:701634/1‑140 (MQ=255)
                                                                                                              cAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGa                                <  1:655171/127‑1 (MQ=255)
                                                                                                              cAGATCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGa                                >  2:655171/1‑127 (MQ=255)
                                                                                                                 aTCGCGAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTg                <  2:67065/140‑1 (MQ=255)
                                                                                                                   cgcgAAAAAGTCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCt              >  1:860120/1‑140 (MQ=255)
                                                                                                                             tCGAGAAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCg    >  1:564603/1‑140 (MQ=255)
                                                                                                                               gagaAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAAGATAGCCGTTGTCTGTTGCTCGACGAACCGAc  <  2:283107/140‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
AGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGAACGTCAGCGGGCGTGGATCGCCATGCTGGTGGCGCAGGATAGCCGTTGTCTGTTGCTCGACGAACCGAC  >  NC_000913/170035‑170301

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