Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA48.1% 67.6 / 23.9 55E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base A (26/0);  total (36/19)
Fisher's exact test for biased strand distribution p-value = 4.12e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508512‑3508811
                                                                                                                                                       |                                                                                                                                                     
ccACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAAt                                                                                                                                                       <  2:215512/151‑1 (MQ=255)
       gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGAAAAAATcacac                                                                                                                                                 >  1:860253/1‑150 (MQ=255)
           ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGT‑GCGCGGGAAAACTcaccatc                                                                                                                                               >  1:736600/1‑144 (MQ=255)
            tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGGCCGGGAAAAAttac                                                                                                                                                    >  1:791402/1‑141 (MQ=255)
             aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTa                                                                                                                                           <  1:346311/150‑1 (MQ=255)
              aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAg                                                                                                                                          >  1:259765/1‑150 (MQ=255)
                    gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCcac                                                                                                                                   <  1:437006/151‑1 (MQ=255)
                      gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACAt                                                                                                                                 <  1:429686/151‑1 (MQ=255)
                        gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAAATCACACACGGTAGCCCCcac                                                                                                                                   >  1:546771/1‑147 (MQ=255)
                                 gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                      <  1:32004/151‑1 (MQ=255)
                                         tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCAAACACGGGAGCCCCCACATCAGAAAGCGTGGGgcgccc                                                                                                              >  1:702503/1‑149 (MQ=255)
                                         tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCCGCGAAAAAAAAACACACCGGAGACCCCCCATCAAAAAGCgggggggg                                                                                                                 >  2:869994/1‑144 (MQ=255)
                                         tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGG‑GGGCGGAAAAAATCACACACGGGGGACCCCACATAAGAAAACGGGGGGCGca                                                                                                               >  1:28009/1‑150 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGCGGGAAAATTCACACACAGGAGCCCCCACATCAGAAAGAGTGGGGGGca                                                                                                               >  2:631039/1‑149 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGcac                                                                                                              >  2:400112/1‑150 (MQ=255)
                                          cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGca                                                                                                               <  2:364310/149‑1 (MQ=255)
                                                  aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGCGCACAcc                                                                                                           >  1:866877/1‑145 (MQ=255)
                                                  aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCAAACACAGGAGACCCCACATCAGAAAGAGGGGGGGGCGCAcc                                                                                                           >  2:558125/1‑145 (MQ=255)
                                                   gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                    <  2:572593/151‑1 (MQ=255)
                                                        tAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GGGCGGGAAAAATCACACACGGTAGGCCCCCCATCAGAAAGCGGGGGGCGCACAcc                                                                                                           >  1:93083/1‑139 (MQ=255)
                                                        tAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGG‑GCGCGGGAAAAATTACACCCCGGAGGCCCCCCATCCGAAAACGGGGTGCGCCCACCGGGCcgcg                                                                                                   >  2:784342/1‑147 (MQ=255)
                                                               tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGGGCGCACACCCAGCTGCGTAGCGatca                                                                                          >  1:171791/1‑149 (MQ=255)
                                                                  gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACAACCAGCTGCGGAGCGATCATccc                                                                                      >  2:352215/1‑150 (MQ=255)
                                                                     cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                                                                   <  1:689798/150‑1 (MQ=255)
                                                                       aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATTAAACACGGTAGCCCCCACATCAGAAAGCGGGGGGCGCCCACCCAGCTGCGTAGCGATCATcccct                                                                                    >  1:146054/1‑145 (MQ=255)
                                                                           tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGC‑GCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGAGTGGGGGGCACACCGAGCTGCGGAGCGGTCAACCCctgctggt                                                                              >  2:255574/1‑147 (MQ=255)
                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                         <  1:151576/150‑1 (MQ=255)
                                                                                gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAg                                                                       <  1:560079/151‑1 (MQ=255)
                                                                                    gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                                                                    <  2:449000/150‑1 (MQ=255)
                                                                                       tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAAAGGGGGGGGGACACCGAGCCGCGTAGCGATCATCCCTTGCTGATAAACCAGCAccg                                                                  >  2:678111/1‑147 (MQ=255)
                                                                                       tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAACACTgg                                                                 >  1:810865/1‑150 (MQ=255)
                                                                                          tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGG‑GCGGGGGAAAAATTACAAAAGGGAGCCCCCCCAACCAAAAGCGTGGTGGGGCCACCCCGCCGCGGAGCGGTCCTCCCCTTCTGGTAAAACCGCAcctgcccc                                                             >  2:73626/1‑144 (MQ=255)
                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGC‑GCGCGGGAAAAATCACACACGGGAGCCCCCCCCACAGAAAACGGGGGGCGCCCACCGAGCCGCGGAGCGATCATCCCCTGCTTATAAACCAacacccggacc                                                             >  2:81356/1‑143 (MQ=255)
                                                                                             aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGC‑GCGCGGGAAAAATCACACACGGTAACCCCCACATCAGAAAGCGGGGGGCGCACACCGCGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACTGGcccca                                                            >  2:760408/1‑149 (MQ=255)
                                                                                              aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGC‑GCGCGGGAAAATTTACACACGGTAGCCCCCCCACCAGAAAACGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGcaccac                                                           >  2:789132/1‑149 (MQ=255)
                                                                                               aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACCGTAGCCCCCCCATCAGAAAGCGCGGGGGGCACCCCCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGcaccccct                                                         >  1:49679/1‑146 (MQ=255)
                                                                                                  ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg                                                     <  1:313539/151‑1 (MQ=255)
                                                                                                   gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                    <  2:231037/151‑1 (MQ=255)
                                                                                                       gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAAGGTGGGGGGGACACCCAGCTGCGTAGCGGTCATCCCCTTCTGATAAAACAGCACCGGCACCACTTCGCggg                                                   >  2:189500/1‑148 (MQ=255)
                                                                                                         ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                              <  2:383633/151‑1 (MQ=255)
                                                                                                          gaaagaaaAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGTGGGCCCACCGCGCTGCGTAGGGGTCATCCCTTTCTGATAAAACAACACTGGCACCACTTCGCggg                                                   >  1:552476/1‑145 (MQ=255)
                                                                                                              gaaaagaGCGACCATTTTGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCCCCACCTCCCGGGGAtggggg                                           >  1:859404/1‑149 (MQ=255)
                                                                                                                     gCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                                  <  1:69351/151‑1 (MQ=255)
                                                                                                                     gCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                                  <  1:352223/151‑1 (MQ=255)
                                                                                                                       gACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAAt                                <  1:771746/151‑1 (MQ=255)
                                                                                                                        aCCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                               <  1:253852/151‑1 (MQ=255)
                                                                                                                            ttttGTGGCGGTGGC‑GCGCGGGAAAAATCACACACGGGAGACCCCCCCCCAGAAAACGGGGGGGGCACCACGAGCTGCGTCGCGATCATCCCCTGCTGATAAAACCGCACCGGCCCCACTTCGCCGGGATTGTGGCTGTGGCCAATgg                              >  2:478078/1‑148 (MQ=255)
                                                                                                                             tttGTGGCGGTTGG‑GCGCGGGGAAAATCACACAAGGTAGCCCCCACATCAGAAAGCGGGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACTTCGCGGGGATGGGGGCGGTGGCCCATTGGCgg                          >  2:879936/1‑151 (MQ=255)
                                                                                                                                 tggcggtggc‑gCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCg                      >  1:662030/1‑151 (MQ=255)
                                                                                                                                  ggcggtggc‑gCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt                     >  2:292527/1‑151 (MQ=255)
                                                                                                                                        ggc‑gCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt               <  2:866858/151‑1 (MQ=255)
                                                                                                                                            gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt            <  1:442572/151‑1 (MQ=255)
                                                                                                                                            gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt            >  2:50657/1‑151 (MQ=255)
                                                                                                                                               cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg         >  1:387234/1‑151 (MQ=255)
                                                                                                                                                  aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc       >  2:218004/1‑150 (MQ=255)
                                                                                                                                                    aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc    >  2:819862/1‑151 (MQ=255)
                                                                                                                                                      aTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAc  >  1:657320/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGC‑GCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCAC  >  NC_000913/3508512‑3508811

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