Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 900,535 C→G 11.1% intergenic (+27/+18) PP_0784 → / ← cysZ hypothetical protein/sulfate transport protein CysZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947900,5350CG11.1% 141.4 / 8.4 53intergenic (+27/+18)PP_0784/cysZhypothetical protein/sulfate transport protein CysZ
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base G (4/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCTT  >  NC_002947/900387‑900681
                                                                                                                                                    |                                                                                                                                                   
cAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGccc                                                                                                                                                     >  4:225699/1‑149 (MQ=255)
cAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGCGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGccc                                                                                                                                                     >  3:236785/1‑149 (MQ=255)
              ggcggGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCt                                                                                                                                      >  4:172701/1‑149 (MQ=255)
                   gTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTc                                                                                                                                 <  7:261599/149‑1 (MQ=255)
                    tCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGACCCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCtt                                                                                                                                  >  6:143856/1‑147 (MQ=255)
                    acAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGGTTTCTCAGAATGCCACCCTCGAACCCTACCGTTACTCCGCCTTCTAAACAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCGCGCTCCtt                                                                                                                                  >  4:284387/2‑147 (MQ=255)
                       aGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGGACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTc                                                                                                                                     <  1:48442/141‑1 (MQ=255)
                       aGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGGACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTc                                                                                                                                     >  2:48442/1‑141 (MQ=255)
                        gCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCATTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGCGCTGCCCTGCAGCCCCTCGCTCCTTCAgctc                                                                                                                            >  1:262502/1‑149 (MQ=255)
                           aCCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAgctccgc                                                                                                                         >  2:35880/1‑149 (MQ=255)
                                gCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGc                                                                                                                                               >  8:136761/1‑123 (MQ=255)
                                    aCCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCAcc                                                                                                                <  3:225699/149‑1 (MQ=255)
                                    aCCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCAc                                                                                                                 <  3:172701/148‑1 (MQ=255)
                                     ccTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAgctccgctc                                                                                                                       >  1:329869/1‑141 (MQ=255)
                                      cTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCAcc                                                                                                                >  5:125392/1‑147 (MQ=255)
                                      cTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCa                                                                                                              <  6:347079/149‑1 (MQ=255)
                                      cTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCAGGGCAGCCCCtttt                                                                                                                                        <  4:90959/124‑4 (MQ=255)
                                      cTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCAGGGCAGCCCCtttt                                                                                                                                        >  3:90959/1‑121 (MQ=255)
                                       tCGTGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCTAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTccc                                                                                                                     <  6:125173/141‑1 (MQ=255)
                                       tCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCgacatcg                                                                                                              >  5:125173/1‑142 (MQ=255)
                                          agaAGGTTAAGGCAGCGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGaac                                                                                                          <  4:236785/149‑1 (MQ=255)
                                            aaGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTg                                                                                                                                                   <  8:254108/107‑1 (MQ=255)
                                              gtttAAGGAAGTGTAACATCTAATACTGCTTTCTCAGATTGCCACCCTCTAACCCTACAGTTTCTCCTACTTCTAAAAAGAAATCCCCAAAAAGGGGCTCCCGTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGaacaac                                                                                                       <  3:321753/146‑1 (MQ=255)
                                                     gCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGCGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGcccc                                                                                               <  2:262502/149‑1 (MQ=255)
                                                            aaCATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCTAGTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGATCGGCTCCGGCACCCAGAACAACTTAGCCACGGcacat                                                                                        >  8:88289/1‑145 (MQ=255)
                                                                tCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTgccg                                                                                    >  7:14326/1‑149 (MQ=255)
                                                                    aTACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCgg                                                                                <  2:111957/149‑1 (MQ=255)
                                                                             tctcAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCAtcatca                                                                        <  6:284204/148‑1 (MQ=255)
                                                                               tcAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCa                                                                     >  8:193416/1‑149 (MQ=255)
                                                                                    ttGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGtt                                                                 >  1:221269/1‑148 (MQ=255)
                                                                                       ccACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCAc                                                              <  8:14326/148‑1 (MQ=255)
                                                                                                   ccTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGTCATCATCAGCACGTTCACGCCTGGAatcatc                                                 <  1:35880/149‑1 (MQ=255)
                                                                                                   ccTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAatcat                                                  >  4:127327/1‑148 (MQ=255)
                                                                                                    cTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAAtcatca                                                <  7:88289/149‑1 (MQ=255)
                                                                                                     tACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAg                                               >  7:212061/1‑149 (MQ=255)
                                                                                                         gTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAgcgc                                            >  8:149634/1‑148 (MQ=255)
                                                                                                              tcCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCagg                                        >  2:124479/1‑147 (MQ=255)
                                                                                                                 gACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCaggtaggt                                   >  1:150682/1‑149 (MQ=255)
                                                                                                                    ttCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTGCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGa                                 >  1:68140/1‑149 (MQ=255)
                                                                                                                    ttCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAAtcatca                                                <  5:125993/133‑1 (MQ=255)
                                                                                                                    ttCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAAtcatca                                                >  6:125993/1‑133 (MQ=255)
                                                                                                                         aaaaGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgccgc                           >  1:23836/1‑149 (MQ=255)
                                                                                                                         aaaaGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgccgc                           <  3:127327/149‑1 (MQ=255)
                                                                                                                         aaaaGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCCCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGCAATCATCAGCGCCAGCTAGGAGATgccgc                           >  4:216238/1‑149 (MQ=255)
                                                                                                                         aaaaGCAATCCGAAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATgccgc                           <  3:176710/149‑1 (MQ=255)
                                                                                                                               aaTCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatcat                                                                          >  7:267547/1‑96 (MQ=255)
                                                                                                                               aaTCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGcatcat                                                                          <  8:267547/96‑1 (MQ=255)
                                                                                                                                   cGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCa                  <  7:149634/148‑1 (MQ=255)
                                                                                                                                    gCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAgg                >  5:153432/1‑149 (MQ=255)
                                                                                                                                    gCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAgg                >  2:113044/1‑149 (MQ=255)
                                                                                                                                     cAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGc               <  2:23836/149‑1 (MQ=255)
                                                                                                                                      aaaaaGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAGAGCCTGGGct              <  8:212061/149‑1 (MQ=255)
                                                                                                                                       gcgaGGGGCTGCAGGGCAGCCCCTTGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGctc             >  7:242973/4‑149 (MQ=255)
                                                                                                                                       gcgaGGGGCTGCAGGGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCGGGACTCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGcgc             >  4:28945/4‑147 (MQ=255)
                                                                                                                                          aGGGGCTGCAGGGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGc          >  3:108043/1‑149 (MQ=255)
                                                                                                                                                  cccTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCtt  <  1:124479/149‑1 (MQ=255)
                                                                                                                                                  cccTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCtt  <  1:113044/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCCTTGGCGTTCGGCGGGTCAAGCGACCTTGCTAACCTCGAGAAGGTTAAGGCAGTGGAACATCTAATACTGCTTTCTCAGATTGCCACCCTCGAACCCTACAGTTTCTCCGACTTCTAAAAAGCAATCCGCAAAAAGGGGCTGCCCTGCAGCCCCTCGCTCCTTCAGCTCCGCTCCCGCACCCAGAACAACGTAGCCCCGGCCACTGCCGCCGGCATCATCAGCACGTTCACGCCTGGAATCATCAGCGCCAGGTAGGTGATGCCGCCAAAGCCCAGGCTCTGCCAGCGCTT  >  NC_002947/900387‑900681

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