Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,111,977 T→A 8.2% H118Q (CAT→CAA)  PP_5362 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,111,9770TA 8.2% 159.7 / 3.5 57H118Q (CAT→CAA) PP_5362hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base A (2/3);  total (26/31)
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.57e-01

CACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCC  >  NC_002947/6111832‑6112119
                                                                                                                                                 |                                                                                                                                              
cacaAACCAACAGCCGCTATGGCACATTGACTTTTTTTACAGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTg                                                                                                                                              >  5:222092/1‑148 (MQ=255)
   aaaCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTCGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGc                                                                                                                                          >  5:69880/1‑149 (MQ=255)
     aCCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCtt                                                                                                                                        <  6:335828/149‑1 (MQ=255)
        aaCAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGaa                                                                                                                                     <  2:70604/149‑2 (MQ=255)
          cAGCGGCTTTGGCACATTTAGTTTTTTTACCGGCGTGGGCAGCGATGGTGAGGCCTGGTTTGGCTACTGCCGCCTCGTAGGATGGACCACCTTGAACCGCGTATCCTGTTCGAGCGGGGAGCGCTTCTTGCTGCAATGCAGCTTGaa                                                                                                                                     <  8:249774/147‑2 (MQ=255)
            gCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                 <  2:193783/149‑1 (MQ=255)
                cTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGt                                                                                                                             <  7:285761/149‑1 (MQ=255)
                cTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGt                                                                                                                             >  7:343397/1‑149 (MQ=255)
                   tGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGcc                                                                                                                          <  1:153138/149‑1 (MQ=255)
                     gCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCg                                                                                                                         <  3:89847/148‑1 (MQ=255)
                     gCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGt                                                                                                                        >  2:136417/1‑149 (MQ=255)
                       acaTTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGAGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCATGAACTAGGGGTGCCGtt                                                                                                                       >  8:253606/1‑148 (MQ=255)
                              cTTTTTTTACCGGCCTGGCCAGCGATGGTGTGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACGTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcct                                                                                                                <  4:289512/148‑1 (MQ=255)
                              cTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctc                                                                                                               <  8:343838/149‑1 (MQ=255)
                               tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                 <  6:8088/130‑1 (MQ=255)
                               tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAg                                                                                                                                 >  5:8088/1‑130 (MQ=255)
                               tttttttACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGcctcc                                                                                                              >  1:381199/1‑149 (MQ=255)
                                        cGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtttt                                                                                                     >  8:271684/1‑149 (MQ=255)
                                        cGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCAtttt                                                                                                     <  8:343397/149‑1 (MQ=255)
                                                ccAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttt                                                                                             >  8:251892/1‑149 (MQ=255)
                                                ccAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGCGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGttt                                                                                             <  2:315026/149‑1 (MQ=255)
                                                 cAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTTGACCAGCTTGATCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttt                                                                                            <  5:98522/149‑1 (MQ=255)
                                                        ggTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctg                                                                                     >  5:268024/1‑149 (MQ=255)
                                                          tGAGGCCTTGGTGGGCTACAGCCGCCGCGTAGGCTGGTCCACGTTGAAGCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTctgc                                                                                    <  5:108110/148‑1 (MQ=255)
                                                            tggCCTGGTTGGGCTACAGCCGCCGCGTAGGCTGGAACACCTTGAAGCGCGTTTCCAGTTCCAGCGGGGAGCGCTTCAATCTTCATTTCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTTAACGTTTTTTCTCtgctgc                                                                                 <  6:222092/148‑1 (MQ=255)
                                                              gCCTGGGTGGGCTACAGCCGCCGAGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGcc                                                                                >  6:199993/1‑148 (MQ=255)
                                                                cTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAt                                                                             >  4:364657/1‑149 (MQ=255)
                                                                cTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAt                                                                             >  4:247579/1‑149 (MQ=255)
                                                                   ggtgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTAtttt                                                                          <  7:251892/149‑1 (MQ=255)
                                                                     tgggCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGt                                                                        <  2:381199/149‑1 (MQ=255)
                                                                          tACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGAGGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGt                                                                    >  4:256435/1‑148 (MQ=255)
                                                                            cAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCa                                                                  <  5:199993/148‑1 (MQ=255)
                                                                               ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                          <  8:20882/121‑1 (MQ=255)
                                                                               ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGtttttt                                                                                          >  7:20882/1‑121 (MQ=255)
                                                                               ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                              >  3:16167/1‑149 (MQ=255)
                                                                               ccgccgCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCtt                                                              >  3:17952/1‑149 (MQ=255)
                                                                                ctccgcGTCGGCTGGACCACCGTGAAACGCGTATCCAGTTCCAGCGGGGAGCGCGTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGGTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAACTCAGCTTc                                                             <  6:391222/147‑1 (MQ=255)
                                                                                cgccgcGTAGGCTTGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGTAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGATCGTTTTTTctct                                                                                      <  8:23902/124‑1 (MQ=255)
                                                                                cgccgcGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTctct                                                                                      >  7:23902/1‑124 (MQ=255)
                                                                                        agccTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCAGCCTATTTTGTAAGTCAGCTTCGAGCAcgc                                                     >  2:284305/4‑149 (MQ=255)
                                                                                          gCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACgcgc                                                   >  6:203012/1‑149 (MQ=255)
                                                                                           cTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCt                                                  <  4:17952/149‑1 (MQ=255)
                                                                                              gACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGaa                                               <  7:233990/149‑1 (MQ=255)
                                                                                                 caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                           <  2:174245/134‑1 (MQ=255)
                                                                                                 caccTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGa                                                           >  1:174245/1‑134 (MQ=255)
                                                                                                       gAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCa                                        <  3:247579/147‑1 (MQ=255)
                                                                                                        aaCCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGc                                      <  5:203012/148‑1 (MQ=255)
                                                                                                                aTCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCaa                             >  1:139143/1‑149 (MQ=255)
                                                                                                                 tCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCaa                             <  4:16167/148‑1 (MQ=255)
                                                                                                                     gTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGA‑GGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTATGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCttt                       <  6:69880/149‑1 (MQ=255)
                                                                                                                       tCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGTGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCttt                       <  4:115572/148‑1 (MQ=255)
                                                                                                                          aGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGTTGCCGTTGACGCCTCCTTTCATTATGAAAGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCa                   <  7:253606/149‑1 (MQ=255)
                                                                                                                          aGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGt                     >  2:315947/1‑147 (MQ=255)
                                                                                                                              gggAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGAcc               >  6:385106/1‑149 (MQ=255)
                                                                                                                                  gcgcTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGGCCAACCCTTTGTCAGCCCGaa            >  8:60559/1‑148 (MQ=255)
                                                                                                                                  ctcgatCAAGCTGCAATGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGaa            >  5:247574/6‑148 (MQ=255)
                                                                                                                                           gCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGcc  <  6:268024/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CACAAACCAACAGCCGCTATGGCACATTGACTTTTTTTACCGGCCTGGCCAGCGATGGTGAGGCCTGGGTGGGCTACAGCCGCCGCGTAGGCTGGACCACCTTGAACCGCGTATCCAGTTCCAGCGGGGAGCGCTTCAAGCTGCATTGCAGCTTGATCGAGGGGTGCCGTTGACGCCTCCTTTCATTTTGAACGTTTTTTCTCTGCTGCCTATTTTGTAAGTCAGCTTCGAGCACGCGCTGAAAACACCAGCCTTGACCAACCCTTTGTCAGACCGAAAAGCTATGCC  >  NC_002947/6111832‑6112119

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: ‑