Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,241,918 C→T 7.0% V538V (GTC→GTT PP_2838 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,241,9180CT 7.0% 167.1 / 2.6 57V538V (GTC→GTTPP_2838hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base T (2/2);  total (28/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 = 8.44e-01

ACCGTGCCGGCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTATCGAAGGGCATCAGCG  >  NC_002947/3241776‑3242066
                                                                                                                                               |                                                                                                                                                     
aCCGTGCCGGCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGgtg                                                                                                                                                 >  4:88376/1‑149 (MQ=255)
aCCGTGCCGGCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGgtg                                                                                                                                                 >  6:242386/1‑149 (MQ=255)
 ccGTGCCGGCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGgtgt                                                                                                                                                <  5:252834/149‑1 (MQ=255)
      ccGGCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCAc                                                                                                                                           <  3:88376/149‑1 (MQ=255)
        ggCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGt                                                                                                                                         >  7:40789/1‑149 (MQ=255)
          cATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTgc                                                                                                                                       <  5:1959/149‑1 (MQ=255)
               ccGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTg                                                                                                                                   >  6:3243/1‑148 (MQ=255)
                             ttCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTg                                                                                                                    <  8:40789/149‑1 (MQ=255)
                               cTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTTGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAg                                                                                                                  >  7:224944/1‑149 (MQ=255)
                               cTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTTGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAg                                                                                                                  >  7:230985/1‑149 (MQ=255)
                                    ccGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTa                                                                                                               >  5:18953/1‑147 (MQ=255)
                                              ccGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGt                                                                                                   >  1:188391/1‑149 (MQ=255)
                                                    tCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCGGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTgcgcgc                                                                                             <  5:381466/149‑1 (MQ=255)
                                                      gCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCt                                                                                           >  2:286064/1‑149 (MQ=255)
                                                       cATGGTCGAGGTCAAGGGGCCCGGTGAACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCAc                                                                                                                                           >  1:4473/1‑101 (MQ=255)
                                                       cATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGTCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCtt                                                                                          >  8:369471/1‑149 (MQ=255)
                                                       cATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGTCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCtt                                                                                          >  8:374864/1‑149 (MQ=255)
                                                       cATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCAc                                                                                                                                           <  2:4473/100‑1 (MQ=255)
                                                       cATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCtt                                                                                          <  5:3243/149‑1 (MQ=255)
                                                         tGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTtgt                                                                                        <  8:301779/149‑1 (MQ=255)
                                                         tGGTCGAGGTAACGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGGCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTtgt                                                                                        <  2:263548/149‑1 (MQ=255)
                                                           gTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTtgtg                                                                                       <  5:34241/148‑1 (MQ=255)
                                                            tCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTgcgc                                                                                               <  4:195520/139‑1 (MQ=255)
                                                            tCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTgcgc                                                                                               >  3:195520/1‑139 (MQ=255)
                                                            tCGAGGTCAAGGGGCCCGGTGACAGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTtgtgt                                                                                      >  2:62994/1‑148 (MQ=255)
                                                               aGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGt                                                                                  <  6:18953/149‑1 (MQ=255)
                                                                 gTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGtt                                                                                 >  3:95841/1‑148 (MQ=255)
                                                                    aaGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAgc                                                                             >  4:7105/1‑149 (MQ=255)
                                                                     aGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAgc                                                                             >  6:46316/1‑148 (MQ=255)
                                                                     aGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAgc                                                                             <  5:46316/148‑1 (MQ=255)
                                                                                    cGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTTGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgac                                                             <  8:224944/149‑1 (MQ=255)
                                                                                    cGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTTGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgac                                                             <  8:230985/149‑1 (MQ=255)
                                                                                        tGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGTCCTGCCGGTCGCGGTGTGCAACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCACGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgacctg                                                          >  5:331173/1‑148 (MQ=255)
                                                                                                aaccaacTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTgcgc                                                 >  6:81913/1‑149 (MQ=255)
                                                                                                   caacTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTc                                              <  3:7105/149‑1 (MQ=255)
                                                                                                     acTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCAc                                            <  4:95841/149‑1 (MQ=255)
                                                                                                              gctggcGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGcc                                   >  6:201414/1‑149 (MQ=255)
                                                                                                                tggcGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgacctggac                                                       >  5:204007/1‑127 (MQ=255)
                                                                                                                tggcGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgacctggac                                                       <  6:204007/127‑1 (MQ=255)
                                                                                                                  gcGTTCTGCAGGCAACACGGCCTGTCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCAc                                <  7:374864/148‑1 (MQ=255)
                                                                                                                  gcGTTCTGCAGGCAACACGGCCTGTCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCAc                                <  7:369471/148‑1 (MQ=255)
                                                                                                                   cGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCg                              >  6:90703/1‑149 (MQ=255)
                                                                                                                    gTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgacctgg                                                         <  4:274865/121‑1 (MQ=255)
                                                                                                                    gTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCgacctgg                                                         >  3:274865/1‑121 (MQ=255)
                                                                                                                         gCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCag                         >  2:128234/1‑148 (MQ=255)
                                                                                                                          cAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCag                         <  1:286064/148‑1 (MQ=255)
                                                                                                                            ggCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCaggagg                     <  5:81913/149‑1 (MQ=255)
                                                                                                                             gCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGc                                                                <  7:333409/105‑1 (MQ=255)
                                                                                                                             gCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGc                                                                >  8:333409/1‑105 (MQ=255)
                                                                                                                               aaCACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTAGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTa                  >  4:273231/1‑149 (MQ=255)
                                                                                                                                 cacGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTtgtgt                                                                                      <  3:120223/80‑1 (MQ=255)
                                                                                                                                 cacGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTtgtgt                                                                                      >  4:120223/1‑80 (MQ=255)
                                                                                                                                  acGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTATCg               <  5:90703/149‑1 (MQ=255)
                                                                                                                                      ccTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTATCGAAgg           <  5:242386/149‑1 (MQ=255)
                                                                                                                                       cTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTATCGAAggg          <  2:188391/149‑1 (MQ=255)
                                                                                                                                              tCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCagga                       >  2:198013/1‑129 (MQ=255)
                                                                                                                                              tCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCagga                       <  1:198013/129‑1 (MQ=255)
                                                                                                                                               cgcgGTCTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGCACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTATCGAAGGGCATCAGCg  >  6:32256/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
ACCGTGCCGGCATGCCCGACCTGATTCAGTTCTGGCCCGAGCAGGGCCGTTATCGCATGGTCGAGGTCAAGGGGCCCGGTGACCGCCTGCAGGACAACCAACTGCGCTGGCTGGCGTTCTGCAGGCAACACGGCCTGCCGGTCGCGGTGTGCCACGTGCGCTGGAGCGAAACACCTTGAGCTACAGCGTAGCGGTGCGCGCCTTGTGTGAGTTCAGCGCCAAGGTCGGCGACCTGGACCTGCGCTTCACGCCATCGCCCACCGCCCAGGAGGGTATCGAAGGGCATCAGCG  >  NC_002947/3241776‑3242066

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: