Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,835,219	0	A	C	41.1%	102.7 / 10.8	56	V736G (GTT→GGT)	hypF	carbamoyl phosphate phosphatase and [NiFe] hydrogenase maturation protein
Reads supporting (aligned to +/- strand):  ref base A (24/9);  new base C (0/23);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 9.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATACC  >  NC_000913/2835070‑2835364                                                                                                                                                        |                                                                                                                                                  ggATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAAc                                                                                                                                                  <  1:1056011/151‑1 (MQ=255)   aTGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAAccc                                                                                                                                                >  1:546300/1‑151 (MQ=255)      cGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCt                                                                                                                                             <  2:966639/151‑1 (MQ=255)        aCGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGc                                                                                                                                           <  2:612284/151‑1 (MQ=255)          gCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGccc                                                                                                                                         <  1:452058/151‑1 (MQ=255)            tGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGcccc                                                                                                                                        >  1:499205/1‑150 (MQ=255)            tGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGcccc                                                                                                                                        >  1:470061/1‑150 (MQ=255)             ggCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCa                                                                                                                                       <  1:961615/150‑1 (MQ=255)               cgcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCag                                                                                                                                      >  1:742194/1‑149 (MQ=255)               cgcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCag                                                                                                                                      >  1:582998/1‑149 (MQ=255)                gcgTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCagag                                                                                                                                    <  2:3628/150‑1 (MQ=255)                  gTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGAc                                                                                                                                  >  1:1062807/1‑150 (MQ=255)                                 cGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCAccc                                                                                                                   <  2:739995/150‑1 (MQ=255)                                  ggggACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACccgc                                                                                                                 >  1:949535/1‑151 (MQ=255)                                     gACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCgg                                                                                                              <  1:44450/151‑1 (MQ=255)                                        gCATGTGTAGGCCGGATAAGGCTTTTACCCCGCACCCGGCAAGGGTGTCCAAATGCAACCCTTTTTATCCTTTCTGGCCTCCCCCCGCT‑ACCCAC‑CGCCCCGCAGCAATACCCCCCTCCCCCAGAAACAACCCGCCATCACCCGCCGGTa                                                                                                            <  2:135526/150‑1 (MQ=255)                                        gCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTa                                                                                                            >  2:933168/1‑150 (MQ=255)                                                   cggaaaagggtttttACCCCGCACCCGAAATTGGTGCCAAAATGCAACACGTTTTACCCGTTCTGGCCTCCCCCCGCT‑AACCAA‑CGCCCCGCAAAAAAACCCCCCTCCCCCAGAAACAACCCGCCATCACCCGCCGGTAAACTCTGTGGaa                                                                                                <  2:1130574/139‑1 (MQ=255)                                                           gggttttCCGCCGCATCGGGCAATGGTTTCAAAATAAACCCCTTTTTATCCGTTCGGGCCTCCCCCCGCA‑AACCAC‑CGCGCCGCAGCAATACCCCCCTCCCCCAGAAAAAACCCCCCACCCCCCGCCGGTAAACTCTGTGGAAAGAGCAAt                                                                                        <  1:235571/147‑1 (MQ=255)                                                           ggCGTTTCCCCCGCACCCGGAAATGGTGCCAAAATGCAACCCGTTTTATCCGTTCTGGCCTTCCCCCCCA‑ACCCAC‑CGCCCCGCAGAAAAACCCCCCTCCCCAAGAGAAAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAAt                                                                                        <  1:117714/151‑1 (MQ=255)                                                                     ccgcTTCGGGAAAGGTTGTCAAAATCCACCCCTTTTTACCCTTTCTGGCCTTC‑CCCCCC‑AACCAACCGCCCCGCAGAAATACCCCCCGGCCCAAGAGAAAAACCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCa                                                                               <  2:373826/150‑1 (MQ=255)                                                                     ccgcATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATc                                                                                >  1:37449/1‑149 (MQ=255)                                                                     ccgcATCCGGCAATGGTGTCCAAAGCCAACCCTTTTTATCCGTTCTGGACTCCCCCCGCA‑AACCAC‑CGCCCCGCAGCAAAACCCCCCTGCCCAAGAAACAACCCCCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCa                                                                               <  1:563905/150‑1 (MQ=255)                                                                     ccccacccGGCAAGGGTGTCAAAATCCAACACTTTTTATCCGTTCGGGCCTCCCCCCGCA‑AACCAC‑CGCGCCGCAGAAATACCCCCCGCCCCAAGAGAAAACCCGCCACCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATc                                                                                <  2:267653/143‑1 (MQ=255)                                                                      cgcccccGGAAAGGGTGCCAAAATGCAACCCTTTTTACCCGTTCTGGACTCCCCCCGCA‑ACCCAC‑CGCCCCGCAGAAATACCCCCCCCCCCAAGAGACAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGc                                                                             <  1:839037/146‑1 (MQ=255)                                                                          tCCGCCAATGTTGACCAAAGGCACCTCGTGTTATCCTTTCGGGACTCCCCCCGTT‑ACCCAC‑GGCGCGGCAGAAATACCCCCCGGCCCAAAAGAAAAGCCGCCCCCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCgag                                                                          <  2:360613/150‑1 (MQ=255)                                                                          ccgggaaaTGGTGTCAAAATGAAACGCTTTTTATCCGTTCTGGCCTTCCCCCGCA‑ACCCAA‑CGCGCCGAAGAAAAACCCCCCTGCCCAAGAGAAAACCCCCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCgag                                                                          <  1:496033/144‑1 (MQ=255)                                                                            cGGCAATGGTGTCCAAATGCCAAACGTTTTATCCGTTCTGGCCTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGCAAATAGCAATGAGAATTCAg                                                                              >  2:1126207/1‑144 (MQ=255)                                                                             ggaattGGTGCCAAAATGCAACATTTTTTACCCTTCCTGGCCTCCCCCCCCA‑ACCCAC‑CGCCCCGCAGAAAAACCCCCCCCCCCAAGAAAAAACCCGCCACCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCgag                                                                          <  2:646425/142‑1 (MQ=255)                                                                             ggCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATa                                                                       >  1:472288/1‑150 (MQ=255)                                                                             ggCAATGGTGCCAAAATGCAACACGTTTTCCCCGTTCTGGCCTCCCCCCCAT‑ACCCAC‑CGCCCCGCAGAAATACCCCCCTGCCCAAGAAAAAACCCGCCACCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATa                                                                       <  2:891440/150‑1 (MQ=255)                                                                                aaTGGTGTCAAAATCCAACACTTTTTACCTTTTCTGGTTTCCCCCCCTT‑ACCCAC‑CGCGCCGCAAAAATACCCCCCTCCCCCAGAAACAACCCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAAt                                                                     <  1:1159808/149‑1 (MQ=255)                                                                                 gtGGTGCCCAAATCCACCCCTTTTTACCCGTTCTGGACTCCCCCCGCT‑ACCCAC‑CCCGCCGCAGAAAAACCCCCCGGCCCCAGAGAAAAACCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATgc                                                                   <  2:560493/149‑1 (MQ=255)                                                                                   ggTGTCCAAATGCAACACGTTTTATCCGTTCTAGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGCAAGAGCAATGTGAAATCAGCGAGATAATGCGcc                                                                >  1:902540/1‑151 (MQ=255)                                                                                     tgtCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGcc                                                                >  2:618796/1‑149 (MQ=255)                                                                                     tgtCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGcc                                                                >  2:237916/1‑149 (MQ=255)                                                                                      gtCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGcc                                                                >  2:991378/1‑148 (MQ=255)                                                                                        caaaaGGCAACCCTTTTTCCCCTTTCTGGCTTCCCCCCGCA‑ACCCAC‑CGCGCCGAAGAAAAACCCCCCTCCCCAAGAGACAACCCCCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCa                                                               <  1:215829/145‑1 (MQ=255)                                                                                          aatgacAACACTTTTTATCCGTCCGGGTCTTC‑CCCGCTAAACCAC‑CGCCCGGCACCAAAAACCCCCCCCCCCAGAGAAAACCCTCCTTCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACg                                                           <  1:532869/144‑1 (MQ=255)                                                                                          aaaTGCAACTCGTTTTACCCTTTCGGGACTCCCCCCGC‑‑AACCAACCGCCCCGCAGCAAAACCCCCCTCCCCAAGAAACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACg                                                           <  2:260166/149‑1 (MQ=255)                                                                                                aaccctttttcccctttCTGGCCTCCCCCCGCA‑ACCCAC‑CGCGCCGCAGAAATACCCCCCGCCCCAAAAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACgca                                                   <  2:565891/136‑1 (MQ=255)                                                                                                aaCACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACgc                                                    >  1:951728/1‑150 (MQ=255)                                                                                                  cgttttttACCCGTTCGGGACTCCCCCCGCA‑ACCCAC‑CGCGCCGCAGAAAAACCCCCCGCCCCCAGAAACAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCAc                                                      <  2:61177/142‑1 (MQ=255)                                                                                                      ttttATCCTTTCTGGACTTCCCCCGCT‑ACCCAC‑CGCGCGGCAGCATTACCCCCCTGCCCCAGAGAAAACCCCCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCaaa                                              <  2:364343/150‑1 (MQ=255)                                                                                                      ttttATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAAc                                             >  1:1025546/1‑151 (MQ=255)                                                                                                      ttttACCCTTTTTGGACTTCCCCCGCT‑ACCCAC‑CGCGCCGCACAAATACCCCCCGCCCCAAGAGAAAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACgcag                                                  <  1:547895/146‑1 (MQ=255)                                                                                                            cctttCTGGCCTCCCCCCCCT‑AACCAC‑CGCCCCGCAGCAATACCCCCCTCCCCCAGAGAAAACCCGCCATCCCCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGt                                          <  1:464632/145‑1 (MQ=255)                                                                                                               ttCTGGACTTCCCCCGCA‑AACCAA‑CGCGCCGCAACAATACCCCCCTCCCCCAGAAAAAACCCGCCACCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTAt                                       <  1:817223/148‑1 (MQ=255)                                                                                                                 cTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAAt                                   >  2:864753/1‑150 (MQ=255)                                                                                                                  tGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATa                                  >  2:157613/1‑150 (MQ=255)                                                                                                                           cccGCA‑ACCCAC‑CGCGCCGCAGAAATACCCCCCCCCCCCAGAGAAAACCCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCccg                            <  1:946117/147‑1 (MQ=255)                                                                                                                                t‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAAt                   >  1:15225/1‑151 (MQ=255)                                                                                                                                               gcagcaATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTg       >  1:425807/1‑150 (MQ=255)                                                                                                                                                 agcaATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGAt     >  1:586044/1‑150 (MQ=255)                                                                                                                                                   caATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATAcc  >  2:76939/1‑151 (MQ=255)                                                                                                                                                   caATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATAcc  >  2:605825/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                  GGATGCGACGCTGGCGCGTCTTATCCGACCTACGGGGACGCATGTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGTCCAAATGCAACACGTTTTATCCGTTCTGGACTTCACCCGCT‑AACCAA‑CGCGCCGCAGCAATAACCCCCTGCCCCAGAGACAAACCGCCATCACCCGCCGGTAAACTCTGTGGAAAGAGCAATGTGAAATCAGCGAGATAATGCGCCAGACGTGCACGCAGCAAACGGTTATGAATAACCCCGCCGCTAAATACCAGCGTAGTGATACC  >  NC_000913/2835070‑2835364

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