Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,727,120:1 +T 10.5% coding (333/645 nt) PP_3292 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,727,1201.T10.5% 200.7 / 20.9 57coding (333/645 nt)PP_3292membrane protein
Reads supporting (aligned to +/- strand):  ref base . (22/29);  new base T (2/4);  total (24/33)
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.68e-01

CGCTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGC  >  NC_002947/3726972‑3727261
                                                                                                                                                     |                                                                                                                                             
cgcTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTg                                                                                                                                                >  6:121115/1‑149 (MQ=255)
cgcTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTg                                                                                                                                                <  1:96456/149‑1 (MQ=255)
           gCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTggg                                                                                                                                    <  6:13611/149‑1 (MQ=255)
           gCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTggg                                                                                                                                    <  2:289529/149‑1 (MQ=255)
             ggTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGca                                                                                                                                  <  1:243834/149‑1 (MQ=255)
             ggTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGca                                                                                                                                  <  5:121115/149‑1 (MQ=255)
             ggTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGca                                                                                                                                  <  1:122310/149‑1 (MQ=255)
                    cGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGCAc                                                                                                                           <  4:148988/149‑1 (MQ=255)
                    cGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGCAc                                                                                                                           <  8:253741/149‑1 (MQ=255)
                     gaggCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCAcc                                                                                                                          <  1:50673/147‑1 (MQ=255)
                        gCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCttt                                                                                                                       <  6:249791/149‑1 (MQ=255)
                            gAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGc                                                                                                                             <  7:285123/139‑1 (MQ=255)
                            gAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGCCAGGc                                                                                                                             >  8:285123/1‑139 (MQ=255)
                             cgctgcACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTtggct                                                                                                                  <  5:125266/145‑1 (MQ=255)
                                           cAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGAGTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGcac                                                                                                                                 >  7:209830/1‑120 (MQ=255)
                                           cAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGcac                                                                                                                                 <  8:209830/120‑1 (MQ=255)
                                            aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGt                                                                                                   >  3:30089/1‑149 (MQ=255)
                                            aaaGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGt                                                                                                   >  5:58047/1‑149 (MQ=255)
                                               gCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCGCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCCCCTTTTGGCTGGCAGGGCGCTCAGTGAc                                                                                                >  5:160149/1‑149 (MQ=255)
                                                  aTTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACg                                                                                               >  6:310760/1‑147 (MQ=255)
                                                        aaCGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGt                                                                                       >  8:258250/1‑149 (MQ=255)
                                                                  agaCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGc                                                                             <  5:310760/149‑1 (MQ=255)
                                                                    aCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGc                                                                             >  1:67064/1‑147 (MQ=255)
                                                                    aCAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGa                                                                           <  7:258250/149‑1 (MQ=255)
                                                                     cAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAg                                                                          <  8:209923/149‑1 (MQ=255)
                                                                      aGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCtt                                                                                                                        <  4:328821/102‑1 (MQ=255)
                                                                      aGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCtt                                                                                                                        >  3:328821/1‑102 (MQ=255)
                                                                        gCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGcc                                                                       <  2:67064/149‑1 (MQ=255)
                                                                          tGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCaa                                                                     >  1:84193/1‑149 (MQ=255)
                                                                           ggCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCa                                                                                                     <  6:268728/116‑1 (MQ=255)
                                                                           ggCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCa                                                                                                     >  5:268728/1‑116 (MQ=255)
                                                                            gCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAAcc                                                                   >  7:268452/1‑149 (MQ=255)
                                                                            gCAGGCCAGACTGGCGTTCATTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGATGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAAATGCGAGGCCAAcc                                                                   >  2:259986/1‑149 (MQ=255)
                                                                             cAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCt                                                                  <  2:84193/149‑1 (MQ=255)
                                                                               ggCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGc                                                                <  6:58047/149‑1 (MQ=255)
                                                                                gCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGc                                                                                                                                   <  4:180134/82‑1 (MQ=255)
                                                                                gCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGTGCCCAGCAGGGc                                                                                                                                   >  3:180134/1‑82 (MQ=255)
                                                                                gCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCa                                                               <  8:268452/149‑1 (MQ=255)
                                                                                         gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTc                                                          >  7:285849/1‑145 (MQ=255)
                                                                                         gCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTc                                                          <  8:285849/145‑1 (MQ=255)
                                                                                           gTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTgg                                                    <  4:30089/149‑1 (MQ=255)
                                                                                            ttCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGca                                                                                                                                  <  6:197667/70‑1 (MQ=255)
                                                                                            ttCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGca                                                                                                                                  >  5:197667/1‑70 (MQ=255)
                                                                                                  cGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGAtata                                               <  6:263303/147‑1 (MQ=255)
                                                                                                        gTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCg                                           >  4:33428/1‑145 (MQ=255)
                                                                                                           gCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATg                                     >  5:105939/1‑148 (MQ=255)
                                                                                                           gCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATg                                     >  4:127895/1‑148 (MQ=255)
                                                                                                               gggAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaaaa                                >  7:250300/1‑149 (MQ=255)
                                                                                                               gggAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGaaa                                  <  3:33428/147‑1 (MQ=255)
                                                                                                                   aTCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCa                            <  8:250300/149‑1 (MQ=255)
                                                                                                                       aaCATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACcgc                        >  1:282383/1‑149 (MQ=255)
                                                                                                                                ctgGTGGTGCTGTTGTGCCTGGCCCTGCTTGGCACAGGAACCTTTTGGCTGGCAGGGCGCTCAGTGAAGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGTCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTtctc                 <  1:259986/147‑1 (MQ=255)
                                                                                                                                 tggtggTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                         <  7:20851/122‑1 (MQ=255)
                                                                                                                                 tggtggTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGaa                                         >  8:20851/1‑122 (MQ=255)
                                                                                                                                        gCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCgg       >  3:175287/1‑149 (MQ=255)
                                                                                                                                          tgttgtGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGgtg     <  6:105939/149‑1 (MQ=255)
                                                                                                                                           gttgtGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAgg                                                                         >  4:76604/1‑80 (MQ=255)
                                                                                                                                           gttgtGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAgg                                                                         <  3:76604/80‑1 (MQ=255)
                                                                                                                                              gtgCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGc  <  6:160149/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CGCTAGGCTGGGCGGTGCAGCGGGGCGTGAAGCGCACCGTTGGCAAAGCCATTTTCAACGGCCAGGAGACAGGCTGGCAGGCCAGACTGGCGTTCCTTCGCTCAGTGGCTCGGGAATCGAACATCCTGCTGGTGGTGCTGTTGTGCCTGGCCCTGCTGGGCACAGGCACCTTTTGGCTGGCAGGGCGCTCAGTGACGCAGATCGTGGCGAACTGCGAGGCCAACCTGCAGGCTCAAGTGGATATACCCGAAGATGAAAAACTCACCGCGCTTCTCGGCAAACCGGTGTGC  >  NC_002947/3726972‑3727261

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