Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,928,591 T→C 6.5% M515V (ATG→GTG)  PP_1728 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,928,5910TC 6.5% 180.3 / 2.2 63M515V (ATG→GTG) PP_1728membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base C (2/2);  total (31/32)
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.18e-01

TCCACATAGTCGTCGATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCATGGTGTACGGCGGTG  >  NC_002947/1928452‑1928733
                                                                                                                                           |                                                                                                                                              
tCCACATAGTCGTCGATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTg                                                                                                                                      <  1:1169128/150‑1 (MQ=255)
   acaTAGTCGTCGATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGTTGACGTTCAGGCCCATGTTGGcc                                                                                                                                   >  2:1231866/1‑150 (MQ=255)
       aGTCGTCGATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCAAGGCCGATCACACCGGCACCGGCCAGCAGTGCGGTGACGTTCATGCCCATGTTGGCCAggg                                                                                                                               <  1:997553/150‑1 (MQ=255)
         tcgtcgATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGcg                                                                                                                             >  2:149765/1‑150 (MQ=255)
          cgtcgATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGcga                                                                                                                            <  2:1172139/150‑1 (MQ=255)
                tGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGa                                                                                                                      <  2:1251898/150‑1 (MQ=255)
                     aGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACcgcg                                                                                                                 <  1:469204/150‑1 (MQ=255)
                      gCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCga                                                                                                                >  1:1203642/1‑150 (MQ=255)
                       cGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCgat                                                                                                               <  2:988963/150‑1 (MQ=255)
                                 atgatgatgAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACAACCCCACCGGCAGCCTACGTTGTTGCGGTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATg                                                                                                     >  1:988010/1‑150 (MQ=255)
                                  tgatgatgaACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATTACCGCGATGATGAAGATGa                                                                                                    <  2:1076986/150‑1 (MQ=255)
                                     tgatgaACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACTCCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCa                                                                                                 >  2:906717/1‑150 (MQ=255)
                                          aaCAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGACCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGACCATGTTGGCCAGGGCAACGATGACCGCGATGATGAAGATGACCACGAAc                                                                                            >  2:135209/1‑150 (MQ=255)
                                            cAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACAt                                                                                          >  1:108045/1‑150 (MQ=255)
                                             aGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATTCTATTGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATc                                                                                         <  1:732664/150‑1 (MQ=255)
                                             aGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAAGCAGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGAGCGACGATGACCGCGATGATGAAGATGACCACGAACATc                                                                                         >  2:633051/1‑150 (MQ=255)
                                              ggCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCa                                                                                        >  1:640993/1‑150 (MQ=255)
                                               gCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCAc                                                                                       <  2:1115458/150‑1 (MQ=255)
                                                 cGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGt                                                                                     <  1:705504/150‑1 (MQ=255)
                                                  ggTGATCAGGTCGGCCACCAGCGACTACGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGa                                                                                                    <  1:267259/134‑1 (MQ=255)
                                                  ggTGATCAGGTCGGCCACCAGCGACTACGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGa                                                                                                    >  2:267259/1‑134 (MQ=255)
                                                       tCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGaa                                                                               <  2:1203642/150‑1 (MQ=255)
                                                        cAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAAt                                                                              <  2:640993/150‑1 (MQ=255)
                                                          ggTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAAGTTACCACGAACATCACGTTGCGAATCa                                                                            >  2:1126212/1‑150 (MQ=255)
                                                               gCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAAc                                                                          <  2:817059/147‑1 (MQ=255)
                                                               gCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAAc                                                                          >  1:817059/1‑147 (MQ=255)
                                                                 caccaGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCACGCCCAt                                                                                                                                          <  1:523536/81‑1 (MQ=255)
                                                                 caccaGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCACGCCCAt                                                                                                                                          >  2:523327/1‑81 (MQ=255)
                                                                 caccaGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCACGCCCAt                                                                                                                                          >  2:523536/1‑81 (MQ=255)
                                                                 caccaGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCACGCCCAt                                                                                                                                          <  1:523327/81‑1 (MQ=255)
                                                                    caGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGcatcat                                                                  >  1:1152350/1‑150 (MQ=255)
                                                                      gCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATgg                                                                >  1:593147/1‑150 (MQ=255)
                                                                         aCTGCGCGCCAAAACCGATGGCCAGGCCGATCACCACGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGttt                                                             <  1:1231866/150‑1 (MQ=255)
                                                                           tGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTgc                                                           <  1:693441/150‑1 (MQ=255)
                                                                              gcgcCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTgc                                                           >  1:656706/1‑147 (MQ=255)
                                                                              gcgcCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTgc                                                           <  2:656706/147‑1 (MQ=255)
                                                                               cgcCAAAACCGATGGCCAGGCCGATCACTCCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCAc                                                       <  1:906717/150‑1 (MQ=255)
                                                                                gcCAAAACCGATGGCCAGGCCGATCACACCGGCACCGTCCAGCAGCGGGGTGACGTTCAt                                                                                                                                                <  1:670533/60‑1 (MQ=255)
                                                                                gcCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCAt                                                                                                                                                >  2:670533/1‑60 (MQ=255)
                                                                                   aaaaCCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGc                                                   >  1:1139582/1‑150 (MQ=255)
                                                                                   aaaaCCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGc                                                   <  2:988010/150‑1 (MQ=255)
                                                                                        cGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGttcg                                              >  2:110531/1‑148 (MQ=255)
                                                                                          aTGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACgg                                                                        >  1:22146/1‑122 (MQ=255)
                                                                                          aTGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACgg                                                                        <  2:22146/122‑1 (MQ=255)
                                                                                          aTGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGcc                                            >  2:762174/1‑150 (MQ=255)
                                                                                             ccctcgacGACCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCAGGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGACCCAACGGCAGCATGGTTTGCGCACGGGCGTTGGCCAGc                                         <  2:593147/143‑1 (MQ=255)
                                                                                               cAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGccc                                       >  1:827950/1‑150 (MQ=255)
                                                                                                aGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGcccc                                      >  2:696807/1‑150 (MQ=255)
                                                                                                 ggCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGccccc                                     >  1:1271591/1‑150 (MQ=255)
                                                                                                         acacCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGAt                                                                                                                        >  1:1027654/1‑59 (MQ=255)
                                                                                                         acacCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGAt                                                                                                                        <  2:1027654/59‑1 (MQ=255)
                                                                                                         acacCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATTACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAg                             >  1:596770/1‑150 (MQ=255)
                                                                                                           accggcaccggcCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTg                                                                                   >  1:608609/1‑94 (MQ=255)
                                                                                                           accggcaccggcCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTg                                                                                   <  2:608609/94‑1 (MQ=255)
                                                                                                               gcaccggcCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCgatgat                                                                                                            >  1:51959/1‑65 (MQ=255)
                                                                                                               gcaccggcCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCgatgat                                                                                                            <  2:51959/65‑1 (MQ=255)
                                                                                                                 accggcCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAg                     <  1:653172/150‑1 (MQ=255)
                                                                                                                         gcagcGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCATGgt             <  2:662803/150‑1 (MQ=255)
                                                                                                                            gcGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCAt                >  1:874999/1‑144 (MQ=255)
                                                                                                                            gcGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCAt                <  2:874999/144‑1 (MQ=255)
                                                                                                                             cGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCATGGTGTAc         <  2:1152350/150‑1 (MQ=255)
                                                                                                                                   gACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCATGGTGTACGGCGgt   <  1:149765/150‑1 (MQ=255)
                                                                                                                                    aCGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGAGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCATGGTGTACGGCGgtg  >  1:1029545/1‑150 (MQ=255)
                                                                                                                                      gTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGc                  >  2:251850/1‑132 (MQ=255)
                                                                                                                                      gTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGc                  <  1:251850/132‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
TCCACATAGTCGTCGATGGCCAGCGAATCTTCGATGATGATGAACAGGCCGGTGATCAGGTCGGCCACCAGCGACTGCGCGCCAAAACCGATGGCCAGGCCGATCACACCGGCACCGGCCAGCAGCGGGGTGACGTTCATGCCCATGTTGGCCAGGGCGACGATGACCGCGATGATGAAGATGACCACGAACATCACGTTGCGAATCAACGGCATCATGGTTTGCGCACGGGCGTTGGCCAGCCCCCGGCGCGAGCGCACCAGGGCATGGTGTACGGCGGTG  >  NC_002947/1928452‑1928733

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