Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,366,654 C→T 7.0% G16S (GGC→AGC)  PP_1191 ← S4 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,366,6540CT 7.0% 162.4 / 2.6 57G16S (GGC→AGC) PP_1191S4 domain‑containing 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 = 9.09e-01

ACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGACACATAG  >  NC_002947/1366509‑1366796
                                                                                                                                                 |                                                                                                                                               
aCCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTCAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAg                                                                                                                                              <  7:6872/149‑1 (MQ=255)
  cGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAgct                                                                                                                                            >  1:94965/1‑149 (MQ=255)
        cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTAGCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTc                                                                                                                                        <  7:292398/147‑1 (MQ=255)
        cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcc                                                                                                                                       <  8:20115/148‑2 (MQ=255)
        cAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGa                                                                                                                                      >  1:193632/1‑149 (MQ=255)
         aGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGGCCTGCGCGAt                                                                                                                                     >  3:1369/1‑149 (MQ=255)
            gTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCTCGGCTGCTcccgccg                                                                                                                                  <  3:81137/149‑7 (MQ=255)
                    aCGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTg                                                                                                                           >  8:92609/1‑148 (MQ=255)
                          gcCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGc                                                                                                                    <  5:273383/149‑1 (MQ=255)
                           cccGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCg                                                                                                                   <  8:273677/149‑1 (MQ=255)
                             cGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGa                                                                                                                 <  2:193632/149‑1 (MQ=255)
                              ggCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGAc                                                                                                                >  5:123555/1‑149 (MQ=255)
                               gcagcaGCGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACg                                                                                                               >  2:287336/1‑149 (MQ=255)
                                 agcagcGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACggg                                                                                                             >  5:150433/1‑149 (MQ=255)
                                  gcagcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt                                                                                                            >  1:233835/1‑149 (MQ=255)
                                  gcagcGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGt                                                                                                            <  1:287336/149‑1 (MQ=255)
                                    agcGCGTCGAGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTc                                                                                                          <  6:150433/149‑1 (MQ=255)
                                     gcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg                                                                                                         >  5:233314/1‑149 (MQ=255)
                                     gcTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCg                                                                                                         <  1:135052/149‑1 (MQ=255)
                                                   tgctCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAAtt                                                                                          <  2:94965/149‑1 (MQ=255)
                                                        aaCCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATgg                                                                                                 >  6:38800/1‑136 (MQ=255)
                                                        aaCCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATgg                                                                                                 <  5:38800/136‑1 (MQ=255)
                                                         aCCTTGAACTGCGGCTGCTCGACCACCACGACATCGACCGTTACCCAGCCTCCTTCAATGTTCTGCTCTGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCTCAAg                                                                                    <  1:194465/149‑1 (MQ=255)
                                                          ccTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAgg                                                                                   >  6:110476/1‑149 (MQ=255)
                                                              gAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACACAGCCTCCATCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGCCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGAAAGGCAgg                                                                               >  4:69752/1‑149 (MQ=255)
                                                              gAACTGCGGCTGCTCGACCACCACGCCATAGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGTCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAg                                                                                <  7:17221/148‑1 (MQ=255)
                                                                 cTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGaa                                                                            <  1:106498/149‑1 (MQ=255)
                                                                             gACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGATGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTg                                                                <  7:70475/149‑1 (MQ=255)
                                                                             gACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCAtt                                                                 <  5:110476/148‑1 (MQ=255)
                                                                             gACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCAt                                                                  <  2:233835/147‑1 (MQ=255)
                                                                              accaccACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGt                                                               <  7:143035/149‑1 (MQ=255)
                                                                                   cacGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCAt                                                                  <  1:197600/141‑1 (MQ=255)
                                                                                   cacGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCAt                                                                  >  2:197600/1‑141 (MQ=255)
                                                                                            gACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGc                                                 >  6:171958/1‑149 (MQ=255)
                                                                                              ccGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGAGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCgg                                               <  3:69752/149‑1 (MQ=255)
                                                                                                    cccAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCCAGCTGCTCCATGAGGCGTTTTGACAGGCGGACGGGTTCGGACATGGGGCAATTAGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAgat                                          >  1:165570/1‑148 (MQ=255)
                                                                                                     ccAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGGGCCGCCCCCCGGGGGGATGTGGCGGTTGGTGAGGCGGCCGGGGGGGGGCAGTGGGGGAATATCAAAGCAGGGGAAAGCGACGATAGTAACAGCCGGCAAGCGGTAgatta                                        >  2:319868/1‑147 (MQ=255)
                                                                                                      cAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAgatgat                                       <  7:92609/149‑1 (MQ=255)
                                                                                                      cAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAgatgat                                       >  4:74093/1‑149 (MQ=255)
                                                                                                        gCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATgg                                     <  5:171958/149‑1 (MQ=255)
                                                                                                                   tGTACTGCTCGGCCTCGGGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGGTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGGGGt                                              <  5:290144/129‑1 (MQ=255)
                                                                                                                   tGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTc                          <  4:1369/149‑1 (MQ=255)
                                                                                                                   agTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCACGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAACGCAGGGAAAGGCCTCCATTTTACCAGCCAGCAAGCCGt                                              >  6:290144/2‑129 (MQ=255)
                                                                                                                         gCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATgc                    >  5:18069/1‑149 (MQ=255)
                                                                                                                          cTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTAGCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATgc                    <  2:135850/148‑1 (MQ=255)
                                                                                                                          cTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATgcg                   >  6:302967/1‑149 (MQ=255)
                                                                                                                          cTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGAAGAGGCGTTTGGACAGGCGGACGGGTTCGGAGATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATgcg                   >  4:297686/1‑149 (MQ=255)
                                                                                                                             ggCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCa                >  8:145166/1‑149 (MQ=255)
                                                                                                                              gCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCAt               >  4:213680/1‑149 (MQ=255)
                                                                                                                                 tCGCGGCGGGGGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCAtttt            >  1:321246/1‑149 (MQ=255)
                                                                                                                                   gcggcgGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCAtttt            <  5:302967/147‑1 (MQ=255)
                                                                                                                                   gcggcgGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTg           <  6:123555/148‑1 (MQ=255)
                                                                                                                                     tcatcgAGCAGCTCGGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGa         >  5:79925/6‑148 (MQ=255)
                                                                                                                                       atcgAGCAGCTCGGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGaca       >  8:114581/4‑148 (MQ=255)
                                                                                                                                         ggggcagcCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGACACAt    <  2:321246/146‑1 (MQ=255)
                                                                                                                                          gAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGAcaca     >  3:109001/1‑147 (MQ=255)
                                                                                                                                           agcagcCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGACACATAg  >  2:223739/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACCGGCTCCAGGGTTTCGGCACGGGCGCCCGGCAGCAGCTCGACGCGCTGCTGCTCAACCTTGAACTGCGGCTGCTCGACCACCACGCCATCGACCGTTACCCAGCCTCCTTCAATGTACTGCTCGGCCTCGCGGCGGGAGCAGCCGAGCTGCTCGATGAGGCGTTTGGACAGGCGGACGGGTTCGGACATGGGGCAATTCGCAAGGCAGGGAAAGGCCTCCATTGTACCTGCCTGCAAGCGGTAGATGATGGGAATCTGTGTCAGATGCGCCATTTTGGACACATAG  >  NC_002947/1366509‑1366796

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