Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,564,768 C→T 9.7% intergenic (‑864/‑34) vgrG‑III ← / → glgA type VI secretion system protein/glycogen synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,564,7680CT 9.7% 162.3 / 8.3 62intergenic (‑864/‑34)vgrG‑III/glgAtype VI secretion system protein/glycogen synthase
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base T (3/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

CCCAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCACAA  >  NC_002947/4564619‑4564917
                                                                                                                                                     |                                                                                                                                                     
cccAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGcc                                                                                                                                                       >  2:927518/1‑150 (MQ=255)
  cAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAg                                                                                                                                                     <  1:1330846/150‑1 (MQ=255)
             accaACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtggc                                                                                                                                           <  2:1230697/149‑5 (MQ=255)
             accaACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtggc                                                                                                                                           >  1:1230697/1‑145 (MQ=255)
             accaACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                          >  2:1220607/1‑150 (MQ=255)
                     aacAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAggg                                                                                                                                  >  2:821555/1‑150 (MQ=255)
                               tCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGaa                                                                                                                        <  1:1158768/150‑1 (MQ=255)
                               tCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGaa                                                                                                                        <  1:1220607/150‑1 (MQ=255)
                               tCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGaa                                                                                                                        <  1:463213/150‑1 (MQ=255)
                                 gCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGt                                                                                                                      <  1:821555/150‑1 (MQ=255)
                                          cGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGAACCCGACCAAGCATAGGGACGGAGAGAAGTATGATAAGt                                                                                                             >  1:75687/1‑150 (MQ=255)
                                             aaCGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGcc                                                                                                          <  1:710498/150‑1 (MQ=255)
                                                gCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                >  1:1047017/1‑141 (MQ=255)
                                                gCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                <  2:1047017/141‑1 (MQ=255)
                                                  cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     <  2:1217495/150‑1 (MQ=255)
                                                  cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     <  1:927518/150‑1 (MQ=255)
                                                  cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     >  1:1189181/1‑150 (MQ=255)
                                                  cGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGt                                                                                                     <  1:40887/150‑1 (MQ=255)
                                                       gcgcCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                                <  2:1189181/150‑1 (MQ=255)
                                                          tcgtgcgcgCACGCACGGCACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACGTCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                                <  1:25561/146‑1 (MQ=255)
                                                          tcgtgcgcgCACGCACGGCACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACGTCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                                >  2:25561/2‑147 (MQ=255)
                                                          ccGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTc                                                                                             <  1:5632/150‑1 (MQ=255)
                                                                  gcgcgcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGAc                                                                                     <  2:714553/150‑1 (MQ=255)
                                                                      gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCCGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                >  1:1381028/1‑119 (MQ=255)
                                                                      gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATc                                                                                                                <  2:1381028/119‑1 (MQ=255)
                                                                        aCGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                      <  2:327360/79‑1 (MQ=255)
                                                                        aCGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                      >  1:327360/1‑79 (MQ=255)
                                                                             acaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                    <  1:1328333/92‑1 (MQ=255)
                                                                             acaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                    >  2:1328333/1‑92 (MQ=255)
                                                                             acaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAAcc                                                                          <  2:240974/150‑1 (MQ=255)
                                                                                     aTGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                                <  1:1173236/136‑1 (MQ=255)
                                                                                     aTGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                                >  2:1173236/1‑136 (MQ=255)
                                                                                       ggAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcgc                                                                <  2:62447/150‑1 (MQ=255)
                                                                                          aggCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAAGAGCGCCTCGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAg                                                             >  2:764506/3‑150 (MQ=255)
                                                                                             cACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcg                                                                 <  1:1251892/143‑1 (MQ=255)
                                                                                             cACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcg                                                                 >  2:1251892/1‑143 (MQ=255)
                                                                                             cACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACAAAGCAAAGGGACGGGGAGAAGTATGATCAGGACGGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCg                                                          >  2:691854/1‑150 (MQ=255)
                                                                                                        aTTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCg                                                                         <  1:240930/124‑1 (MQ=255)
                                                                                                        aTTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCg                                                                         >  2:240930/1‑124 (MQ=255)
                                                                                                         ttCAGCCAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGACCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGAt                                              >  2:1384244/1‑150 (MQ=255)
                                                                                                               aagcaGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGcc                                        >  1:1088588/1‑150 (MQ=255)
                                                                                                                 gcaGAGAGAACGAGCACGAAACATCGCCAAGGGTGCCAGCACCCGGCCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGACTCACGTAGACGCTTTCAAAGCGGATAACTGCGAGCCGCTCACCCCGGGTTTCGCCAc                                      >  2:202207/1‑150 (MQ=255)
                                                                                                                     agagaACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACa                                  >  1:113201/1‑150 (MQ=255)
                                                                                                                     agagaACCAGCACGAAACATCGCCACGGGGGCCAGAACCGGAGGATGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCAGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACa                                  >  2:275766/1‑150 (MQ=255)
                                                                                                                           ccAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCaaa                                                                            >  1:691769/1‑102 (MQ=255)
                                                                                                                           ccAGCACGAAACATCGCCAAGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCaaa                                                                            <  2:691769/102‑1 (MQ=255)
                                                                                                                            cAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGt                                                                                                             <  2:191558/68‑1 (MQ=255)
                                                                                                                            cAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGt                                                                                                             >  1:191558/1‑68 (MQ=255)
                                                                                                                             aGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGc                                                            <  1:265783/116‑1 (MQ=255)
                                                                                                                             aGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGc                                                            >  2:265783/1‑116 (MQ=255)
                                                                                                                               cACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                          <  2:1369874/36‑1 (MQ=255)
                                                                                                                               cACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                          >  1:1369874/1‑36 (MQ=255)
                                                                                                                                 cGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGc                                         <  2:623386/131‑1 (MQ=255)
                                                                                                                                 cGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGc                                         >  1:623386/1‑131 (MQ=255)
                                                                                                                                  gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCg                     >  1:306093/1‑150 (MQ=255)
                                                                                                                                  gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCg                     >  1:51808/1‑150 (MQ=255)
                                                                                                                                   aaaCATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCgg                    >  2:508647/1‑150 (MQ=255)
                                                                                                                                           ccACGGGTGCCAGCACCCGACCAAGCATAGGGACGGTGTGACGAAACTTGCGTTAGCCAGGCGATCCTCACTGATACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAgc            >  1:739121/1‑150 (MQ=255)
                                                                                                                                            cACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAgcg           >  2:617140/1‑150 (MQ=255)
                                                                                                                                             tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGa                                                                       <  1:1018340/88‑1 (MQ=255)
                                                                                                                                             tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGa                                                                       >  2:1018340/2‑89 (MQ=255)
                                                                                                                                             tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCg                                                 >  1:1269886/2‑111 (MQ=255)
                                                                                                                                             tcGGGTGCTGGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCg                                                 <  2:1269886/110‑1 (MQ=255)
                                                                                                                                                     cAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGGAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCACaa  >  1:116956/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCCAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCACAA  >  NC_002947/4564619‑4564917

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