Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,317,956 C→T 6.3% A212V (GCG→GTG)  U712_06715 → Putative mannosyltransferase ykcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,317,9560CT 6.3% 209.7 / 2.1 63A212V (GCG→GTG) U712_06715Putative mannosyltransferase ykcB
Reads supporting (aligned to +/- strand):  ref base C (20/39);  new base T (2/2);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGCAGCATTTTTCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGTTTTGGAATTGGCTTTCGGCTATAACGGAATC  >  CP006881/1317819‑1318089
                                                                                                                                         |                                                                                                                                     
agcagcATTTTTCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTCCCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCg                                                                                                                                      >  6:67135/1‑139 (MQ=255)
  cagcaTTTTTCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGtt                                                                                                                                    <  3:132982/139‑1 (MQ=255)
  cagcaTTTTTCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGtt                                                                                                                                    <  3:127474/139‑1 (MQ=255)
   agcaTTTTTCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTa                                                                                                                                   <  6:126878/139‑1 (MQ=255)
        ttttCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGccct                                                                                                                              >  1:366810/1‑137 (MQ=255)
           tCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTa                                                                                                                                   <  2:318020/131‑1 (MQ=255)
           tCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTa                                                                                                                                   >  1:318020/1‑131 (MQ=255)
                gttgGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTagca                                                                                                                       <  2:134486/138‑1 (MQ=255)
                  tgGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGcc                                                                                                                    >  2:319887/1‑139 (MQ=255)
                    gAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCg                                                                                                                   <  5:34323/138‑1 (MQ=255)
                    gAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCg                                                                                                                   <  7:85086/138‑1 (MQ=255)
                    gAGTAGGTGTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCg                                                                                                                   <  6:276336/138‑1 (MQ=255)
                                 tACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAATTGTGTTAGCAGTATTAGCAGCAGTTTAGTTATAGttg                                                                                                     >  2:275992/1‑137 (MQ=255)
                                 tACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGTGTTATCCGTATTAGCAGCCTTTTCGTTATCGttt                                                                                                     >  2:281029/1‑137 (MQ=255)
                                      aaaTGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCt                                                                                                <  8:15420/139‑1 (MQ=255)
                                       aaTGCTGCAGGCTTACAAGAGATTCCCGGCTTTCTTATTCTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGGTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTc                                                                                               <  8:189469/139‑1 (MQ=255)
                                        aTGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              <  8:68319/139‑1 (MQ=255)
                                        aTGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              <  1:335096/139‑1 (MQ=255)
                                             gcAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCGAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCAttgtt                                                                                         <  7:287082/139‑1 (MQ=255)
                                               aGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGa                                                                                       <  7:185496/139‑1 (MQ=255)
                                                     aCATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAata                                                                                  >  3:206647/1‑138 (MQ=255)
                                                      cATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAatat                                                                                 <  7:97535/138‑1 (MQ=255)
                                                      cATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAatat                                                                                 <  6:288401/138‑1 (MQ=255)
                                                      cATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAatat                                                                                 <  2:225710/138‑1 (MQ=255)
                                                           ttttACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTc                                                                                               >  4:274293/1‑119 (MQ=255)
                                                           ttttACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTc                                                                                               <  3:274293/119‑1 (MQ=255)
                                                             ttACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTGATTGTTTACCatatacatga                                                                           >  2:224730/1‑131 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTGATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  7:133756/139‑1 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  3:82575/139‑1 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  5:67135/139‑1 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  6:205626/139‑1 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  4:117500/139‑1 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  2:42483/139‑1 (MQ=255)
                                                                  ggCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGcaagcaag                                                                    <  1:319887/139‑1 (MQ=255)
                                                                   gCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCaagcaagc                                                                   >  2:298686/1‑139 (MQ=255)
                                                                   gCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCaagcaagc                                                                   <  7:64256/139‑1 (MQ=255)
                                                                               ttCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAAACGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCg                                                       >  4:247135/1‑139 (MQ=255)
                                                                                 cTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              >  7:52800/1‑98 (MQ=255)
                                                                                 cTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              <  8:52800/98‑1 (MQ=255)
                                                                                  tATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTCGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCt                                                              >  6:267265/1‑129 (MQ=255)
                                                                                   aTCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              >  1:254958/1‑96 (MQ=255)
                                                                                   aTCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              <  2:254958/96‑1 (MQ=255)
                                                                                   aTCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAag                                                   >  4:125451/1‑139 (MQ=255)
                                                                                   aTCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAag                                                   >  1:258179/1‑139 (MQ=255)
                                                                                    tCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGa                                                     <  4:51326/136‑1 (MQ=255)
                                                                                    tCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAagc                                                  >  4:108281/1‑139 (MQ=255)
                                                                                        aTTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAagcagc                                               <  1:224730/138‑1 (MQ=255)
                                                                                         ttGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              <  1:274805/90‑1 (MQ=255)
                                                                                         ttGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCa                                                                                              >  2:274805/1‑90 (MQ=255)
                                                                                         ttGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGATGCAGcc                                              >  3:125495/1‑138 (MQ=255)
                                                                                                 aaCGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGTGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  1:281029/139‑1 (MQ=255)
                                                                                                 aaCGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  8:315687/139‑1 (MQ=255)
                                                                                                 aaCGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  1:298686/139‑1 (MQ=255)
                                                                                                 aaCGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  2:151969/139‑1 (MQ=255)
                                                                                                 aaCGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  7:247372/139‑1 (MQ=255)
                                                                                                 aaCGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  2:72938/139‑1 (MQ=255)
                                                                                                 aaCGCCAAAATCAAGAAAAAAATCGTTTCATTGGTAAGTGTGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAAct                                     <  1:275992/139‑1 (MQ=255)
                                                                                                       acaATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGttt                               >  7:290530/1‑139 (MQ=255)
                                                                                                       acaATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTAGATCGGAAGCAGCCAAACCAACTCTGttt                               >  2:4982/1‑139 (MQ=255)
                                                                                                              aGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGTTTTGGAAtt                        >  8:83091/1‑139 (MQ=255)
                                                                                                                     aaTCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGTTTTGGAATTGGCTTTc                 <  3:125451/139‑1 (MQ=255)
                                                                                                                              aTTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGTTTTGGAATTGGCTTTCGGCTATaa         <  4:169898/138‑1 (MQ=255)
                                                                                                                                     aGTGCGTTAGCCGTATTAGCAGCCGGTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGTTTTGGAATTGGCTTTCGGCTATAACGGAATc  <  2:233527/138‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                     
AGCAGCATTTTTCGTTGTTGGAGTAGGTTTTAATACAAAAATGCTGCAGGCTTACATGATTTTACCGGCTTTCTTATTGTTCTATCTGATTGCGGCAAACGCCACAATCAAGAAAAAAATCGTTTCATTGGTAAGTGCGTTAGCCGTATTAGCAGCCGTTTCGTTATCGTGGCCGCTCATTGTTGACAATATCCCTGCAAGCAAGCGGCCTTATATCGGAAGCAGCCAAACCAACTCTGTTTTGGAATTGGCTTTCGGCTATAACGGAATC  >  CP006881/1317819‑1318089

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: