Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 36,174 A→C 7.6% intergenic (+59/+24) PP_0032 → / ← PP_0033 hypothetical protein/undecaprenyl‑glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294736,1740AC 7.6% 129.6 / 4.8 53intergenic (+59/+24)PP_0032/PP_0033hypothetical protein/undecaprenyl‑glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (22/27);  new base C (2/2);  total (24/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 = 1.96e-01

GAGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGATGT  >  NC_002947/36025‑36320
                                                                                                                                                     |                                                                                                                                                  
gAGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTaaaa                                                                                                                                                    >  1:1010990/1‑150 (MQ=255)
gAGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTaaaa                                                                                                                                                    >  1:1146393/1‑150 (MQ=255)
    gATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTg                                                                                                                                                <  2:730052/150‑1 (MQ=255)
           ctTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCttttt                                                                                                                                         <  1:1216974/150‑1 (MQ=255)
           ctTCGGAACAGAATAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCttttt                                                                                                                                         <  2:1055777/150‑1 (MQ=255)
            ttCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                        <  1:226276/150‑1 (MQ=255)
              cGGACCAGAACAAACAACAAGTCGCGCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTtgt                                                                                                                                      <  2:519/150‑1 (MQ=255)
               ggaacAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTtgtg                                                                                                                                     <  2:1055149/146‑1 (MQ=255)
                     gAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                        <  2:405383/141‑1 (MQ=255)
                     gAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                        >  1:405383/1‑141 (MQ=255)
                      aacaaacaACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGaa                                                                                                                              >  2:120905/1‑150 (MQ=255)
                        caaacaaCAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAAtt                                                                                                                            >  2:673240/1‑150 (MQ=255)
                        caaacaaCAAGCCGCTCAAGAGCAAAAGAAAGAAAGCAAGGAGGGTACGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAAtt                                                                                                                            <  1:496081/150‑1 (MQ=255)
                          aacaacAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCtttt                                                                                                                                          <  1:198879/134‑1 (MQ=255)
                          aacaacAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCtttt                                                                                                                                          >  2:198879/1‑134 (MQ=255)
                          aacaacAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTa                                                                                                                          <  2:1010990/150‑1 (MQ=255)
                           acaacaAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAg                                                                                                                         >  2:1112205/1‑150 (MQ=255)
                            caacaaGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGc                                                                                                                        >  1:994604/1‑150 (MQ=255)
                                          agagCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGa                                                                                                          <  2:1146393/150‑1 (MQ=255)
                                                   gaaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTcc                                                                                                 <  1:1112205/150‑1 (MQ=255)
                                                        aaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCAtt                                                                                            <  1:277661/150‑1 (MQ=255)
                                                             aaGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCttt                                                                                       >  1:521169/1‑150 (MQ=255)
                                                             aaGAATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCCGCAAAAAAGCAAGGGTAAAACCTGGGTTTTTTGTGCTTAGAATTTAGCCTTCGGTCTTGCGAGGCGTTTCCGCATTGCttt                                                                                       >  2:40119/1‑150 (MQ=255)
                                                                gATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTtctc                                                                                    <  1:1263657/150‑1 (MQ=255)
                                                                         gAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGAc                                                                           >  2:828078/1‑150 (MQ=255)
                                                                         gAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGAc                                                                           >  1:1233585/1‑150 (MQ=255)
                                                                               gCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCgg                                                                     <  2:1175467/150‑1 (MQ=255)
                                                                                ccGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGc                                                                                                       <  1:138964/115‑1 (MQ=255)
                                                                                ccGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGc                                                                                                       >  2:138964/1‑115 (MQ=255)
                                                                                ccGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTtc                                                                                      <  1:232189/132‑1 (MQ=255)
                                                                                ccGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTtc                                                                                      >  2:232189/1‑132 (MQ=255)
                                                                                        tAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTtgtg                                                            <  2:994604/150‑1 (MQ=255)
                                                                                            cggcggTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCggg                                                        <  1:120905/150‑1 (MQ=255)
                                                                                            ccgcggTTAAATAGGACTTCAAATGACGAAAGAAAGAGACAAAAAAGCCAGGGTGAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCggg                                                        <  2:864507/148‑1 (MQ=255)
                                                                                             ggcggTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCgggg                                                       <  2:1265135/150‑1 (MQ=255)
                                                                                              gcggTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGc                                                      <  2:521169/150‑1 (MQ=255)
                                                                                                     aaTAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCa                                               >  2:1154339/1‑150 (MQ=255)
                                                                                                     aaTAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCa                                               >  2:70983/1‑150 (MQ=255)
                                                                                                       tAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAAt                                             <  2:1233585/150‑1 (MQ=255)
                                                                                                          gACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGc                                          >  1:440473/1‑150 (MQ=255)
                                                                                                          gACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGACACGCCCTCGGCACCTTGTGCGGGGCGTTTCTAATGGc                                          >  1:67269/1‑150 (MQ=255)
                                                                                                            cTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGaac                                        >  1:135625/1‑150 (MQ=255)
                                                                                                             ttCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGcaca                                       >  1:969369/1‑150 (MQ=255)
                                                                                                                       gaaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctg                             >  1:6936/1‑150 (MQ=255)
                                                                                                                        aaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGCCACAGCAGGAgctgg                            >  2:119453/1‑150 (MQ=255)
                                                                                                                           gaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctggctg                         <  2:440473/150‑1 (MQ=255)
                                                                                                                           gaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctggctg                         >  2:442496/1‑150 (MQ=255)
                                                                                                                               tCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGt                     <  1:442496/150‑1 (MQ=255)
                                                                                                                               tCGACAAAAAACCCACGGTAAAACCTGGCGTTTTTGTGCTTGGAATTTAGCCTTCGCTCTTGCCAGGCCGCTCCGCACTTCCTTCTCCCGCCATACCCTCGGCACCTTGTGCGGGCCGTTTCCAATGGCACAGCACGAGCTGGCTGGAGt                     >  2:1200547/1‑150 (MQ=255)
                                                                                                                                cGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGtt                    <  1:828078/150‑1 (MQ=255)
                                                                                                                                 gACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTc                   >  1:216826/1‑150 (MQ=255)
                                                                                                                                  aCAAAAAAGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCg                  >  1:1242354/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTc              <  1:7154/150‑1 (MQ=255)
                                                                                                                                         aGCCAGGGTAAAACCTGGCTTTTTTGTGGTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGGCCCTCGGCACCTTGTGCGGGGCGTTTCCAAGGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGt           <  1:451317/150‑1 (MQ=255)
                                                                                                                                                ttAAAACCTCGCTTTTTTGTTCTTGGATTTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGAACCTCGGCACCTTGTGCGGCGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTTGCTTCGGTGATCTAt    <  1:40119/149‑1 (MQ=255)
                                                                                                                                                  aaaaCCTGATTTGTTGGTGCTGGGAATTTAGCCTTCGTTCTTGAGAGTCATTTCCGCATTGCTGTCTCCGGCCCGTCGACCGGCACCTTGTGCGGGGCGTTTCCAAGGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGAtgt  <  1:927403/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GAGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGATCGATGT  >  NC_002947/36025‑36320

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: