Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,919,639 G→T 7.3% K75N (AAG→AAT PP_5190 → type II secretion system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,919,6390GT 7.3% 172.4 / 3.8 55K75N (AAG→AATPP_5190type II secretion system protein
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base T (2/2);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

ATACCCTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTT  >  NC_002947/5919497‑5919787
                                                                                                                                                   |                                                                                                                                                    
aTACCCTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGc                                                                                                                                                <  4:128341/149‑1 (MQ=255)
    ccTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAAc                                                                                                                                            >  3:193146/1‑149 (MQ=255)
    ccTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAAc                                                                                                                                            <  3:77743/149‑1 (MQ=255)
    ccTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAAc                                                                                                                                            <  4:139566/149‑1 (MQ=255)
     cTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAAcc                                                                                                                                           <  5:6053/149‑1 (MQ=255)
         aaGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTAc                                                                                                                                       <  7:256767/149‑1 (MQ=255)
              cTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTgcg                                                                                                                                  >  5:285397/1‑149 (MQ=255)
                     acGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGAc                                                                                                                           <  5:68588/149‑1 (MQ=255)
                         cgcTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCt                                                                                                                                     >  4:81212/1‑135 (MQ=255)
                         cgcTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCt                                                                                                                                     <  3:81212/135‑1 (MQ=255)
                               ccACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAAc                                                                                                                 >  2:206101/1‑149 (MQ=255)
                                cACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTAGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGTTGGCCGAC‑‑CTTGTCGGCCAGcc                                                                                                                                           >  1:218990/1‑123 (MQ=255)
                                cACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGTTGGCCGAC‑‑CTTGTCGGCCAGcc                                                                                                                                           <  2:218990/122‑1 (MQ=255)
                                cACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACt                                                                                                                >  3:49881/1‑149 (MQ=255)
                                 aCCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGc                                                                                                                    >  1:289406/1‑144 (MQ=255)
                                 aCCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGc                                                                                                                    <  2:289406/144‑1 (MQ=255)
                                 aCCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCATGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTc                                                                                                               >  8:88204/1‑149 (MQ=255)
                                    ccGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCATGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGAc                                                                                                            <  7:88204/149‑1 (MQ=255)
                                          gCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGc                                                                                                        <  6:285397/147‑1 (MQ=255)
                                          gCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGc                                                                                                        <  5:72630/147‑1 (MQ=255)
                                          gCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGc                                                                                                        <  4:193146/147‑1 (MQ=255)
                                            cACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCAACCCGGCCAGCCCCTGGACTTGGAGCAACTTTGCCAGTGGCTGGCCGAC‑‑AAGGTAGGCCAACCCTACCTGCGCATAGACCCCATGCAACTCGCCCTGccc                                                                                                      >  7:244654/1‑147 (MQ=255)
                                                           aTCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGAt                                                                                     >  1:174850/1‑149 (MQ=255)
                                                                cgccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCCTGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTccc                                                                                <  7:125810/149‑1 (MQ=255)
                                                                cgccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCa                                                                                            <  5:88791/137‑1 (MQ=255)
                                                                cgccTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCa                                                                                            >  6:88791/1‑137 (MQ=255)
                                                                  ccTGCAGCCTGGAAAGCCGCCAGCACCCCGGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTccc                                                                                >  1:22012/1‑148 (MQ=255)
                                                                      cAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGTC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCt                                                                                        <  8:42186/135‑1 (MQ=255)
                                                                      cAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCa                                                                              <  7:244510/145‑1 (MQ=255)
                                                                      cAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCCACCTGCCCCAGGTCACCAGCCTGATCTCCCCa                                                                              >  8:244510/1‑145 (MQ=255)
                                                                       aGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCATGTCACCAGcc                                                                                         <  3:64948/133‑1 (MQ=255)
                                                                       aGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGcc                                                                                         >  4:64948/1‑133 (MQ=255)
                                                                        gCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTc                                                                        <  4:49881/149‑1 (MQ=255)
                                                                         ccTGGAAAGCCGCCAGGACCCCGGCCAGCCCCTTGACTGTGAGAACCTTTGCCAGTGGCTGGCCGAC‑‑AATGTCGGCCAACCCTACCTGCGCATCGACCCCATTCAACTCGACCTGCCCCAGGTCACCATCCTGATCTCCCCAGCCTTCg                                                                       <  8:244654/149‑1 (MQ=255)
                                                                         ccTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCGCCAGGTCACCAGCCTGATCTCCCCAGCCTTCg                                                                       <  6:284952/149‑1 (MQ=255)
                                                                           tGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGcc                                                                     >  8:85805/1‑149 (MQ=255)
                                                                                           aCCCCGGCCAGCCCCTGGCCTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGAGCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCAGGGCATc                                                      <  6:156777/149‑1 (MQ=255)
                                                                                             cccGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATcc                                                     >  4:192174/1‑148 (MQ=255)
                                                                                             cccGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCt                                                    >  3:172628/1‑149 (MQ=255)
                                                                                             cccGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCt                                                    >  1:232052/1‑149 (MQ=255)
                                                                                              ccGGCCAGCCCCTGGTCTTGGAGCACCTTTGCCATTGGGTGGGCGACTTACAGTTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATc                                                      <  1:98770/149‑1 (MQ=255)
                                                                                               cGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCa                                                  >  7:172913/1‑149 (MQ=255)
                                                                                                       ccccTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTcgc                                           >  4:252608/1‑149 (MQ=255)
                                                                                                         ccTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCgg                                         >  3:206030/1‑149 (MQ=255)
                                                                                                          cTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGc                                        <  1:206101/149‑1 (MQ=255)
                                                                                                                   gAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTc                               <  7:85805/149‑1 (MQ=255)
                                                                                                                        ccTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGACCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCg                          >  1:21356/1‑149 (MQ=255)
                                                                                                                              ccAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCAcc                     <  2:232052/148‑1 (MQ=255)
                                                                                                                               cAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGt                   >  1:18335/1‑149 (MQ=255)
                                                                                                                                aGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTc                  <  2:22012/149‑1 (MQ=255)
                                                                                                                                    gctggcCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGAt                                    >  7:52821/1‑127 (MQ=255)
                                                                                                                                    gctggcCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGAt                                    <  8:52821/127‑1 (MQ=255)
                                                                                                                                     ttggcCGAC‑‑CTTGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAg             >  2:106826/2‑149 (MQ=255)
                                                                                                                                         ccGAC‑‑AAGGTCGGCCAACCCTACATGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGCCATCCTCATCGTCCCGGCCGAAGCCTCCAGCGTCACCGACGCCAGcccc         >  8:119675/1‑149 (MQ=255)
                                                                                                                                          cGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGccc        <  2:18335/149‑1 (MQ=255)
                                                                                                                                                   ggTCGGCCAACCCTACCTGCGCATCGCCCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCtt  >  8:67471/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATACCCTGCAAGTACTCGAACACGCCGCTGCCCACCCCGCTAGCCACCCGCTGGAGCACATCGCCGCCTGCAGCCTGGAAAGCCGCCAGCACCCCGGCCAGCCCCTGGACTTGGAGCACCTTTGCCAGTGGCTGGCCGAC‑‑AAGGTCGGCCAACCCTACCTGCGCATCGACCCCATGCAACTCGACCTGCCCCAGGTCACCAGCCTGATCTCCCCAGCCTTCGCCCAGCGCCACGGCATCCTCATCGTCGCGGCCGATGCCTCCAGCGTCACCGTCGCCAGCGCCCAGCCTT  >  NC_002947/5919497‑5919787

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: