Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,886,454 C→T 7.0% G468D (GGC→GAC)  PP_3430 ← two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,886,4540CT 7.0% 165.1 / 2.6 57G468D (GGC→GAC) PP_3430two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGCCCCCAGCGCCCGCAGCAACAGGGG  >  NC_002947/3886312‑3886602
                                                                                                                                              |                                                                                                                                                    
gAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGtt                                                                                                                                                <  7:201565/149‑1 (MQ=255)
gAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGtt                                                                                                                                                >  1:318665/1‑149 (MQ=255)
gAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGtt                                                                                                                                                <  5:44266/149‑1 (MQ=255)
    gCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAg                                                                                                                                            <  8:29319/149‑1 (MQ=255)
        tGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCAc                                                                                                                                        >  3:259337/1‑149 (MQ=255)
          gCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACg                                                                                                                                       >  3:124599/1‑148 (MQ=255)
          gCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACg                                                                                                                                       <  4:124599/148‑1 (MQ=255)
            gATATCCGCCGCGGCCATGCCATGGACATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATCGTGACATTGCCCTGTTCCAGGccagctg                                                                                                                                    >  8:142448/1‑143 (MQ=255)
                tCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCAc                                                                                                                                 <  1:131794/148‑1 (MQ=255)
                 ggatgatggCCATGGCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCAcgc                                                                                                                               >  5:247559/8‑149 (MQ=255)
                      gcgGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCATGCGCAtt                                                                                                                          <  4:259337/149‑1 (MQ=255)
                      gcgGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCAtt                                                                                                                          <  5:384818/149‑1 (MQ=255)
                       cgGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCAGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTg                                                                                                                         <  1:72570/149‑1 (MQ=255)
                            atgccatgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGcaca                                                                                                                    >  7:287992/1‑149 (MQ=255)
                               ccatgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACgcgc                                                                                                                             >  3:239653/1‑137 (MQ=255)
                               ccatgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGcacac                                                                                                                   >  7:153203/1‑147 (MQ=255)
                               ccatgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACaa                                                                                                                 >  7:153200/1‑149 (MQ=255)
                               ccatgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGGCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACgcgc                                                                                                                             <  4:239653/137‑1 (MQ=255)
                                 atgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGt                                                                                                               <  2:322312/149‑1 (MQ=255)
                                  tgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGca                                                                                                                      <  4:265254/141‑1 (MQ=255)
                                  tgGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGATGCTGATGGTGACATTGCCCTCTTCCAGTCACGCTTCACGCGCATTGGcacatatcg                                                                                                               >  3:265254/1‑143 (MQ=255)
                                       atcatcCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGTCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATc                                                                                                          >  7:286313/1‑148 (MQ=255)
                                          atcCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAg                                                                                                       <  8:225684/148‑1 (MQ=255)
                                                 aCATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACTCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATtctc                                                                                               >  2:228433/1‑149 (MQ=255)
                                                       aaGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTg                                                                                         <  2:318665/149‑1 (MQ=255)
                                                        aGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTgg                                                                                        <  1:144291/149‑1 (MQ=255)
                                                          gCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTgggg                                                                                      <  5:171648/149‑1 (MQ=255)
                                                                 aGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTg                                                                                         >  8:174197/1‑139 (MQ=255)
                                                                 aGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTg                                                                                         <  7:174197/139‑1 (MQ=255)
                                                                  gTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTc                                                                                >  2:36942/1‑147 (MQ=255)
                                                                   tCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGTCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAg                                                                             <  8:286313/149‑1 (MQ=255)
                                                                   tCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAg                                                                             >  6:310021/1‑149 (MQ=255)
                                                                      cccGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCt                                                                          <  8:153203/149‑1 (MQ=255)
                                                                      cccGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCt                                                                          <  8:153200/149‑1 (MQ=255)
                                                                                   cACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGt                                                             <  1:36942/149‑1 (MQ=255)
                                                                                    aCCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGTCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAg                                                               <  5:28394/146‑1 (MQ=255)
                                                                                    aCCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGTCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAg                                                               >  6:28394/1‑146 (MQ=255)
                                                                                      cTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTc                                                                                                                                               >  3:1597/1‑64 (MQ=255)
                                                                                      cTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTc                                                                                                                                               <  4:1597/64‑1 (MQ=255)
                                                                                             cGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTc                                                     >  1:282667/1‑147 (MQ=255)
                                                                                              gggAAAAACGCGGTCAGCCGTTTGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTAACGCGCATTGGCACACAAGTGGATCAAGGCAATCTCTAGCAGGGGTACGTCCAGGCTCACTGCCCAAGGt                                                             >  7:236913/1‑138 (MQ=255)
                                                                                                gAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAg                                                               >  2:235251/1‑134 (MQ=255)
                                                                                                gAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAg                                                               <  1:235251/134‑1 (MQ=255)
                                                                                                    aaCGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCa                                            <  6:247559/149‑1 (MQ=255)
                                                                                                    aaCGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCa                                            <  8:15735/149‑1 (MQ=255)
                                                                                                      cgcgGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGt                                          <  8:287992/149‑1 (MQ=255)
                                                                                                      cgcgGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGt                                          >  7:195634/1‑149 (MQ=255)
                                                                                                          gTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCaa                                      <  8:195634/149‑1 (MQ=255)
                                                                                                           tCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGAACAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCa                                                   <  4:6328/135‑1 (MQ=255)
                                                                                                           tCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGAACAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCa                                                   >  3:6328/1‑135 (MQ=255)
                                                                                                            cAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAAc                                    >  4:131469/1‑149 (MQ=255)
                                                                                                                      ttGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGc                          >  2:77264/1‑149 (MQ=255)
                                                                                                                        gACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGccc                        >  1:354386/1‑149 (MQ=255)
                                                                                                                        gACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGcc                         <  2:282667/148‑1 (MQ=255)
                                                                                                                                  aTGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCAtct                              <  1:227868/133‑1 (MQ=255)
                                                                                                                                  aTGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCAtct                              >  2:227868/1‑133 (MQ=255)
                                                                                                                                             gCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGCGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGCCCCAAGCGCCCGCAGCAACAgg    <  7:142448/148‑1 (MQ=255)
                                                                                                                                              cccTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGCCCCCAGCGCCCGCAGCAACAggg   <  2:354386/148‑1 (MQ=255)
                                                                                                                                              cccTGTTCCAGGCACGCTTCACGAGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTTCCCCCAGCGCCCCCAGCAACAggcg  >  8:143691/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GAAAGCCCTGGCGATATCCGCCGCGGCCATGCCATGGCCATCATCCTCGACATGCAAGGCCACATAGTCGCCCGGCGGCAAGCCACCTTCATCCGGGAAAAACGCGGTCAGCCGTTCGTTGACGCTGCTGATGGTGACATTGCCCTGTTCCAGGCACGCTTCACGCGCATTGGCACACAAGTGGATCAAGGCATTCTCTAGCTGGGGTACGTCCAGGCTCACTGCCCAAGGTGCAACGTCCAGTTGCCAGTGCAAACGCATCTCTGCCCCCAGCGCCCGCAGCAACAGGGG  >  NC_002947/3886312‑3886602

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