Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 39,378 C→A 100% D586Y (GAT→TAT)  mecA ← penicillin binding protein 2 prime

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073039,3780CA98.2% 149.9 / ‑3.5 57D586Y (GAT→TAT) mecApenicillin binding protein 2 prime
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (29/26);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCAT  >  CP000730/39241‑39512
                                                                                                                                         |                                                                                                                                      
aaTAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAGATAtattc                                                                                                                                     >  1:61833/1‑141 (MQ=255)
aaTAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAtattc                                                                                                                                     >  1:658844/1‑141 (MQ=255)
 aTAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAtattct                                                                                                                                    >  2:463786/1‑141 (MQ=255)
 aTAGCCATCATCATGTTTGGATTATCTTTATAATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATAtattct                                                                                                                                    >  1:468989/1‑141 (MQ=255)
  tAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCTATTTGTCTTCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATattctt                                                                                                                                   <  1:738932/141‑1 (MQ=255)
     tcctcatcatGTTTGGATTATCTTTATCATATTATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTAt                                                                                                                                <  2:194700/138‑1 (MQ=255)
     ccATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTATGCAGTACCGGATTGGCCAATTAAGTTTGCATAAGATCTATAAATCTATTCtttct                                                                                                                                >  2:279858/1‑139 (MQ=255)
         catcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATgtgt                                                                                                                            >  2:235088/1‑141 (MQ=255)
         caccatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATgtgt                                                                                                                            <  2:389378/138‑1 (MQ=255)
          atcatGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGtt                                                                                                                           >  2:245718/1‑141 (MQ=255)
             atGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTttta                                                                                                                        >  2:116374/1‑141 (MQ=255)
             atGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTATTTATGTGTttta                                                                                                                        >  2:528122/1‑141 (MQ=255)
             atGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCATTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTATTTATGGGttttt                                                                                                                        >  2:685769/1‑140 (MQ=255)
             atGTTTGGATTATCTTTATAATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTttta                                                                                                                        >  2:637108/1‑141 (MQ=255)
                 ttGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTtttattta                                                                                                                    <  1:334048/141‑1 (MQ=255)
                    gATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTCTGTGTTTTATTTACGa                                                                                                                 >  2:711344/1‑141 (MQ=255)
                         tCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTATGACttat                                                                                                            >  2:617005/1‑139 (MQ=255)
                            ttATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGc                                                                                                         <  1:672900/141‑1 (MQ=255)
                             tATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCa                                                                                                        <  1:724320/141‑1 (MQ=255)
                             tATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCGCATTGGTTCAGTTTGAGTTCTGCAGTACCGGATTGGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCa                                                                                                        <  1:609326/141‑1 (MQ=255)
                              ttCATCTGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATATCTTATATTTACGACTTGTTGCAt                                                                                                       <  2:762182/140‑1 (MQ=255)
                              ttCATATGGTTGAAACCACCCAATTTGTCTGCCAGTTTCGCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCAt                                                                                                       <  2:47697/140‑1 (MQ=255)
                               tCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATa                                                                                                      >  2:47041/1‑141 (MQ=255)
                                    tGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTTCATACCAta                                                                                                 >  1:254044/1‑140 (MQ=255)
                                      atatAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTCTGTGGTGTATTTACGACTTGTTGCATACCATCAg                                                                                               <  2:328361/141‑1 (MQ=255)
                                        ataAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATCAGATCTCTAAGTAAATTATTTATGTCTTTTATTTACACATTGTTGCATAACGTCAAtt                                                                                             >  1:355607/1‑141 (MQ=255)
                                        ataAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGtt                                                                                             <  1:548520/141‑1 (MQ=255)
                                          aaaCCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTCGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTaa                                                                                           <  2:701463/141‑1 (MQ=255)
                                            accaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACTATCAGTTAATa                                                                                         >  1:498070/1‑141 (MQ=255)
                                            accaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATa                                                                                         >  1:767155/1‑141 (MQ=255)
                                             ccaccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAg                                                                                        >  2:609129/1‑141 (MQ=255)
                                              caccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAgg                                                                                       >  2:684241/1‑140 (MQ=255)
                                              caccCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGa                                                                                       >  2:93550/1‑141 (MQ=255)
                                                 ccAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCATTTAATAGATTg                                                                                    >  2:724491/1‑141 (MQ=255)
                                                 ccAATTTGTCTGCCAGTTCCGGCCGCTGTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCCTCAGTTAATAGATTg                                                                                    >  2:364470/1‑141 (MQ=255)
                                                  cAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCCGTTAATAGATTGa                                                                                   >  2:131825/1‑141 (MQ=255)
                                                  cAATTTGTCGTCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGAACCGTCGCATACCATCAGTTAATAGATTGa                                                                                   <  2:170117/141‑1 (MQ=255)
                                                   aaTTTGTCTGCCAGTTTCTCCTTGTTTAATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCCTCAGTTAATAGATTGat                                                                                  >  2:327788/1‑141 (MQ=255)
                                                    aTTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAACTCTTTTCCTGGTGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGata                                                                                 >  1:537807/1‑141 (MQ=255)
                                                      ttGTCTGCCAGTGTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAACTATATTCTTTGTGTGTTTTATTTAAGAATTTTTGGATCCAATCAGTTAAAAGCTTGattat                                                                               >  1:335890/1‑138 (MQ=255)
                                                        gTCTGTCATTTTCTTCTTGTTTAATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAtttt                                                                             <  1:51832/141‑1 (MQ=255)
                                                        gTCTGCCAGTTTCTTCTTGTTTCATTTTGAGTTCTGAAGTACCGGATTTGCCAACTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCGGTTAATAGATTGATAtttt                                                                             <  1:442019/141‑1 (MQ=255)
                                                                gTTTTTCCTTGTTTCATTTTGAGTTATGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATCTTCTGTGGaa                                                                     <  2:93423/141‑1 (MQ=255)
                                                                  tcctcCTTGGTGCATTGTGGGTTGCGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGCaat                                                                   <  2:463081/139‑1 (MQ=255)
                                                                   tctcCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAaata                                                                  <  2:268223/141‑1 (MQ=255)
                                                                                ttttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATCTTTCTTcc                                                     >  2:549935/1‑141 (MQ=255)
                                                                                 tttGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCa                                                    >  1:48693/1‑141 (MQ=255)
                                                                                   tGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCaaa                                                  >  2:81748/1‑141 (MQ=255)
                                                                                    gAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAAc                                                 >  2:606749/1‑141 (MQ=255)
                                                                                     aGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATACTTTTCTTCCAAACt                                                <  1:612877/141‑1 (MQ=255)
                                                                                       ttCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACttt                                              <  2:290390/141‑1 (MQ=255)
                                                                                        tCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCAGACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTg                                             <  1:713939/141‑1 (MQ=255)
                                                                                          tGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGtt                                           <  2:209726/141‑1 (MQ=255)
                                                                                                   ggATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTGCGACTTGTTGTATCCCATCAGTTAATATATTGATATTTTTTTTGGAAATTATATTTTTGTTTCAATCTTTGTTTTTCgtgtt                                  >  2:19555/1‑140 (MQ=255)
                                                                                                         gCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTATTTGGAAATAATATTTTTCTTCCACACTTTGTTTTTCGTCTCttttca                            >  2:178771/1‑139 (MQ=255)
                                                                                                               ttaGTTTTCATAAGATCTATAAATATATTCTTTCTTTGTTTTATTTATGACTTCTCGCATACCATTAGTTAATAGTTTGATATTTTCTTTGGAAATTATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTg                      <  2:434311/139‑1 (MQ=255)
                                                                                                                     tGCATAATATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATGTTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATACGTGAGGTGc                <  2:507130/141‑1 (MQ=255)
                                                                                                                          aaGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTaa           <  1:203841/141‑1 (MQ=255)
                                                                                                                             aTCTATAAATATATTCTTTATGTGTTTTATTTACGCCTTGTCGGATAGCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAatat        <  1:609129/141‑1 (MQ=255)
                                                                                                                             aTCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAatat        <  1:463786/141‑1 (MQ=255)
                                                                                                                                  taAATTTATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTATAATAAGTGAGGTGCGTTAATATTGCCa   <  1:509351/141‑1 (MQ=255)
                                                                                                                                  aaattataTTCTCTATGTGTTTTATTTACGACTTGTTGCATACCAGCAGGTAATAGATTGATATTTTCTTTGGAAATAATATTTTACTTTCCAACTTTGTTTTTCGTGTTCTTTAATACGTGAGGTGCGTTAATATTGCCa   <  2:479388/137‑1 (MQ=255)
                                                                                                                                   aaaTATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCat  >  1:703268/1‑141 (MQ=255)
                                                                                                                                         |                                                                                                                                      
AATAGCCATCATCATGTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCCAT  >  CP000730/39241‑39512

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