Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,389,508 T→A 10.8% T411T (ACA→ACT)  PP_3865 ← tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,389,5080TA10.8% 152.7 / 2.7 57T411T (ACA→ACT) PP_3865tail protein
Reads supporting (aligned to +/- strand):  ref base T (20/30);  new base A (3/3);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGC  >  NC_002947/4389361‑4389650
                                                                                                                                                   |                                                                                                                                              
ggCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGACCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACtg                                                                                                                                               >  8:110207/1‑149 (MQ=255)
   cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGGGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                            <  6:1311/149‑1 (MQ=255)
   cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGCTAAGGCCAATCGGCGTGGCGGCGAAGGCCCCACCGCGCACCCGCGGGGGCCCTGTgg                                                                                                                                            >  3:262452/1‑149 (MQ=255)
   cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                            >  2:268833/1‑149 (MQ=255)
   cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACAGTgg                                                                                                                                            <  8:252324/149‑1 (MQ=255)
       aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgc                                                                                                                                        <  4:122275/149‑1 (MQ=255)
       aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcg                                                                                                                                         <  1:268833/148‑1 (MQ=255)
          aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGaa                                                                                                                                     <  2:299068/149‑1 (MQ=255)
          aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGAGCAGGCCAAACGGCGTGGCCGCCAAGGCCACACCCAGCACCCCCGCGGCCCCCGGGGCcccgcg                                                                                                                                     >  4:221872/1‑147 (MQ=255)
          aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCCCGCTCCCCGCGGCGAGCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCCCCCCCGCCGCCACACTGGCcacgcg                                                                                                                                     >  4:217434/1‑147 (MQ=255)
           gCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgcgg                                                                                                                                    <  3:217434/149‑4 (MQ=255)
           gCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgcgg                                                                                                                                    <  3:221872/149‑4 (MQ=255)
            cTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACAATCAGACCGGCCGCCCCCCCCAAACCCCCGGCGAAAAGGCCAAAACGCGTGGCCGCGAACGCCACACCGCGCACCCCCCTTGCCAAAGTGGCCGCGAaca                                                                                                                                   >  8:180081/1‑149 (MQ=255)
                 acaGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGGCGCGTTCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGGCACTGTGGCCGCGAACAccgcc                                                                                                                              <  6:262979/149‑1 (MQ=255)
                    gggCGCGGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACAATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGGGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCcgca                                                                                                                           <  7:180081/149‑1 (MQ=255)
                          tGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGCTCCCCGCGGCGATCAGGCCAAACGGCCTGGCCGCGAAGGCCCCACCCAGCCCCCGCCTGGCCACTGTGGCCCCGAACACCGCCGCACGCAgc                                                                                                                     >  7:9985/1‑149 (MQ=255)
                          tGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCCATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCCCGAACACCGCCGCACGCAgc                                                                                                                     >  7:4735/1‑149 (MQ=255)
                            aagaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACAGCCGCACGCAgcgc                                                                                                                   >  5:360166/1‑149 (MQ=255)
                                   gCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAg                                                                                                            <  3:25984/149‑1 (MQ=255)
                                    cccGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAgg                                                                                                           >  2:58545/1‑149 (MQ=255)
                                          cTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCCCTCCCCCCGGCGCTCAGGCCAATCGGCGTGGCCGCGCAGGCCACACCGAGCACCCCCGTGGCCACAGTGGCCGCGAACACCCCCGCACGCAGCGCGGTGAAGGCCGccc                                                                                                     >  3:69933/1‑149 (MQ=255)
                                                ccAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCaa                                                                                               <  1:58545/149‑1 (MQ=255)
                                                 cAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATACCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAt                                                                                              >  6:74682/1‑149 (MQ=255)
                                                            gctCAGACCGCCCGCCGCCGCGATCCCCCCGGCGATCACGCCAATCGGCCTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCcgc                                                                                   <  4:262452/147‑1 (MQ=255)
                                                            gaTCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAac                                                                                                                                    <  5:294840/100‑1 (MQ=255)
                                                                aGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCa                                                                                <  3:166852/148‑1 (MQ=255)
                                                                 gACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCAc                                                                               >  1:231495/1‑148 (MQ=255)
                                                                         gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                      <  5:314556/149‑1 (MQ=255)
                                                                         gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                      <  8:97577/149‑1 (MQ=255)
                                                                         gccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                      >  6:314556/1‑149 (MQ=255)
                                                                            gccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                   <  1:117682/149‑1 (MQ=255)
                                                                                     cccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                                                           <  4:69933/148‑1 (MQ=255)
                                                                                      tcgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatc                                                         <  7:210386/148‑1 (MQ=255)
                                                                                      tcgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatc                                                         <  7:204958/148‑1 (MQ=255)
                                                                                      ccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGat                                                          >  4:305552/1‑148 (MQ=255)
                                                                                       cgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatca                                                        <  3:305552/149‑1 (MQ=255)
                                                                                          ggcgATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCt                                                     <  6:360166/149‑1 (MQ=255)
                                                                                            cgATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTgc                                                   <  7:128396/149‑1 (MQ=255)
                                                                                             gATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTgcg                                                  <  7:110207/149‑1 (MQ=255)
                                                                                                   gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat                                            >  6:324373/1‑149 (MQ=255)
                                                                                                    ccAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTg                                                            >  5:67584/1‑132 (MQ=255)
                                                                                                    ccAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTg                                                            <  6:67584/132‑1 (MQ=255)
                                                                                                        tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACg                                                                                          <  8:76851/98‑1 (MQ=255)
                                                                                                        tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACg                                                                                          >  7:76851/1‑98 (MQ=255)
                                                                                                          ggCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGAt                                     <  5:156450/149‑1 (MQ=255)
                                                                                                            cGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATcc                                   <  5:74682/149‑1 (MQ=255)
                                                                                                              tGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                 <  2:231495/149‑1 (MQ=255)
                                                                                                                  cgcgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGcc                             >  4:353607/1‑149 (MQ=255)
                                                                                                                       aGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGt                        <  8:4735/149‑1 (MQ=255)
                                                                                                                       aGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGt                        <  8:9985/149‑1 (MQ=255)
                                                                                                                          ccACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGc                     >  5:17666/1‑149 (MQ=255)
                                                                                                                           cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGAt                                     <  4:136293/132‑1 (MQ=255)
                                                                                                                              acCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCTCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                   >  2:53791/1‑147 (MQ=255)
                                                                                                                                cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                 >  5:248510/1‑147 (MQ=255)
                                                                                                                                cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                 <  6:248510/147‑1 (MQ=255)
                                                                                                                                       ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAt        >  1:177259/1‑149 (MQ=255)
                                                                                                                                          cgTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGa     >  4:205221/1‑149 (MQ=255)
                                                                                                                                             ggCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGTTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGc  >  3:268468/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GGCCTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGC  >  NC_002947/4389361‑4389650

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: ‑