Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,061,893 G→T 10.0% V270L (GTG→TTG)  PP_1838 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,061,8930GT10.0% 148.0 / 5.4 56V270L (GTG→TTG) PP_1838membrane protein
Reads supporting (aligned to +/- strand):  ref base G (30/20);  new base T (4/2);  total (34/22)
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.94e-01

ATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGCCGGCC  >  NC_002947/2061747‑2062034
                                                                                                                                                   |                                                                                                                                             
aTGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGTGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAgtgt                                                                                                                                              <  1:49415/148‑1 (MQ=255)
aTGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGtgtg                                                                                                                                             >  2:104790/1‑149 (MQ=255)
  gCCTTTACCACGGACTCCAACGTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGCGGACCGTGCCAACATCTGGAAGCCAGTCATGCAGCAGGCCGTGGCCACGCAAAGGTTGCg                                                                                                                                           >  7:328039/1‑149 (MQ=255)
   ccTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGt                                                                                                                                          <  7:3592/149‑1 (MQ=255)
            cGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAActg                                                                                                                                 >  1:54087/1‑149 (MQ=255)
                   cAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCg                                                                                                                          >  3:251307/1‑149 (MQ=255)
                   cAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCg                                                                                                                          <  8:204305/149‑1 (MQ=255)
                    aaCTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGc                                                                                                                         >  6:85103/1‑149 (MQ=255)
                     aCTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGcc                                                                                                                        <  5:85103/149‑1 (MQ=255)
                     aCTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGcc                                                                                                                        >  8:343537/1‑149 (MQ=255)
                        tCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGCACGACCCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCACGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCa                                                                                                                      <  1:103851/149‑1 (MQ=255)
                           tCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGCCCGCCCACGac                                                                                                                  >  6:114630/1‑149 (MQ=255)
                           tCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGac                                                                                                                  <  1:227374/149‑1 (MQ=255)
                                 aGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTg                                                                                                            >  1:353466/1‑149 (MQ=255)
                                 aGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTg                                                                                                            <  4:251307/149‑1 (MQ=255)
                                  gCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTgg                                                                                                           >  7:68703/1‑149 (MQ=255)
                                       gCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGAc                                                                                                      <  7:321228/149‑1 (MQ=255)
                                       gCCTGAACCGCACCCTGACCCTGATCGACCCCGCCAAGAGCCCCTTTGGCGAGGACCGTGCCAACATCTGCAAGCCACTCCTGCACCACGACGTGGCCACGCAAACTGTGCGAGAGCAACTGCTGACCGCCCACGACACCATCGTCGa                                                                                                       >  5:383740/1‑148 (MQ=255)
                                        ccTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACg                                                                                                     >  7:252151/1‑149 (MQ=255)
                                        ccTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACg                                                                                                     >  7:247212/1‑149 (MQ=255)
                                                 cACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAg                                                                                            <  3:12226/149‑1 (MQ=255)
                                                 cACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGAGTGAGCAACTGCGGACCGCCCACGACACCCTGCTCGACGCCGACGAAg                                                                                            >  4:246831/1‑149 (MQ=255)
                                                    cctgaccctgaTCGACGCCGCCAAAAGCCGCATTGGCGAGGACCGTGACAACATCTCGAAGCCAGTCCTGCCGCAGGACGACGCCACGCAAAGAGTCCGTGAGAAACTGCTGACCGCCCACGACACCCTGGTCGAAGCCGACGAAGCCg                                                                                         >  2:112784/1‑149 (MQ=255)
                                                             gaTCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATccgcc                                                                                >  4:80324/1‑149 (MQ=255)
                                                                ctaCGCCGCCATGAGCCGCTTTGGCGAGGACCGTGCCAAGATCTTGAAGCCAGTCCTGGAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGCCCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCcgccgc                                                                              <  1:278381/146‑1 (MQ=255)
                                                                    gccgccAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACtt                                                                         >  4:166257/1‑149 (MQ=255)
                                                                          aaGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTgcc                                                                   >  8:191694/1‑149 (MQ=255)
                                                                                cGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCgg                                                              <  8:68703/148‑1 (MQ=255)
                                                                                cGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGc                                                             >  5:174231/1‑149 (MQ=255)
                                                                                     tGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGAc                                                         <  5:165689/148‑1 (MQ=255)
                                                                                        cGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTa                                                     >  5:272262/1‑149 (MQ=255)
                                                                                              ccGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGc                                               <  3:80324/149‑1 (MQ=255)
                                                                                                   ccAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGCCCGTCCCCTGCCGAt                                           >  4:209028/1‑148 (MQ=255)
                                                                                                   ccAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATc                                          <  2:54087/149‑1 (MQ=255)
                                                                                                      aCATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATc                                          >  4:174649/1‑146 (MQ=255)
                                                                                                           tGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgt                                  >  6:74036/1‑149 (MQ=255)
                                                                                                           tGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgt                                  >  8:105940/1‑149 (MQ=255)
                                                                                                                  cAGTCCTGCTGTAGGACGTGGCGACGCTTAGTGTGCGTGAGCAAGTGCTGACCGCCTCCGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACgtggtggtg                           <  2:353466/149‑1 (MQ=255)
                                                                                                                        cTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGTTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATcc                       <  1:112784/148‑1 (MQ=255)
                                                                                                                        cTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCt                      >  3:290448/1‑149 (MQ=255)
                                                                                                                             gcagGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATg                  >  2:324053/1‑148 (MQ=255)
                                                                                                                               agGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGGCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGaa               <  3:246831/149‑1 (MQ=255)
                                                                                                                                 gACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGTTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAg             <  3:40758/149‑1 (MQ=255)
                                                                                                                                 gACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAg             >  2:231587/1‑149 (MQ=255)
                                                                                                                                     tGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTc         <  7:191694/149‑1 (MQ=255)
                                                                                                                                      ggCCACGCAAAGTGTGTGTGAGCAACTGCTGACCGCCCACGACACCCTGGTAGACGCAGACGAAGCAGCTATCCGCCGCGACTTTCTGCCGCCGGCTGACAGTACCCTGCCGATCAAGAACGTGGTGTTGATCCTGATGGAACGCTTCg        >  2:289688/1‑149 (MQ=255)
                                                                                                                                      gctcACGCACACTTTGCGTGAGCAACTGCTGACCACCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCg        >  8:2418/4‑149 (MQ=255)
                                                                                                                                       gCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGa                    <  1:247337/136‑1 (MQ=255)
                                                                                                                                       gCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACACTGCCGATCAAGAACGTGGTGGTGATCCTGa                    >  2:247337/1‑136 (MQ=255)
                                                                                                                                       gCCACGCAAAGTGTGCGAGAGCAACTGCTGACCGCCCACGCCACCCTGGTAGACGCCGACGAAGCAGCTATCCTCCGCCACTATGTGCAGCAGGCTGAGCGTAACATGACGATCAAGAACGTCGAGGAGATCCTGATGGACAGCTTccc       >  4:40559/1‑147 (MQ=255)
                                                                                                                                        tcACGCACACTTTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGc                                               >  5:14996/2‑108 (MQ=255)
                                                                                                                                        tcACGCACACTTTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGc                                               <  6:14996/107‑1 (MQ=255)
                                                                                                                                        tcACGCACACTTTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGcc      >  7:207806/2‑149 (MQ=255)
                                                                                                                                        ccACGCAAAGTTGGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGcc      <  4:339132/149‑1 (MQ=255)
                                                                                                                                        ccACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGcc      >  3:250917/1‑149 (MQ=255)
                                                                                                                                         cACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGCCg     >  6:30705/1‑149 (MQ=255)
                                                                                                                                            gCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGccaacc  >  2:351348/1‑145 (MQ=255)
                                                                                                                                            gCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGCCGGcc  <  6:50360/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ATGCCTTTACCACGGACTCCAACTTCGTCAACCAGCTGGGCCTGAACGGCACCCTGACCCTGATCGACGCCGCCAAGAGCCGCTTTGGCGAGGACCGTGCCAACATCTGGAAGCCAGTCCTGCAGCAGGACGTGGCCACGCAAAGTGTGCGTGAGCAACTGCTGACCGCCCACGACACCCTGGTCGACGCCGACGAAGCCGCTATCCGCCGCGACTTTGTGCCGCCGGCTGACCGTACCCTGCCGATCAAGAACGTGGTGGTGATCCTGATGGAAAGCTTCGCCGGCC  >  NC_002947/2061747‑2062034

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