Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,416 A→G 6.9% N122S (AAC→AGC)  PP_3450 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,908,4160AG 6.9% 162.0 / 2.1 58N122S (AAC→AGC) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base G (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAG  >  NC_002947/3908276‑3908563
                                                                                                                                             |                                                                                                                                                   
cacaGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGGCGCGGTTGGCCATGGCCAGCcgcgtct                                                                                                                                             <  8:305254/149‑4 (MQ=255)
        aGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTc                                                                                                                                     <  3:362619/149‑1 (MQ=255)
                   cAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccac                                                                                                                            >  2:46446/1‑147 (MQ=255)
                     ggACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACcag                                                                                                                         >  1:410704/1‑148 (MQ=255)
                        cTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCAtt                                                                                                                                           >  2:224645/1‑127 (MQ=255)
                        cTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCAtt                                                                                                                                           <  1:224645/127‑1 (MQ=255)
                           gCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCGACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgcg                                                                                                                                                >  4:1013/1‑119 (MQ=255)
                           gCAACGCCTGGCCAATCCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGTTGGCCATGGCCAGCcgcg                                                                                                                                                <  3:1013/119‑1 (MQ=255)
                             aaCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCACTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGccc                                                                                                                  >  1:84660/1‑148 (MQ=255)
                              aCGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAAc                                                                                                               <  3:81385/149‑1 (MQ=255)
                                     gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggct                                                                                                        <  2:410704/149‑1 (MQ=255)
                                     gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                                          <  1:46446/147‑1 (MQ=255)
                                     gCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                                          <  6:332029/147‑1 (MQ=255)
                                       cAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggc                                                                                                         <  5:244750/146‑1 (MQ=255)
                                       cAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggc                                                                                                         >  6:244750/1‑146 (MQ=255)
                                        aaTGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctggctg                                                                                                       >  5:318683/1‑147 (MQ=255)
                                           gCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTc                                                                                                  >  2:360787/1‑149 (MQ=255)
                                            ccTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTAGCCTTGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTc                                                                                                  <  2:371720/148‑1 (MQ=255)
                                              tataTCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGc                                                                                                >  5:387790/1‑148 (MQ=255)
                                                   ccAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGcc                                                                                                                   >  6:10883/1‑124 (MQ=255)
                                                   ccAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGcc                                                                                                                   <  5:10883/124‑1 (MQ=255)
                                                    cAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGcc                                                                                                                                 <  8:326093/109‑1 (MQ=255)
                                                    cAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGcc                                                                                                                                 >  7:326093/1‑109 (MQ=255)
                                                    cAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGcc                                                                                         <  6:387790/149‑1 (MQ=255)
                                                      gCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTg                                                                                       <  1:360787/149‑1 (MQ=255)
                                                              acgaAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTg                                                                               >  8:182/1‑149 (MQ=255)
                                                                gaAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGCTCAACTgcc                                                                             >  3:81100/1‑148 (MQ=255)
                                                                        ccTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgca                                                                             >  6:168883/1‑141 (MQ=255)
                                                                        ccTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgca                                                                             <  5:168883/141‑1 (MQ=255)
                                                                             aCTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                                  <  5:171901/148‑1 (MQ=255)
                                                                                      ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGc                                                         <  6:318683/148‑1 (MQ=255)
                                                                                         gCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACCCTCGCTGCAGTCGCTACCGCCAGATCAACTGCAGCAACCCGAACAGCGCCAGCAACt                                                     >  2:193679/1‑149 (MQ=255)
                                                                                           ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                                   >  6:336811/1‑149 (MQ=255)
                                                                                           ccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                                   <  2:72469/149‑1 (MQ=255)
                                                                                            cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                                  >  4:317215/1‑149 (MQ=255)
                                                                                            cACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                                  <  2:84660/149‑1 (MQ=255)
                                                                                             aCGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                                 <  2:221595/149‑1 (MQ=255)
                                                                                                    cTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgca                                                                             <  3:195684/114‑1 (MQ=255)
                                                                                                    cTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgca                                                                             >  4:195684/1‑114 (MQ=255)
                                                                                                       ccccGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAAc                                                                        >  4:6022/1‑116 (MQ=255)
                                                                                                         ccGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCcgc                                     >  7:379340/1‑149 (MQ=255)
                                                                                                         ccGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCcgc                                     >  5:2508/1‑149 (MQ=255)
                                                                                                           gCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGcc                                                                                                                                              >  4:15165/1‑42 (MQ=255)
                                                                                                           gCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGcc                                                                                                                                              <  3:15165/42‑1 (MQ=255)
                                                                                                                 gCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGcga                             <  7:182/149‑1 (MQ=255)
                                                                                                                   cAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctg                                                                                                           <  5:211068/69‑1 (MQ=255)
                                                                                                                   cAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGACGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGcgac                            >  1:268961/1‑148 (MQ=255)
                                                                                                                   cAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGACGTCGGCCACCACCAGCAGGCCCAACgctg                                                                                                           >  6:211068/1‑69 (MQ=255)
                                                                                                                      ccAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGGAACTGGCCCGCGAAGGCGACGATg                        >  7:281978/1‑149 (MQ=255)
                                                                                                                           cgcgGTTGGCCATGGCCAGCCGCGTCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCAGCGAAGGCGACGATGCAgcg                   >  4:68850/1‑149 (MQ=255)
                                                                                                                            gcgGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGc                       >  5:357210/1‑144 (MQ=255)
                                                                                                                            gcgGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGc                       <  6:357210/144‑1 (MQ=255)
                                                                                                                               gCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTgg               >  4:340057/1‑149 (MQ=255)
                                                                                                                                 tggccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                                 >  4:287729/1‑113 (MQ=255)
                                                                                                                                 tggccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTgg                                                 <  3:287729/113‑1 (MQ=255)
                                                                                                                                    ccatggccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcgg          >  1:236726/1‑149 (MQ=255)
                                                                                                                                          ccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcggcgagat    >  3:59635/1‑145 (MQ=255)
                                                                                                                                          ccaACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTgcggc         <  4:59635/144‑1 (MQ=255)
                                                                                                                                            aaCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCt     >  4:343513/1‑146 (MQ=255)
                                                                                                                                            aaCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGa   <  3:343513/148‑1 (MQ=255)
                                                                                                                                            aaCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAg  >  3:151605/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCGCTGGTGCGGCGCCTGAG  >  NC_002947/3908276‑3908563

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