Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,100,833 G→C 7.0% W236S (TGG→TCG)  PP_3608 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,100,8330GC 7.0% 177.4 / 3.9 57W236S (TGG→TCG) PP_3608LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (25/28);  new base C (2/2);  total (27/30)
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.16e-01

CTGCCTGCCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAGCCTGGC  >  NC_002947/4100685‑4100979
                                                                                                                                                    |                                                                                                                                                  
ctgcctgcCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTg                                                                                                                                                    >  3:370663/1‑149 (MQ=255)
     tgccAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTg                                                                                                                                               >  8:155403/1‑149 (MQ=255)
           tGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAg                                                                                                                                         <  3:88026/149‑1 (MQ=255)
               gATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCGCACTCCAGTGTGGCaaggcc                                                                                                                                     <  4:137645/149‑6 (MQ=255)
                  atGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGt                                                                                                                                  >  7:4446/1‑149 (MQ=255)
                    ggCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGc                                                                                                                                >  3:290954/1‑149 (MQ=255)
                               aCTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCAcc                                                                                                                     >  4:5391/1‑149 (MQ=255)
                                   gctggcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCGCACTCCAGTGTGGCaa                                                                                                                                         <  8:186515/125‑2 (MQ=255)
                                   gctggcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCGCACTCCAGTGTGGCaa                                                                                                                                         >  7:186515/1‑124 (MQ=255)
                                   gctggcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGAcc                                                                                                                 <  2:296482/149‑1 (MQ=255)
                                       gcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGa                                                                                                              <  4:96538/148‑1 (MQ=255)
                                       gcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGa                                                                                                              <  4:290954/148‑1 (MQ=255)
                                       gcACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGAt                                                                                                             >  8:59093/1‑149 (MQ=255)
                                         aCGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATcc                                                                                                           >  3:163688/1‑149 (MQ=255)
                                            cccACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGa                                                                                                        <  7:263793/149‑1 (MQ=255)
                                             ccACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGa                                                                                                        >  2:9392/1‑148 (MQ=255)
                                                 aTGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTg                                                                                                   <  6:127835/149‑1 (MQ=255)
                                                   gTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGc                                                                                                 >  3:341933/1‑149 (MQ=255)
                                                    tATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGc                                                                                                 <  7:342873/148‑1 (MQ=255)
                                                      ttGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCa                                                                                               <  4:163688/148‑1 (MQ=255)
                                                         aTGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                          <  3:137552/118‑1 (MQ=255)
                                                         aTGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                          >  4:137552/1‑118 (MQ=255)
                                                         aTGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCcg                                                                                           >  2:54309/1‑149 (MQ=255)
                                                          tGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                          >  3:140502/1‑117 (MQ=255)
                                                          tGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCg                                                                                                                          <  4:140502/117‑1 (MQ=255)
                                                              cATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGc                                                                                                                            >  6:199713/1‑111 (MQ=255)
                                                                tGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCtgtg                                                                                    <  5:39386/149‑1 (MQ=255)
                                                                  cgcCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTggg                                                                                  <  4:341933/149‑1 (MQ=255)
                                                                        gATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGtt                                                                            <  4:370663/149‑1 (MQ=255)
                                                                          tGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgc                                                                          <  1:361087/149‑1 (MQ=255)
                                                                             tttCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgcgc                                                                        >  3:22703/1‑148 (MQ=255)
                                                                              ttCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGcc                                                                       >  6:331301/1‑148 (MQ=255)
                                                                                 gATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAaca                                                                   >  3:239171/1‑149 (MQ=255)
                                                                                   tGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAg                                                                                                     >  4:201617/1‑113 (MQ=255)
                                                                                   ggAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAAcaca                                                                 <  7:338157/148‑1 (MQ=255)
                                                                                      aaTCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                                <  1:231748/147‑1 (MQ=255)
                                                                                      aaTCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                                >  2:231748/1‑147 (MQ=255)
                                                                                      aaTCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACATGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGa                                                              >  3:218885/1‑149 (MQ=255)
                                                                                              gTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgcgc                                                                        <  5:250827/131‑1 (MQ=255)
                                                                                              gTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTgcgc                                                                        >  6:250827/1‑131 (MQ=255)
                                                                                                   cTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGtt                                                 <  1:9392/149‑1 (MQ=255)
                                                                                                      gccgccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTcac                                              <  7:155403/149‑1 (MQ=255)
                                                                                                      gccgccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTca                                               <  1:341700/148‑1 (MQ=255)
                                                                                                         gccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGtt                                                 <  3:226867/143‑1 (MQ=255)
                                                                                                         gccCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGtt                                                 >  4:226867/1‑143 (MQ=255)
                                                                                                           ccGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                                >  6:32541/1‑126 (MQ=255)
                                                                                                           ccGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACaa                                                                <  5:32541/126‑1 (MQ=255)
                                                                                                                 cTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTgg                                   <  5:331301/149‑1 (MQ=255)
                                                                                                                        cgcCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATccc                            <  3:5391/149‑1 (MQ=255)
                                                                                                                         gcCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATcccc                           <  1:54309/149‑1 (MQ=255)
                                                                                                                              aTTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGtt                      >  5:50366/1‑149 (MQ=255)
                                                                                                                               ttGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTc                     >  2:117973/1‑149 (MQ=255)
                                                                                                                                   ccTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGAc                  >  6:383133/1‑148 (MQ=255)
                                                                                                                                   ccTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACt                 <  4:239171/149‑1 (MQ=255)
                                                                                                                                         tctGCCGCACTCCAGTGTGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCACTTCACACCGCCACTGGCGATCCCCTCGATCGACTATCAg            >  8:296248/3‑148 (MQ=255)
                                                                                                                                            gCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAg        >  4:342764/1‑149 (MQ=255)
                                                                                                                                                  cTGGTGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAGCCTGGc  <  4:55639/149‑1 (MQ=255)
                                                                                                                                                  cTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAGCCTGGc  <  6:50366/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGCCTGCCAGTGGCGATATGGCCCTGGAAGACTGGCTGGCACGCCCACATGTATTGATGGCCATGCGCCCGGATGCTTTCGATGAAATCGAGCGTGCACTGGCCGCCCGGGGCTTGCAGCGCCGTATTGCCCTGGCCTTGCCACACTGGAGTGCGGCAGTCGAAGTGCTGGCCGGCACCGACCTGATCCTGACAGTGGCCAGCCGCGCTGTGGGGTCGTTGCGCCAACACAAGACACTGCGCCAGTTCACACCGCCACTGGCCATCCCCTCGTTCGACTATCAGCAAGCCTGGC  >  NC_002947/4100685‑4100979

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