Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,428,464 C→A 9.2% G101V (GGC→GTC)  sohB ← peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,428,4640CA 9.2% 163.9 / 4.7 53G101V (GGC→GTC) sohBpeptidase
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base A (2/3);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCA  >  NC_002947/4428317‑4428602
                                                                                                                                                   |                                                                                                                                           
gcCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGc                                                                                                                                             >  4:336302/1‑148 (MQ=255)
       gCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGGCGACATCaaa                                                                                                                                      <  5:181244/148‑1 (MQ=255)
         gACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                   <  7:263954/149‑1 (MQ=255)
          aCGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTc                                                                                                                                   >  4:205744/1‑148 (MQ=255)
           cGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGGCGACATCAAAgcctc                                                                                                                                 <  4:70140/149‑5 (MQ=255)
              aGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTa                                                                                                                              <  1:327741/149‑1 (MQ=255)
                  aaGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGCGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGaa                                                                                                                         >  6:272206/1‑149 (MQ=255)
                     ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATag                                                                                                                      <  4:21836/149‑2 (MQ=255)
                     ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAc                                                                                                                      <  3:26285/149‑1 (MQ=255)
                     ccGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGAACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATAc                                                                                                                      <  5:296092/149‑1 (MQ=255)
                            tgtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                               >  3:4566/1‑149 (MQ=255)
                             gtgCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctct                                                                                                              >  4:98151/1‑149 (MQ=255)
                                cACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGt                                                                                                           <  4:4566/149‑1 (MQ=255)
                                cACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGt                                                                                                           >  5:161109/1‑149 (MQ=255)
                                 aCCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTc                                                                                                          <  3:98151/149‑1 (MQ=255)
                                 aCCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTc                                                                                                          >  3:241499/1‑149 (MQ=255)
                                   cAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCt                                                                                                        >  2:127639/1‑149 (MQ=255)
                                           cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGTCTCTTGTCCTCGGCCTTg                                                                                                <  6:161109/149‑1 (MQ=255)
                                           cgccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTg                                                                                                >  7:308083/1‑149 (MQ=255)
                                            gccgcTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGc                                                                                               <  4:241499/149‑1 (MQ=255)
                                                     cAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTtcttct                                                                                      >  4:67246/1‑149 (MQ=255)
                                                                 cACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCg                                                                          >  6:50102/1‑149 (MQ=255)
                                                                  aCTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGc                                                                         >  5:284444/1‑149 (MQ=255)
                                                                  aCTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGc                                                                         >  3:126599/1‑149 (MQ=255)
                                                                       gTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCGTGCCCTTCTTCTGCTGTTTTTCCGCCTTgg                                                                    <  5:272206/149‑1 (MQ=255)
                                                                        tCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGATTTTCCGCCTTGGc                                                                   <  1:127639/149‑1 (MQ=255)
                                                                         cACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGcc                                                                  >  3:316430/1‑149 (MQ=255)
                                                                          aCGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                               <  6:9745/103‑1 (MQ=255)
                                                                          aCGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGctc                                                                                                               >  5:9745/1‑103 (MQ=255)
                                                                           cgcgcgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCctg                                                                <  6:20169/149‑1 (MQ=255)
                                                                               cgGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCtgctgc                                                            <  8:308083/149‑1 (MQ=255)
                                                                                     cgagcgTCAGCAGTGCGGTGATTTCGGTGCGCAGGCTTTCGCGGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTgcgc                                                      <  5:294522/149‑1 (MQ=255)
                                                                                             agcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                              >  2:184807/1‑149 (MQ=255)
                                                                                             agcagTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCtt                                              >  6:263040/1‑149 (MQ=255)
                                                                                                    cGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                       >  5:154558/1‑149 (MQ=255)
                                                                                                    cGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                       <  5:50102/149‑1 (MQ=255)
                                                                                                        gATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCt                                   >  7:33958/1‑149 (MQ=255)
                                                                                                          tttCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCtgc                                 >  2:8816/1‑147 (MQ=255)
                                                                                                             cGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTc                               <  8:33958/148‑1 (MQ=255)
                                                                                                              gTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCg                              <  6:154558/148‑1 (MQ=255)
                                                                                                                  cgcAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAg                         <  3:205744/149‑1 (MQ=255)
                                                                                                                        cTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTCAGCTGGGCCTGCTCGACAAGGcccag                   <  1:8816/149‑4 (MQ=255)
                                                                                                                         tttCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGAc                  <  3:336302/149‑1 (MQ=255)
                                                                                                                         tttCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGAc                  >  7:94728/1‑149 (MQ=255)
                                                                                                                          ttCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTCACCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACt                 >  3:291676/1‑149 (MQ=255)
                                                                                                                           tCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTc                >  1:247379/1‑149 (MQ=255)
                                                                                                                              gTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAg             <  5:263040/149‑1 (MQ=255)
                                                                                                                                   cGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGgcgc         >  1:220686/1‑149 (MQ=255)
                                                                                                                                    gAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGta                               >  1:212150/1‑124 (MQ=255)
                                                                                                                                    gAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTc                               <  2:212150/125‑1 (MQ=255)
                                                                                                                                        cTTTGATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTCCGCAGGGCCTTGAGCTGTGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg    >  5:340161/1‑148 (MQ=255)
                                                                                                                                         tttGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg    >  3:410/1‑147 (MQ=255)
                                                                                                                                         tttGATGGCGACATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTgg                                       >  2:119434/1‑112 (MQ=255)
                                                                                                                                           tGATGTCGACATCAAAGTCCATTACGAATACCCGGCTCTTGTCCTCTGCCTTGCCCTTCTTCTGATTTTTTTCCGCCTTTGCCTGCTGCTTGCGCAGTGCCTTGATATTTTCATTTTCGAGCAGGCCGGACTCCATTCGCTCACGCa  <  7:164340/147‑1 (MQ=255)
                                                                                                                                                   ccATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg    <  1:60610/137‑1 (MQ=255)
                                                                                                                                                   ccATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACg    >  2:60610/1‑137 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
GCCAGTTGCGACGCAGCCAAGCCGTAACTGTGCACCAGGCCACCGCCGCTTTCCAGGCGCAGTACCACTTCGTCACGCGCGGTAGCGAGCGTCAGCAGTGCGGTGATTTCGTTGCGCAGGCTTTCGGTGGCCGAGGCTTTGATGTCGCCATCAAAGTCCAGTACGAATACCCGGCTCTTGTCCTCGGCCTTGCCCTTCTTCTGCTGTTTTTCCGCCTTGGCCTGCTGCTTGCGCAGGGCCTTGAGCTGGGCCTTGTCGAGCAGGCCGGACTCCAGGCGCTCACGCA  >  NC_002947/4428317‑4428602

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: