Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,670 C→T 8.5% intergenic (‑69/‑144) PP_4522 ← / → PP_4523 LysR family transcriptional regulator/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,137,6700CT 8.5% 167.5 / 3.4 58intergenic (‑69/‑144)PP_4522/PP_4523LysR family transcriptional regulator/agmatinase
Reads supporting (aligned to +/- strand):  ref base C (30/23);  new base T (3/2);  total (33/25)
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.75e-01

GTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACC  >  NC_002947/5137527‑5137808
                                                                                                                                               |                                                                                                                                          
gTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGaa                                                                                                                                       >  5:277924/1‑149 (MQ=255)
 ttGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAAc                                                                                                                                      <  2:261658/149‑1 (MQ=255)
    gCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACACGCCTTTAGTTCAGAAATTCCTGAACTaa                                                                                                                                   >  8:218354/1‑149 (MQ=255)
     cAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAg                                                                                                                                  >  5:344876/1‑149 (MQ=255)
     cAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAg                                                                                                                                  >  7:92486/1‑149 (MQ=255)
          cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAttt                                                                                                                             <  4:241156/149‑1 (MQ=255)
          cccTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTAttt                                                                                                                             >  3:241156/1‑149 (MQ=255)
               gtggCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCg                                                                                                                         <  8:7207/148‑1 (MQ=255)
                    gcACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaaggc                                                                                                                               <  7:140882/137‑5 (MQ=255)
                    gcACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGGAATTTCTGAACTaaaggc                                                                                                                               >  8:140882/1‑133 (MQ=255)
                    gcACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGc                                                                                                                    >  5:264148/1‑148 (MQ=255)
                    gcACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCg                                                                                                                   <  2:285005/149‑1 (MQ=255)
                       ccacGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAtt                                                                                                                >  4:196024/1‑149 (MQ=255)
                        cacGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttc                                                                                                               <  7:38676/149‑1 (MQ=255)
                                 gaAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGgc                                                                                                      >  4:168764/1‑149 (MQ=255)
                                     aTGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGcgctc                                                                                                  <  3:168764/149‑1 (MQ=255)
                                           agcagTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCtt                                                                                              >  1:270661/1‑147 (MQ=255)
                                           agcagTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGa                                                                                            <  6:277924/149‑1 (MQ=255)
                                              agTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCa                                                                                         <  6:195520/149‑1 (MQ=255)
                                                tttCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCa                                                                                         <  1:131541/147‑1 (MQ=255)
                                                    aGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTg                                                                                   <  3:196024/149‑1 (MQ=255)
                                                     ggTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTgg                                                                                  >  5:207604/1‑149 (MQ=255)
                                                      gTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttcttc                                                                                                            >  6:22701/1‑122 (MQ=255)
                                                      gTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttcttc                                                                                                            <  5:22701/122‑1 (MQ=255)
                                                            cAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttc                                                                                                               <  4:28692/113‑1 (MQ=255)
                                                            cAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGAttc                                                                                                               >  3:28692/1‑113 (MQ=255)
                                                            cAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAAt                                                                           <  6:207604/149‑1 (MQ=255)
                                                                gTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTcc                                                                       <  2:307261/149‑1 (MQ=255)
                                                                      gCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGccc                                                                 <  6:264148/149‑1 (MQ=255)
                                                                        cATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgc                                                               >  4:358012/1‑149 (MQ=255)
                                                                         aTGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCccgcg                                                              >  2:158421/1‑148 (MQ=255)
                                                                                aaCGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAAc                                                       <  4:294087/149‑1 (MQ=255)
                                                                                aaCGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAAc                                                       <  2:270661/149‑1 (MQ=255)
                                                                                      tCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACg                                                 >  2:283880/1‑149 (MQ=255)
                                                                                         acaAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGa                                              <  1:283880/149‑1 (MQ=255)
                                                                                           aaGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTcc                                                                       >  4:29360/1‑122 (MQ=255)
                                                                                           aaGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTcc                                                                       <  3:29360/122‑1 (MQ=255)
                                                                                             ggcaaagtgCGCGCAGTATAACTCCTCTACACGCCTTTACATCAGAAATTCCTGAACTAAGTATTAGTCGGTAGCGATTCTACCCAGGCGCACGCTTGAGCAGACTTCGCCGAAATCTCCCTGCACGCCGGCCAACCACACGAGAATgg                                          >  8:110878/7‑149 (MQ=255)
                                                                                                     gcgcgcAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCgcggc                                  >  7:230331/1‑149 (MQ=255)
                                                                                                         gcAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGca                              >  2:75166/1‑149 (MQ=255)
                                                                                                           aGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGca                              >  7:237389/1‑147 (MQ=255)
                                                                                                               tAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGtt                        >  4:47024/1‑149 (MQ=255)
                                                                                                                aaCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTc                       >  1:262976/1‑149 (MQ=255)
                                                                                                                aaCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCAGGCAACCAACACGTGAATGGCGCGCGGAGGCACAGGTTc                       >  3:229589/1‑149 (MQ=255)
                                                                                                                 aCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTc                       <  4:173968/148‑1 (MQ=255)
                                                                                                                 aCTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTc                       >  8:42134/1‑148 (MQ=255)
                                                                                                                     ctctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaac                   <  3:358012/148‑1 (MQ=255)
                                                                                                                     ctctACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaac                   >  7:158671/1‑148 (MQ=255)
                                                                                                                             gCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                <  6:158132/143‑1 (MQ=255)
                                                                                                                             gCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaa                >  5:158132/1‑143 (MQ=255)
                                                                                                                              tacttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaaca              <  2:200916/141‑1 (MQ=255)
                                                                                                                              tacttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGaacaaaca              >  1:200916/4‑144 (MQ=255)
                                                                                                                                 ttAGTTCAGGAATTTCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATc      >  1:6316/1‑149 (MQ=255)
                                                                                                                                     ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgacc  >  3:356796/1‑149 (MQ=255)
                                                                                                                                     ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgacc  >  3:262719/1‑149 (MQ=255)
                                                                                                                                     ttCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATcgacc  >  2:254508/1‑149 (MQ=255)
                                                                                                                                          aaaTTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGcac                             <  7:267194/117‑1 (MQ=255)
                                                                                                                                          aaaTTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGcac                             >  8:267194/1‑117 (MQ=255)
                                                                                                                                               |                                                                                                                                          
GTTGGCAAAGCCCTGGTGGCGCACCACGCTGATGAAGATGCGCAGCAGTTTCAGGTCGGGCAAGGTCGAGGCCATGGTGCAACGGTTCCACAAGGGTATGTGCGCGCAGTCTAACTGCTCTACCCGCCTTTAGTTCAGAAATTCCTGAACTAAGTATTTGTCGGTAGCGATTCTTCCCCGGCGCTCGCTTGAGCAGACTTGGCCGAAATGTCCCTGCCCGCCGGCAAACCACACGTGAATGGCGCGCGGCGGCACAGGTTCGAACAAACACAACAATCGACC  >  NC_002947/5137527‑5137808

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