Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,105,532 Δ1 bp 7.4% coding (20/702 nt) arsH ← ArsH protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,105,5320T. 7.4% 163.5 / 7.4 54coding (20/702 nt)arsHArsH protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base . (2/2);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

AGTGCTCCAGGAGGCGCGCGGCTTCTTCTACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGG  >  NC_002947/3105386‑3105679
                                                                                                                                                  |                                                                                                                                                   
aGTGCTCCAGGAGGCGCGCGGCTTCTTCTACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTa                                                                                                                                                   <  2:129260/149‑1 (MQ=255)
    cTCCAGGAGGCGCGCGGCTTCTTCTACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCa                                                                                                                                               <  2:114173/149‑1 (MQ=255)
               gcgcgGCTTCTTCTACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGa                                                                                                                                     >  1:80394/1‑148 (MQ=255)
                       tcttctACCAGCAAGCGACTGAAGGAGAGTTCACGCCTAGACCCGTACAGCAGCAGCATCAGAGGTATGAGTTCGAAGGACGGGGTTGCCGAAGGCGGCAGGTCGAACAGCTTGAGGTCGAGGTTAGCCACTTGCTCGGACATGattac                                                                                                                            >  7:341307/1‑147 (MQ=255)
                         ttctACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTttctt                                                                                                                          >  5:186208/1‑149 (MQ=255)
                         ttctACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTttctt                                                                                                                          >  6:291286/1‑149 (MQ=255)
                              ccAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCt                                                                                                                       <  5:75157/147‑1 (MQ=255)
                              ccAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCt                                                                                                                       <  3:30104/147‑1 (MQ=255)
                               cAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGtt                                                                                                                    >  5:66507/1‑149 (MQ=255)
                               cAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGtt                                                                                                                    <  3:5939/149‑1 (MQ=255)
                                 gCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCa                                                                                                                  <  3:62644/149‑1 (MQ=255)
                                   aaGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGtt                                                                                                                    >  5:238054/1‑145 (MQ=255)
                                   aaGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGtt                                                                                                                    <  6:238054/145‑1 (MQ=255)
                                            aaGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAgg                                                                                                                                                 >  2:127279/1‑107 (MQ=255)
                                            aaGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGa                                                                                                       >  7:139781/1‑149 (MQ=255)
                                            aaGGAGCGTTCTCGGGTAGACCCCTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGa                                                                                                       >  5:397813/1‑149 (MQ=255)
                                              ggAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATg                                                                                                     <  6:172235/149‑1 (MQ=255)
                                                      ctcGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAg                                                                                             >  1:117962/1‑149 (MQ=255)
                                                       tcGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAgc                                                                                            >  6:192234/1‑149 (MQ=255)
                                                         gggTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAgc                                                                                            <  5:336235/147‑1 (MQ=255)
                                                           gTAGTCCCGTACAGCAGCAGGATGCTATGTTTGTGTTCGCAGGTCGGGGTTGCCTAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAATCGAGCCGATGCGATCCAGCGCGt                                                                                        <  7:381577/149‑1 (MQ=255)
                                                           gTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGt                                                                                        >  3:377083/1‑149 (MQ=255)
                                                              gACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGc                                                                                      <  3:286599/148‑1 (MQ=255)
                                                                        gcagcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTc                                                                           >  4:366470/1‑149 (MQ=255)
                                                                         cagcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTcg                                                                          >  8:206368/1‑149 (MQ=255)
                                                                         cagcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTcgc                                                                         >  3:355842/1‑149 (MQ=255)
                                                                         cagcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTcgc                                                                         >  3:361688/1‑149 (MQ=255)
                                                                          agcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATcc                                                                                               >  3:145094/1‑127 (MQ=255)
                                                                          agcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATcc                                                                                               <  4:145094/127‑1 (MQ=255)
                                                                          agcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCtt                                                                                    >  5:116412/1‑138 (MQ=255)
                                                                          agcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCtt                                                                                    <  6:116412/138‑1 (MQ=255)
                                                                          agcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTcgcg                                                                        <  4:361688/149‑1 (MQ=255)
                                                                          agcagGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTcgcg                                                                        <  4:355842/149‑1 (MQ=255)
                                                                              ggATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTc                                                                     <  5:192234/149‑1 (MQ=255)
                                                                               gATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGt                                                                                        <  1:161427/129‑1 (MQ=255)
                                                                                aTGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGa                                                                   <  8:139781/149‑1 (MQ=255)
                                                                                   cGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAgg                                                                 <  4:377083/148‑1 (MQ=255)
                                                                                    gAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGt                                                                >  3:178412/1‑148 (MQ=255)
                                                                                    gAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTg                                                               <  6:397813/149‑1 (MQ=255)
                                                                                      ggTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGAt                                                             >  7:296911/1‑149 (MQ=255)
                                                                                            tgtTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTTATCGAAc                                                        <  5:291286/148‑1 (MQ=255)
                                                                                                        ggggTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCgaag                                            <  1:43689/148‑1 (MQ=255)
                                                                                                          ggTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGaagaag                                         >  2:260556/1‑149 (MQ=255)
                                                                                                          ggTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCAAGGTGATCGAACGGCAGCGCGaagaag                                         >  1:393895/1‑149 (MQ=255)
                                                                                                                cGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTgg                                   <  4:83945/149‑1 (MQ=255)
                                                                                                                 gAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCgaagaa                                          >  5:162935/1‑141 (MQ=255)
                                                                                                                 gAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCgaagaa                                          <  6:162935/141‑1 (MQ=255)
                                                                                                                   aGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAg                                <  2:117962/149‑1 (MQ=255)
                                                                                                                            gTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCt                                     <  2:217801/135‑1 (MQ=255)
                                                                                                                            gTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCt                                     >  1:217801/1‑135 (MQ=255)
                                                                                                                            gTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTc                        >  3:198204/1‑148 (MQ=255)
                                                                                                                                        gAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAgt           >  4:168629/1‑149 (MQ=255)
                                                                                                                                         aGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGTAGGCCTGGCAGCGAGACTCGATACGGGCCAgtg          <  1:260556/149‑1 (MQ=255)
                                                                                                                                            tCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCCTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTgg       >  7:377329/1‑149 (MQ=255)
                                                                                                                                                  ttAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAAAGGCAGCGCGAAGAAGGCCTGGCAGCGAGCCTCGATACGGGCCAGTGTGGCCCgg  >  8:375738/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AGTGCTCCAGGAGGCGCGCGGCTTCTTCTACCAGCAAGCGACTGAAGGAGCGTTCTCGGGTAGACCCGTACAGCAGCAGGATGCGAGGTTTGTGTTCGCAGGTCGGGGTTGCCGAAGGCGGCAGGTCGAACAGCGTGAGGTCGAGGTTAGGCAGTTGCTCGGACATGTTTTCTTCCTGTTCAAAGCGAGCCGATGCGATCCAGCGCGTGCTTGAGCTGCTCGCGGTCGAGGTGATCGAACGGCAGCGCGAAGAAGGCCTGGCAGCGAGACTCGATACGGGCCAGTGTGGCCCGG  >  NC_002947/3105386‑3105679

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: