Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 247,399 G→C 6.5% S526T (AGT→ACT)  yheS → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947247,3990GC 6.5% 191.3 / 3.4 61S526T (AGT→ACT) yheSABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base C (2/2);  total (35/26)
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

GCCGGTGCAGTGGTGGTGGTGTCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGCTGGAAACCGAGCTCAACCA  >  NC_002947/247251‑247534
                                                                                                                                                    |                                                                                                                                       
gCCGGTGCAGTGGTGGTGGTGTCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGAATTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGt                                                                                                                                        >  1:48439/1‑150 (MQ=255)
            gtggtggtgTCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGt                                                                                                                            >  1:192355/1‑150 (MQ=255)
                    gtCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCGGCGCACTGCGCTGGcggtac                                                                                                                           >  2:447583/1‑138 (MQ=255)
                    gtCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCGGCGCACTGCGCTGGcggtac                                                                                                                           <  1:447583/143‑6 (MQ=255)
                     tCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCg                                                                                                                   <  1:501397/150‑1 (MQ=255)
                       ccACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGc                                                                                                                 >  2:1279526/1‑150 (MQ=255)
                       ccACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGc                                                                                                                 >  2:870939/1‑150 (MQ=255)
                       ccACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGc                                                                                                                 <  1:1273397/150‑1 (MQ=255)
                       ccACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGc                                                                                                                 <  2:1005066/150‑1 (MQ=255)
                       ccACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGc                                                                                                                 >  1:809030/1‑150 (MQ=255)
                        cACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCa                                                                                                                              >  2:943867/1‑136 (MQ=255)
                        cACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCa                                                                                                                              <  1:943867/136‑1 (MQ=255)
                                   tttGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACaa                                                                                                     >  2:479368/1‑150 (MQ=255)
                                        tCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCTACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACAGCCAGCGCAGTGCGCCGGCCAGTACTTCACCGGCCAACCCGGACAAGACCg                                                                                                >  2:154558/1‑150 (MQ=255)
                                        tCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCg                                                                                                >  2:1210215/1‑150 (MQ=255)
                                          aaGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCg                                                                                                          <  2:1279206/138‑1 (MQ=255)
                                          aaGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCg                                                                                                          >  1:1279206/1‑138 (MQ=255)
                                              gCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGa                                                                                          >  2:536045/1‑150 (MQ=255)
                                               caccacTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTCCTGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGaa                                                                                         >  1:152403/1‑150 (MQ=255)
                                                  cacTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCGGCGCACTGCGCTGGcgg                                                                                                                              >  1:924591/1‑108 (MQ=255)
                                                  cacTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTg                                                                                                                                       <  1:790095/101‑1 (MQ=255)
                                                  cacTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTg                                                                                                                                       >  2:790095/1‑101 (MQ=255)
                                                  cacTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGc                                                                                      >  1:693040/1‑150 (MQ=255)
                                                  cacTCATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTCCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCGGCGCACTGCGCTGGcgg                                                                                                                              <  2:924591/110‑3 (MQ=255)
                                                   acTGATGACTTCCTGCTGGTGGCCGACGGCAAGGACTACACCTTCGACGGTGACCTGCAGGACTACACCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGcc                                                                                     >  1:930256/1‑150 (MQ=255)
                                                     tgatgaCTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCa                                                                                   <  2:693040/150‑1 (MQ=255)
                                                      gatgaCTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAg                                                                                  <  2:511069/150‑1 (MQ=255)
                                                       atgaCTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAgc                                                                                 >  1:502150/1‑150 (MQ=255)
                                                           cTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCa                                                                             >  1:717303/1‑150 (MQ=255)
                                                            ttCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAg                                                                            >  1:660471/1‑150 (MQ=255)
                                                            ttCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAg                                                                            <  2:809030/150‑1 (MQ=255)
                                                                tgctgGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTg                                                                                                                                       >  1:723098/1‑87 (MQ=255)
                                                                tgctgGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTg                                                                                                                                       <  2:723098/87‑1 (MQ=255)
                                                                tgctgGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAAccc                                                                                                           <  2:1068883/115‑1 (MQ=255)
                                                                tgctgGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAAccc                                                                                                           >  1:1068883/1‑115 (MQ=255)
                                                                        gCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGa                                                                                                        >  1:217233/1‑110 (MQ=255)
                                                                        gCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGa                                                                                                        <  2:217233/110‑1 (MQ=255)
                                                                          cGACGGCAAGGTCGCCACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACACGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCtt                                                              >  2:604542/1‑150 (MQ=255)
                                                                                   ggTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTcagca                                                     <  2:502150/150‑1 (MQ=255)
                                                                                     tCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGt                                                   <  1:479368/150‑1 (MQ=255)
                                                                                     ccaacacCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGt                                                   <  2:930256/147‑1 (MQ=255)
                                                                                            cTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCAcc                                            >  2:1285569/1‑150 (MQ=255)
                                                                                             ttCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCg                                           >  1:1129880/1‑150 (MQ=255)
                                                                                               cGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGAAGTTGGCACCGca                                         <  2:660471/150‑1 (MQ=255)
                                                                                                        ccTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGc                                <  2:1129880/150‑1 (MQ=255)
                                                                                                         cTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCt                               >  2:196348/1‑150 (MQ=255)
                                                                                                                aCTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGACAGCCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAGGGCTGCCGACa                        >  1:288754/1‑150 (MQ=255)
                                                                                                                  tACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGaaa                                                                                        <  2:247671/84‑1 (MQ=255)
                                                                                                                  tACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGaaa                                                                                        >  1:247671/1‑84 (MQ=255)
                                                                                                                   aCAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGc                     <  1:870939/150‑1 (MQ=255)
                                                                                                                     aGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCg                                                                                                                   <  2:538644/54‑1 (MQ=255)
                                                                                                                     aGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCg                                                                                                                   >  1:538644/1‑54 (MQ=255)
                                                                                                                     aGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGCTg                   <  1:536045/150‑1 (MQ=255)
                                                                                                                      tgcGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGCCACCGCAAGAGAAGGCGGGCGAGGAGCTgg                  >  2:73703/3‑150 (MQ=255)
                                                                                                                             gcttgTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGCTGGAAACCGa           >  2:638717/1‑150 (MQ=255)
                                                                                                                              ctggTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGCTGGAAACCGAg          <  1:196348/150‑1 (MQ=255)
                                                                                                                                     aGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTg                                                             >  2:1249505/1‑92 (MQ=255)
                                                                                                                                     aGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTg                                                             <  1:1249505/92‑1 (MQ=255)
                                                                                                                                      gTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGCTGGAAACCGAGCTCAACCa  >  2:1150178/1‑150 (MQ=255)
                                                                                                                                           gCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACaagaa                                   <  1:1170948/112‑1 (MQ=255)
                                                                                                                                           gCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACaagaa                                   >  2:1170948/1‑112 (MQ=255)
                                                                                                                                                    gTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGAc                                                                                                  <  2:1280140/40‑1 (MQ=255)
                                                                                                                                                    gTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGAc                                                                                                  >  1:1280140/1‑40 (MQ=255)
                                                                                                                                                    |                                                                                                                                       
GCCGGTGCAGTGGTGGTGGTGTCCCACGACCGTCATTTGCTCAAGAGCACCACTGATGACTTCCTGCTGGTGGCCGACGGCAAGGTCGACACCTTCGACGGTGACCTGGAGGACTACAGCCGCTGGCTGGTCGAGTACCGCCAGCGCAGTGCGCCGGCCAGTACTGCACCGGCCAACCCGGACAAGACCGACAAGAAAGCCCAGCGCCAGGCTGCCGCAGCCTTGCGTCAGCAGTTGGCACCGCACAAGAAGGCTGCCGACAAGCTGGAAACCGAGCTCAACCA  >  NC_002947/247251‑247534

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