Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,542,895 G→A 6.4% T225T (ACG→ACA PP_1353 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,542,8950GA 6.4% 182.6 / 2.7 62T225T (ACG→ACAPP_1353membrane protein
Reads supporting (aligned to +/- strand):  ref base G (31/27);  new base A (2/2);  total (33/29)
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.80e-01

ACCGCCAGCCGACGCGCAAGGTGGTGTTCGATGTGGGCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTC  >  NC_002947/1542753‑1543040
                                                                                                                                              |                                                                                                                                                 
agcGCCAGCCGTCGCGCAAGGTGGTGTTCGATGTGGGCGTTGACTACGAGGCCGACGTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTTTTCTGCAAGACCCGTCACCGCAGGCGGTCGTGTCGACGCTGGGc                                                                                                                                             <  3:251219/147‑1 (MQ=255)
acCGCCAGCCGACGCGCAAGGTGGTGTTCGATGTGGGCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTgg                                                                                                                                               >  1:275552/1‑147 (MQ=255)
acCGCCAGCCGACGCGCAAGGTGGTGTTCGATGTGGGCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTgg                                                                                                                                               <  2:275552/147‑1 (MQ=255)
                         gtTCGATGTGGGCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTgc                                                                                                                    >  4:319041/1‑149 (MQ=255)
                              aTGTGGGCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTtg                                                                                                               >  5:25808/1‑149 (MQ=255)
                                    gCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGAcc                                                                                                         <  8:196380/149‑1 (MQ=255)
                                       ttGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAg                                                                                                      <  3:319041/149‑1 (MQ=255)
                                           cTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCt                                                                                                  >  7:311267/1‑149 (MQ=255)
                                               gAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGa                                                                                              >  2:67049/1‑149 (MQ=255)
                                                   ccGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTAc                                                                                          <  1:67049/149‑1 (MQ=255)
                                                     gACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTtgtg                                                                                                             <  4:264961/128‑1 (MQ=255)
                                                     gACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTtgtg                                                                                                             >  3:264961/1‑128 (MQ=255)
                                                           cAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGAc                                                                                    <  8:311267/147‑1 (MQ=255)
                                                            aGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGt                                                                                  >  1:64917/1‑148 (MQ=255)
                                                                  cccGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGtt                                                                            >  6:361158/1‑148 (MQ=255)
                                                                   ccGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCa                                                                          >  5:27023/1‑149 (MQ=255)
                                                                                 tGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGc                                                                                                                                                  <  8:211300/63‑1 (MQ=255)
                                                                                         gCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGa                                                                                              <  2:327478/107‑1 (MQ=255)
                                                                                         gCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGa                                                    >  4:282318/1‑149 (MQ=255)
                                                                                         gCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGa                                                    >  4:208806/1‑149 (MQ=255)
                                                                                            aaGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCg                                                 <  2:248759/149‑1 (MQ=255)
                                                                                                 cccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTg                                             >  3:128130/1‑148 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTtgtg                                                                                                             <  8:37405/83‑1 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTtgtg                                                                                                             >  7:37405/1‑83 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACt                                                                                         <  3:94669/103‑1 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACt                                                                                         >  4:94669/1‑103 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTg                                                                                 <  1:65981/111‑1 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTg                                                                                 >  2:65981/1‑111 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGa                                                                <  8:217793/128‑1 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGa                                                                >  7:217793/1‑128 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTg                                             >  8:171512/1‑147 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTg                                             >  4:314873/1‑147 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGaa                                           <  3:314873/149‑1 (MQ=255)
                                                                                                  ccGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGaa                                           <  6:27023/149‑1 (MQ=255)
                                                                                                   cgcgTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAg                                          <  4:128130/149‑1 (MQ=255)
                                                                                                    gcgTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAgg                                         >  5:181135/1‑149 (MQ=255)
                                                                                                         ttCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCg                                      >  4:254364/1‑147 (MQ=255)
                                                                                                         ttCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGaa                                    <  2:64917/149‑1 (MQ=255)
                                                                                                           cTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAg                                          >  6:72400/1‑141 (MQ=255)
                                                                                                           cTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAg                                          <  5:72400/141‑1 (MQ=255)
                                                                                                           cTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAg                                   <  7:171512/148‑1 (MQ=255)
                                                                                                              cAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGAt                                                                               <  2:267090/101‑1 (MQ=255)
                                                                                                              cAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGCACGTGATGTTCATGCTCAACGACTATGCGCGTCACCGGTTGAAGGCCGAAGGGAt                               >  2:327995/1‑149 (MQ=255)
                                                                                                              cAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACAGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGAt                                                                               >  1:267090/1‑101 (MQ=255)
                                                                                                              cAAGACCCGGCACCGCAGGCGGTCGTGTCGACACTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGAt                               >  8:297879/1‑149 (MQ=255)
                                                                                                              cAAGACCCGGCACCGCAGGCGGTCGTGTCGACACTGGGCGACAGCTCGATCAACGTATCGCTGCGCTTGTGGACCTAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTTTTCTCGTGACCGGTTGATGGCCGAAGGGa                                >  8:293193/1‑148 (MQ=255)
                                                                                                                 gACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCg                             >  3:336163/1‑148 (MQ=255)
                                                                                                                   cccGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACa                          >  7:257852/1‑149 (MQ=255)
                                                                                                                    ccGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAg                                          >  3:156608/1‑132 (MQ=255)
                                                                                                                    ccGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAg                                          <  4:156608/132‑1 (MQ=255)
                                                                                                                           cGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCCTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGttt                  <  5:361158/149‑1 (MQ=255)
                                                                                                                           cGCAGGCGGTCGTGTCGACACTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGttt                  <  7:297879/149‑1 (MQ=255)
                                                                                                                           cGCAGGCGGTCGTGTCGACACTGGGCGACAGCTCGATCACCGTATCGATTCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGttt                  <  7:293193/149‑1 (MQ=255)
                                                                                                                            gCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATg                                                           <  1:12722/107‑1 (MQ=255)
                                                                                                                            gCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATg                                                           >  2:12722/1‑107 (MQ=255)
                                                                                                                            gCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTc                 <  3:208806/149‑1 (MQ=255)
                                                                                                                                gCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGt                                                              <  4:190472/100‑1 (MQ=255)
                                                                                                                                gCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGt                                                              >  3:190472/1‑100 (MQ=255)
                                                                                                                                    tCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCgg         >  6:136414/1‑149 (MQ=255)
                                                                                                                                        gtCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAACGGATCGACATTCCGTTTCCCCAGCGGGTTa     >  7:257168/1‑149 (MQ=255)
                                                                                                                                         tCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTAt    >  6:95106/1‑149 (MQ=255)
                                                                                                                                           gACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTc  <  8:257168/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
ACCGCCAGCCGACGCGCAAGGTGGTGTTCGATGTGGGCGTTGACTACGAGGCCGACCTGCAGAAGGCCCGCAATGTGTTGCTGGAACTGGCAAAGGACCCGCGTGTTCTGCAAGACCCGGCACCGCAGGCGGTCGTGTCGACGCTGGGCGACAGCTCGATCACCGTATCGCTGCGCTTGTGGACCAAGACCTCCGATTACTGGGACGTGATGTTCATGCTCAACGAGTATGCGCGTGACCGGTTGAAGGCCGAAGGGATCGACATTCCGTTTCCGCAGCGGGTTATTC  >  NC_002947/1542753‑1543040

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: