Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,002,200 T→G 90.9% M173R (ATG→AGG)  mecA → Adapter protein MecA 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,002,2000TG90.9% 192.7 / 5.2 62M173R (ATG→AGG) mecAAdapter protein MecA 1
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (33/23);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

AAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTAA  >  USA300TCH1516_ALE/1002062‑1002338
                                                                                                                                           |                                                                                                                                          
aaGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGACATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAggg                                                                                                                                          <  1:352112/141‑1 (MQ=255)
               aaaaGTCTTCAGGACGTAAAAATACAAGAATAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttatt                                                                                                                           <  1:503595/141‑1 (MQ=255)
               aaaaGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATACCAATCCAGTAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttatt                                                                                                                           >  1:451158/1‑141 (MQ=255)
                  cgacTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttattatt                                                                                                                        >  2:374900/4‑141 (MQ=255)
                  aGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttattatt                                                                                                                        >  1:386552/1‑141 (MQ=255)
                  aGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACttattatt                                                                                                                        >  1:398803/1‑141 (MQ=255)
                  aGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTGATATAGGGTTGATGGTACttattatt                                                                                                                        >  1:386356/1‑141 (MQ=255)
                         aacACGTAAAGAGACAAGAACAATGATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTTTTATATAGGGTTGATGGTACTTATTATTATGCTGt                                                                                                                 <  1:158824/139‑1 (MQ=255)
                          gCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGta                                                                                                                >  1:532813/1‑141 (MQ=255)
                             cGTAAAAATACAAGAACAAACATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTACGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATAGAGGGTTGATGGTACTTATTATTATGCTGTatat                                                                                                             >  2:341518/1‑141 (MQ=255)
                                aaaaaTACAAGAACACTCCTTGTGCAAGTTAACGTTTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACCTATTATTATGCTGTATAtttt                                                                                                          <  1:317071/141‑1 (MQ=255)
                                aaaaaTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAGTTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATAgttt                                                                                                          <  1:439492/141‑1 (MQ=255)
                                    aTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATa                                                                                                      >  2:229186/1‑141 (MQ=255)
                                    aTACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATa                                                                                                      >  2:532888/1‑141 (MQ=255)
                                          aaaCAATCAGAGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATg                                                                                                <  2:198255/140‑1 (MQ=255)
                                           aaCAATCATTGTGAAATTTAACGTTTTAGAAGTTGTTATTAATTATGCCTATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGt                                                                                               <  1:264/141‑1 (MQ=255)
                                           aaCAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTGTTATTATGCTGTATATTTTGATAGTCATGt                                                                                               >  1:149686/1‑141 (MQ=255)
                                           aaCAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGTTGTATATTTTGATAGTCATGt                                                                                               >  1:149719/1‑141 (MQ=255)
                                               aTCATTGTGAAATTTAACGATTTCGAAGAGGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGAt                                                                                           <  2:552916/141‑1 (MQ=255)
                                                 cATTGTGAAATTTAACGATTTAGACAGCTATTTGAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCa                                                                                         >  1:296528/1‑141 (MQ=255)
                                                 cATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAAGCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCa                                                                                         >  1:394796/1‑141 (MQ=255)
                                                  aTTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaa                                                                                        <  1:122410/141‑1 (MQ=255)
                                                  aTTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCaa                                                                                        <  1:122639/141‑1 (MQ=255)
                                                   ttGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTAATCaag                                                                                       >  2:372197/1‑141 (MQ=255)
                                                        aaaTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTc                                                                                  <  2:532813/141‑1 (MQ=255)
                                                         aaTTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCa                                                                                 <  1:269769/141‑1 (MQ=255)
                                                           tttAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCAtt                                                                               <  2:537609/141‑1 (MQ=255)
                                                               aCGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATg                                                                           >  1:220730/1‑141 (MQ=255)
                                                                    ttAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGt                                                                      >  1:431670/1‑141 (MQ=255)
                                                                        aagaTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACa                                                                  <  1:349909/141‑1 (MQ=255)
                                                                            aTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATAGTTTGATAGTCATGTTTATCAAGAAGTCATTAATGATAGTTACAGTc                                                               <  1:294304/141‑1 (MQ=255)
                                                                            aTGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGGTGGTTAGAGTc                                                               >  1:154441/1‑141 (MQ=255)
                                                                               ttattaATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAAt                                                            <  1:235971/141‑1 (MQ=255)
                                                                                   taattaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCt                                                        >  1:279474/1‑141 (MQ=255)
                                                                                      ttaTGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGa                                                     >  1:193753/1‑141 (MQ=255)
                                                                                        aatcATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACGGTCAATTGCTTGAAt                                                   <  1:308071/138‑1 (MQ=255)
                                                                                                cATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTAt                                           <  1:514652/141‑1 (MQ=255)
                                                                                                cATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGAGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTAt                                           >  1:335539/1‑141 (MQ=255)
                                                                                                 aTAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATc                                          <  2:225794/141‑1 (MQ=255)
                                                                                                    gCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCaa                                       <  1:479596/141‑1 (MQ=255)
                                                                                                       aTCCAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGATTGAATTTGCTTATCAAacag                                    >  2:266643/1‑141 (MQ=255)
                                                                                                          cAATAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAacagaca                                 <  1:68276/141‑1 (MQ=255)
                                                                                                             tAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGGATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaa                              <  2:19438/141‑1 (MQ=255)
                                                                                                             tAACTACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGacagaa                              >  1:546449/1‑141 (MQ=255)
                                                                                                                 tACAGAGTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGaaagaacaga                          >  2:224674/1‑141 (MQ=255)
                                                                                                                    agagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGt                       <  1:229186/141‑1 (MQ=255)
                                                                                                                    agagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGt                       <  1:532888/141‑1 (MQ=255)
                                                                                                                     gagTTTGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtt                      <  1:59459/141‑1 (MQ=255)
                                                                                                                       gTTTGAAGATTTGTTAGATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGttta                    >  2:325143/1‑141 (MQ=255)
                                                                                                                        tttGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtctat                   <  1:439904/141‑1 (MQ=255)
                                                                                                                        tttGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGAGAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttat                   <  2:487328/141‑1 (MQ=255)
                                                                                                                         ttGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttatt                  >  2:192073/1‑141 (MQ=255)
                                                                                                                          tGAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttattt                 >  1:280034/1‑141 (MQ=255)
                                                                                                                           gAAGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGtttattta                >  1:496481/1‑141 (MQ=255)
                                                                                                                            aaGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaa               >  2:71857/1‑141 (MQ=255)
                                                                                                                            aaGATTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaa               >  1:373589/1‑141 (MQ=255)
                                                                                                                            aaGATTTGTTATATAGGGTTGATGGTACGTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTaa               >  2:357499/1‑141 (MQ=255)
                                                                                                                               aTTTGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCGTGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTCTTTaaacg            <  2:451158/141‑3 (MQ=255)
                                                                                                                                 ttGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGAc          >  1:524434/1‑141 (MQ=255)
                                                                                                                                 ttGTTATATAGGGTTGATGGTACTTATTATTATGCTCTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACATAACAGAAGTTTATTTAAATGAc          >  1:423249/1‑141 (MQ=255)
                                                                                                                                  tGTTATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACt         >  2:389732/1‑141 (MQ=255)
                                                                                                                                    ttATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATgtgtat       >  2:349874/1‑136 (MQ=255)
                                                                                                                                    ttATATAGGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTAt       >  2:186289/1‑141 (MQ=255)
                                                                                                                                       tatagGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCt    <  1:303864/141‑1 (MQ=255)
                                                                                                                                         tatGGTTGATGGTACTTATGATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTaa  <  1:9227/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                          
AAGGTTCTAAACGTCAAAAGTCTTCAGCACGTAAAAATACAAGAACAATCATTGTGAAATTTAACGATTTAGAAGATGTTATTAATTATGCATATCATAGCAATCCAATAACTACAGAGTTTGAAGATTTGTTATATATGGTTGATGGTACTTATTATTATGCTGTATATTTTGATAGTCATGTTGATCAAGAAGTCATTAATGATAGTTACAGTCAATTGCTTGAATTTGCTTATCCAACAGACAGAACAGAAGTTTATTTAAATGACTATGCTAA  >  USA300TCH1516_ALE/1002062‑1002338

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