Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,882,813 A→G 10.5% intergenic (+56/‑169) oprN → / → PP_3428 multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,882,8130AG10.5% 142.0 / 8.0 57intergenic (+56/‑169)oprN/PP_3428multidrug RND transporter outer membrane protein OprN/TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (2/4);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

AGGGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGC  >  NC_002947/3882666‑3882956
                                                                                                                                                   |                                                                                                                                               
aGGGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAg                                                                                                                                                >  2:15273/1‑149 (MQ=255)
    tGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGcc                                                                                                                                             >  6:6455/1‑148 (MQ=255)
      aaGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCAt                                                                                                                                           <  8:236141/148‑1 (MQ=255)
           gACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCcg                                                                                                                                     >  1:419639/1‑149 (MQ=255)
                 tATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCgggg                                                                                                                                <  5:384830/148‑1 (MQ=255)
                 tATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCggg                                                                                                                                 <  5:6455/147‑1 (MQ=255)
                 tATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCgggg                                                                                                                                <  7:94643/148‑1 (MQ=255)
                  aTCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCGGGCCATGGCCTGCGGGGc                                                                                                                               <  5:89465/148‑1 (MQ=255)
                  aTCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGc                                                                                                                               <  4:400771/148‑1 (MQ=255)
                  aTCGAGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCt                                                                                                                              >  2:292068/1‑149 (MQ=255)
                     gcgGCCTCGTCGCGCTCTACCAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCCCTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCtttt                                                                                                                           <  6:423893/149‑1 (MQ=255)
                      cgGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCGGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCttttt                                                                                                                          <  7:245022/149‑1 (MQ=255)
                          aTCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGTCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATg                                                                                                                      >  5:228264/1‑149 (MQ=255)
                          aTCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGAGCACTGATTCCACTGCCTCAACGACTCCATTGGTTGGCTCCTCCCCCAACCGTTTTCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATg                                                                                                                      >  5:234158/1‑149 (MQ=255)
                           tagTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCAAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGt                                                                                                                     <  2:110290/147‑1 (MQ=255)
                             gTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGttt                                                                                                                   <  7:51360/149‑1 (MQ=255)
                              tctcGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTg                                                                                                                  >  6:424883/1‑149 (MQ=255)
                                     cTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGACTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATct                                                                                                            >  6:303303/1‑147 (MQ=255)
                                      tACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTa                                                                                                          >  4:128248/1‑149 (MQ=255)
                                        cAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTAt                                                                                                         >  6:8146/1‑148 (MQ=255)
                                           ggCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAc                                                                                                     <  3:26267/149‑1 (MQ=255)
                                            gCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAcc                                                                                                    <  2:419639/149‑1 (MQ=255)
                                               cTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTTGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAcccc                                                                                                  >  5:306381/1‑148 (MQ=255)
                                                   gTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATAc                                                                                             >  3:404975/1‑149 (MQ=255)
                                                        ggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTTGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGt                                                                                        <  6:306381/149‑1 (MQ=255)
                                                        ggcTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGt                                                                                        >  2:264933/1‑149 (MQ=255)
                                                                         ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGcc                                                                                                                <  2:277134/108‑1 (MQ=255)
                                                                         ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGcc                                                                                                                >  1:277134/1‑108 (MQ=255)
                                                                         ggTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCg                                                                        >  5:288461/1‑148 (MQ=255)
                                                                           tGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAg                                                                     <  1:264933/149‑1 (MQ=255)
                                                                                 gcgcgcACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctgg                                                               >  1:241806/1‑149 (MQ=255)
                                                                                  cgcgcACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggc                                                              <  1:292068/149‑1 (MQ=255)
                                                                                        cTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctgtcgcctgg                                                        <  5:303303/149‑1 (MQ=255)
                                                                                           aTTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGgcctggcgcctggc                                                       <  5:424883/147‑1 (MQ=255)
                                                                                             tCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAAc                                                   >  1:164487/1‑149 (MQ=255)
                                                                                                aCTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGc                                                >  4:363036/1‑149 (MQ=255)
                                                                                                aCTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGGTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGTGCCTGGCCAACGGc                                                <  1:15273/149‑1 (MQ=255)
                                                                                                      tCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGc                                                >  2:42404/1‑143 (MQ=255)
                                                                                                      tCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGc                                                <  1:42404/143‑1 (MQ=255)
                                                                                                         aCGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATCTATGCCATCGTATACACCCCAATAAGCAGTGCATGGAATGGCACCGCAGGCCTGGCGCCTGGCCAACGGCGTCCGcgga                                       >  6:30183/1‑148 (MQ=255)
                                                                                                           gtcTCCTTTGGATGGGTCCTCGCCCATCCGTTTGCCGCGCAGGCCATGGCCCGCGGGGCTTTATTATGTTTGCGATCGTATTCACGCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttc                                     <  3:409710/147‑1 (MQ=255)
                                                                                                              tCCTTTGGTTGGCTCCTCCCCCACCCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCAACGTATTCCCCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCCACGGCGGGCGcggttcggt                                  >  6:128770/1‑149 (MQ=255)
                                                                                                                cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcggt                                  <  3:128248/147‑1 (MQ=255)
                                                                                                                cTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCGTTTTTATGTTTGCCATCGTATTGCACCCAAAACGCAGTGCATGGACTGCCACCGCAGGCCTGGCGCCTGGCCACCGGCGTCCCCGGTTCGGGTg                                >  4:271823/1‑149 (MQ=255)
                                                                                                                     gttggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCAGCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCAttt                           <  8:122566/149‑1 (MQ=255)
                                                                                                                        ggCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTTCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCCCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTAcc                        >  8:275372/1‑149 (MQ=255)
                                                                                                                           tcctccCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCTCGGTTCGGTTGCATTTAccccc                     >  2:196888/1‑149 (MQ=255)
                                                                                                                                 cccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTg               >  1:252789/1‑149 (MQ=255)
                                                                                                                                  ccAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTgg              <  2:241806/149‑1 (MQ=255)
                                                                                                                                    aaCCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGcc            <  4:404975/149‑1 (MQ=255)
                                                                                                                                      ccGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGccc           >  8:113932/1‑148 (MQ=255)
                                                                                                                                            gCCCCGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCcggttcgg                                   <  6:170978/118‑1 (MQ=255)
                                                                                                                                            gCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGctct    >  2:332175/1‑149 (MQ=255)
                                                                                                                                             ccccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTTGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTg   <  1:332175/149‑1 (MQ=255)
                                                                                                                                              cccGCGGGCCATGGCCTGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGc  >  7:342211/1‑149 (MQ=255)
                                                                                                                                              cccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGc  <  6:288461/149‑1 (MQ=255)
                                                                                                                                              cccGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGc  >  5:186996/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGGGTGAAGTCGACCTTTATCGCGGCATCGTCTCGATCTACAAGGCCCTGGGTGGCGGCTGGCAACCGGAGACGGTGGCCAGCGCGCACTGATTCCACTGCCTCAACGACTCCTTTGGTTGGCTCCTCCCCCAACCGTTTGCCCCGCAGGCCATGGCCCGCGGGGCTTTTTTATGTTTGCCATCGTATTCACCCCAATACGCAGTGCATGGAATGCCACCGCAGGCCTGGCGCCTGGCCAACGGCGTGCCCGGTTCGGTTGCATTTACCCCCCGGCTGGCCCAAGCTCTGC  >  NC_002947/3882666‑3882956

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: