Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,287,320 G→A 10.0% intergenic (+18/+23) PP_2889 → / ← xthA transmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,287,3200GA10.0% 146.4 / 5.1 57intergenic (+18/+23)PP_2889/xthAtransmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base A (4/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTCGAT  >  NC_002947/3287176‑3287465
                                                                                                                                                |                                                                                                                                                  
gCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCt                                                                                                                                               >  1:294778/1‑150 (MQ=255)
   gCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGAGTGGGGCGCGGGCGCTGCGGACGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTgg                                                                                                                                            >  1:1377895/1‑150 (MQ=255)
          gATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTc                                                                                                                                    >  2:1131914/1‑150 (MQ=255)
           aTGNCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGCGCGAAAGGTATGGTGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCa                                                                                                                                   >  2:1203989/1‑150 (MQ=255)
           aTGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCa                                                                                                                                   >  2:858503/1‑150 (MQ=255)
            tGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCAc                                                                                                                                  >  1:1372179/1‑150 (MQ=255)
            gggccgtcTGGGGGCGGGGTGTTGGTCGGTTGGGGTTTATAATTTGGCCTGGGTCGGGGCGGCGGACGTGGGGCGTGGGGGGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCAc                                                                                                                                  <  1:1657283/143‑1 (MQ=255)
               ccAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCt                                                                                                                               >  1:1363064/1‑150 (MQ=255)
                 aGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTCTCAATTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTaa                                                                                                                             >  1:912414/1‑150 (MQ=255)
                 aGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTaa                                                                                                                             >  2:1561695/1‑150 (MQ=255)
                  gcgggAGGCGCGGTGCTGGTGGCATGGGTGTTATAACTTTGCGTTGGGCGGCGGGGGGGACGGTGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAAt                                                                                                                            <  2:96064/147‑1 (MQ=255)
                                       gcgggggcCTTATAATTTGGCTTGGGGCGGGGGGGGGGAGGGTGGGGGTGGGGGGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATc                                                                                                       <  2:1072267/146‑1 (MQ=255)
                                         aTGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCgg                                                                                                     >  2:335394/1‑150 (MQ=255)
                                                            cTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCttt                                                                                  >  1:1907903/1‑150 (MQ=255)
                                                                gTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCCCGCTAATGCGGTCTCATCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCa                                                                              >  2:1452427/1‑150 (MQ=255)
                                                                   ggcgcGGGGGAGGGTGGGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATcc                                                                           <  2:1363064/149‑1 (MQ=255)
                                                                            gACGATGTGCGTGGGGCGGTGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCACATCGGCGCATTGTCCGACGGCTTTTCCATCCCGCGCAGCt                                                                  <  1:366762/150‑1 (MQ=255)
                                                                               gATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGCGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCAt                                                               >  2:1833831/1‑150 (MQ=255)
                                                                               gATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCAt                                                               >  2:1273325/1‑150 (MQ=255)
                                                                                      ggggggggggggggggggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCAc                                                        <  2:1372179/148‑1 (MQ=255)
                                                                                      gTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCAc                                                        >  1:1194401/1‑150 (MQ=255)
                                                                                           gggggggTGGGGGAGGCGTTGAGGTTCGGGTTGTAAGTGCGGCCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCg                                                   <  1:1561695/148‑1 (MQ=255)
                                                                                                    ggggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGa                                          >  2:532891/1‑150 (MQ=255)
                                                                                                     gggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGAt                                         <  2:294778/150‑1 (MQ=255)
                                                                                                     gggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGAt                                         >  2:602435/1‑150 (MQ=255)
                                                                                                     gggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGAt                                         >  2:1901510/1‑150 (MQ=255)
                                                                                                     gggAGGCGATGAAGTAGGGGTTGTAAGTGTGGGCCAACGGTGTGCACCTTGGGCCCGACAACCTAATAAGCTCCACTCAAGCCAGAACGGAGCATGGTCCGACGGCTTTTCCATCCAGCACAGCTCCTAGTCCAAGCCCGCTGCCTTGAt                                         >  2:857843/1‑150 (MQ=255)
                                                                                                      ggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATg                                        >  2:1147587/1‑150 (MQ=255)
                                                                                                          gCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTccagcc                                                                                                            >  2:1597119/1‑78 (MQ=255)
                                                                                                          gCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTccagcc                                                                                                            <  1:1597119/78‑1 (MQ=255)
                                                                                                             tggAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTATGACAAAACCTTCGATTACACTTGTCAGCTCATCTCCAGCCAATGCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCa                                 <  1:1869067/148‑1 (MQ=255)
                                                                                                              tGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCAc                                <  2:669144/150‑1 (MQ=255)
                                                                                                               gAGGTGGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCAcc                               <  2:1145785/150‑1 (MQ=255)
                                                                                                                 cgatctggTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAg                             <  1:584619/144‑1 (MQ=255)
                                                                                                                  gTTGGGGTTGTAAGTGTGGGCCCGAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGAGCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGcgc                                                    >  2:949397/1‑124 (MQ=255)
                                                                                                                  gTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGccc                                                    <  1:949397/126‑1 (MQ=255)
                                                                                                                        gTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCAGCCTTGCTGCGCGGCACCGgcagatcg                      >  2:584619/1‑144 (MQ=255)
                                                                                                                          tGTAAGTGTGGGCCCACGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa                    <  1:858503/150‑1 (MQ=255)
                                                                                                                          tGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa                    <  2:1377895/150‑1 (MQ=255)
                                                                                                                          tGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa                    <  1:335394/150‑1 (MQ=255)
                                                                                                                           gTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAg                   <  2:1517203/150‑1 (MQ=255)
                                                                                                                             aaGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGc                 >  1:902858/1‑150 (MQ=255)
                                                                                                                                tgtgGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCAt              >  2:1890059/1‑150 (MQ=255)
                                                                                                                                    ggCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATc          >  2:1129710/1‑150 (MQ=255)
                                                                                                                                    ggCCAAAGGTGGGCACCTTGCGCTCCACAACCTAACCAGCTCAACTCCAGCCAGTTCGGCGCATCGTCCCACTGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATc          <  2:729335/150‑1 (MQ=255)
                                                                                                                                     gCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCATCTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGAGCa         <  1:1616119/150‑1 (MQ=255)
                                                                                                                                      cccAAGGTGTGCACCTTGGGCCCGTCACCCTACAGCGCTCGTCTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAg        <  2:912414/150‑1 (MQ=255)
                                                                                                                                      cccAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa                    <  2:1847751/138‑1 (MQ=255)
                                                                                                                                      cccAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGa                    >  1:1847751/1‑138 (MQ=255)
                                                                                                                                      cccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGgc                                                                                                    <  1:422874/58‑1 (MQ=37)
                                                                                                                                      cccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGgc                                                                                                    >  2:422874/1‑58 (MQ=38)
                                                                                                                                      cccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCTCCAGCTGCTGTGAGGCCATGATCAg        >  1:1704112/1‑150 (MQ=255)
                                                                                                                                      cccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAg        >  2:1400302/1‑150 (MQ=255)
                                                                                                                                       ccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAg                                                             <  1:1566073/96‑1 (MQ=255)
                                                                                                                                       ccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAg                                                             >  2:1566073/1‑96 (MQ=255)
                                                                                                                                          aGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTCg    >  2:794020/1‑150 (MQ=255)
                                                                                                                                            gtggGCACCTTTGGCCCGTCACCCTAGTCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTACACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTCGAt  <  1:761980/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCAGCGGGTCGATGGCCAGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTCGAT  >  NC_002947/3287176‑3287465

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