Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,103,656 Δ1 bp coding (1320/1374 nt) cpxA ← sensory histidine kinase in two‑component regulatory system with CpxR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,103,6540G.98.2% 242.9 / ‑3.3 56coding (1322/1374 nt)cpxAsensory histidine kinase in two‑component regulatory system with CpxR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (25/30);  minor base C (1/0);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGA  >  NC_000913/4103512‑4103791
                                                                                                                                              |                                                                                                                                         
ccTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTCTACAGCGGCAACCAAATCACCAGCCCTCACCCGCCCAGCg                                                                                                                                           >  2:843374/1‑143 (MQ=255)
   gATAAGACGCTATCAGCGTCGCATCAGGTATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCggct                                                                                                                                        <  1:452863/143‑3 (MQ=255)
    aTAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCcgcgcctg                                                                                                                                       >  1:900505/1‑139 (MQ=255)
       aGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGtct                                                                                                                                   <  2:19217/143‑1 (MQ=255)
         aCGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTAAAGTCAAGTTCGAAGTTTAACCCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGtcttc                                                                                                                                 <  1:1092917/143‑1 (MQ=37)
           tcTATCTGCGTCGCTTCAGGCATCCTGCCCAAGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTGCAAATGCGAAGCTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTg                                                                                                                               <  2:643467/142‑1 (MQ=37)
             tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGCTGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAACTCACCCGCCGTACACCGCCCCGC‑GGCCGCCTTCTGcc                                                                                                                             >  1:301015/1‑143 (MQ=37)
             tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGcc                                                                                                                             >  2:543494/1‑143 (MQ=255)
             tATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGcc                                                                                                                             <  2:516349/143‑1 (MQ=255)
              aTCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCt                                                                                                                            >  1:178804/1‑143 (MQ=255)
               tCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCtt                                                                                                                           >  1:1098557/1‑143 (MQ=255)
               tCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCtt                                                                                                                           <  1:18558/143‑1 (MQ=255)
                    gTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCcccc                                                                                                                      >  1:883540/1‑143 (MQ=255)
                    gTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCAccc                                                                                                                      <  2:1004822/143‑1 (MQ=255)
                       gCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGc                                                                                                                   <  1:861168/143‑1 (MQ=255)
                            aGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCCGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                              >  1:79585/1‑143 (MQ=255)
                            aGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                              <  1:610847/143‑1 (MQ=255)
                            aGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGa                                                                                                              <  1:71808/143‑1 (MQ=255)
                              gCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATCCAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTCCTGCCTTCACCCAGCCACGAtg                                                                                                            >  2:839811/1‑143 (MQ=255)
                                   cTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATgctgct                                                                                                       <  1:835848/143‑1 (MQ=255)
                                    tGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGctgctg                                                                                                      >  1:1026463/1‑143 (MQ=255)
                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGa                                                                                                     >  1:82700/1‑143 (MQ=255)
                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGa                                                                                                     <  1:845438/143‑1 (MQ=255)
                                     gCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCACTTATTCAGCGACAACCAAATCACCAGACGTAAAATCATCA‑CCGGCTATCATCTGCCTTCACCACGCCACTATACTGCTGa                                                                                                     >  2:650688/1‑143 (MQ=21)
                                        cAAACGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATg                                                                                                  <  1:978756/143‑1 (MQ=255)
                                          aaTGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGc                                                                                                >  2:306943/1‑143 (MQ=255)
                                               cGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTACCTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGttt                                                                                           <  2:346389/143‑1 (MQ=255)
                                               cGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATGCAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGttt                                                                                           <  2:299062/143‑1 (MQ=255)
                                                  atccggcgTCACCGCCTTATCCCGCCTGCAAATGCGACGTTTACCCCCGCTTCTACAGCGGCAACCAACTCCCCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTcaa                                                                                        <  1:62112/140‑1 (MQ=255)
                                                  aTGCGGCGTAAACGCCTTAGCCTGCCTGCAAATGCGAAGTTCAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTcaa                                                                                        <  1:1002549/143‑1 (MQ=255)
                                                      ggcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCCGCCACGATGCTGCTGAATGGCGGTTTCAACAAt                                                                                    >  1:598113/1‑143 (MQ=255)
                                                       gcgTAAACGCCTTATCCTGCCTGCAACTGCGAAGTTTAACTCCGCTTATACAGCGGCCACCAAATCACAAGCCGTAAACCGCCCCGC‑GGCGGTCTGCTGCCATCACCCAGCAGCGATGCTGCTGACTGCCGGTGTCAACAATc                                                                                   >  2:441961/1‑143 (MQ=255)
                                                       gcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATc                                                                                   <  2:96820/143‑1 (MQ=255)
                                                       gcgTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATc                                                                                   <  1:306943/143‑1 (MQ=255)
                                                                    acgctgcctgcAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAAc                                                                      <  2:797310/140‑1 (MQ=255)
                                                                    aTCCTGCCTGCAAATTCGACGTTGAACTCCGCTTACACAGCGGCAACCAAATCGCCAGCCGTAAACCGCCCCGC‑GGCTGTCTTCTGCCTTCACCCCGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAAc                                                                      <  2:805573/143‑1 (MQ=255)
                                                                     tCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAGcc                                                                     >  1:641141/1‑143 (MQ=255)
                                                                               aaaTGCGAAGGTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCa                                                           <  1:774020/143‑1 (MQ=255)
                                                                                           tAACTCCGCTTATACAGCGGCAACCAAATCACCACCCGTACACCGCCCAGC‑GCCCCTCTTCTGCCTTCACCCACCCACGATGCTGCTGAATGGCGGTCTTACCCATCGCCACCCCCAACCCTCTCACGCCAGACCCACCATcg                                               >  2:392879/1‑143 (MQ=255)
                                                                                               tCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATcgcggt                                           >  1:368023/1‑141 (MQ=255)
                                                                                               tCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACCACcgcgca                                           >  2:562834/1‑141 (MQ=255)
                                                                                                 cGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCt                                         <  1:285574/143‑1 (MQ=255)
                                                                                                 cGCTTAGACAGCGGCGCGCGAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTGCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCt                                         <  1:843374/143‑1 (MQ=255)
                                                                                                          aGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTa                                >  1:817542/1‑143 (MQ=255)
                                                                                                          aGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCCGAATGGCGGTTTCAACAATCCCCAGCCCCACACCAGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTa                                <  1:220300/143‑1 (MQ=255)
                                                                                                             ggCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGa                             <  1:277390/143‑1 (MQ=255)
                                                                                                              gCAACCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTAGCGCCAGCTTCACGATCGCGTGCTCCATCGGCCCGAt                            >  2:508482/1‑143 (MQ=255)
                                                                                                                aaCCAAATCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCCGTCTCAACAATCGCCAGCCCCAAAGCTGGACCGCCAGACTCACGATCGCGTCCTCCCTCGGCACGCTac                          >  2:159438/1‑142 (MQ=255)
                                                                                                                    aaaTCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAacg                      <  1:1035624/143‑1 (MQ=255)
                                                                                                                     aaTCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAacgg                     >  1:985418/1‑143 (MQ=255)
                                                                                                                     aaTCACCAGCCGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATACTGCTCGATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGACAGATTCACGATCGCACGCTTCCTCGGTACGATAGAacgc                     >  2:130684/1‑142 (MQ=255)
                                                                                                                               cGTAAACCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATc           <  1:1040193/143‑1 (MQ=255)
                                                                                                                                  aaaCCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGt        <  2:631227/143‑1 (MQ=255)
                                                                                                                                    aCCGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGCGGGAGAATCGGtgc      >  1:675791/1‑141 (MQ=255)
                                                                                                                                     ccGCCCAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTcg     >  1:1005995/1‑143 (MQ=255)
                                                                                                                                      cGCCCAGC‑GGCTGTCCTCTGCCCCCACCCAGCCACGATGCTGCTGAATGGCTGTTTCACCGACCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTcgc    <  1:405697/143‑1 (MQ=255)
                                                                                                                                        cccAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGa  <  1:439668/143‑1 (MQ=255)
                                                                                                                                        cccAGC‑GGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGa  >  1:644420/1‑143 (MQ=255)
                                                                                                                                        cccAGC‑GGCTGGCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTGCCACAATCGCCAGCCCCAATCCTGTACCGCCAGCTTAACGATCGTGATTGTCATCGGTTCGATAGCACGGGCGGAGCATCTGTGcgcgg  >  1:40687/1‑142 (MQ=255)
                                                                                                                                              |                                                                                                                                         
CCTGATAAGACGCTATCAGCGTCGCATCAGGCATCCTGCTCAAATGCCGGATGCGGCGTAAACGCCTTATCCTGCCTGCAAATGCGAAGTTTAACTCCGCTTATACAGCGGCAACCAAATCACCAGCCGTAAACCGCCCAGCGGGCTGTCTTCTGCCTTCACCCAGCCACGATGCTGCTGAATGGCGGTTTCAACAATCGCCAGCCCCAAACCTGTACCGCCAGATTCACGATCGCGTGCTTCATCGGTACGATAGAACGGACGGAAAATCTGTTCGCGA  >  NC_000913/4103512‑4103791

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑