Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,977,949	G→C	10.9%	P217P (CCC→CCG)	flgH ←	flagellar L‑ring protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,977,949	0	G	C	10.9%	154.8 / 10.6	55	P217P (CCC→CCG)	flgH	flagellar L‑ring protein
Reads supporting (aligned to +/- strand):  ref base G (23/26);  new base C (3/3);  total (26/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01

GCAATGTCCTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGG  >  NC_002947/4977806‑4978097                                                                                                                                                   |                                                                                                                                                     gCAATGTCCTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGggc                                                                                                                                                   <  4:72893/147‑1 (MQ=255)    aTGTCCTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGggctggct                                                                                                                                              <  8:52787/149‑1 (MQ=255)       tCCTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGctggctgg                                                                                                                                            >  3:219078/1‑148 (MQ=255)         cTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCAGCCCGGgtggctggac                                                                                                                                          <  8:253086/148‑3 (MQ=255)         cTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCAGCCCGGgtggctggac                                                                                                                                          <  2:71897/148‑3 (MQ=255)          ttCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATc                                                                                                                                        >  4:58251/1‑149 (MQ=255)             aGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGc                                                                                                                                     <  4:281437/149‑1 (MQ=255)                             cAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCg                                                                                                                     >  4:147805/1‑149 (MQ=255)                                cGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGa                                                                                                                   <  6:150260/148‑1 (MQ=255)                                 gAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAg                                                                                                                  >  6:276482/1‑148 (MQ=255)                                 gAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGt                                                                                                                 <  5:276482/149‑1 (MQ=255)                                  aaGTCGCAGGACAAAACCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCAACCGTACTCAGAAAGGCCAAACCTTCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGa                                                                                                                   <  1:245766/149‑1 (MQ=255)                                       gcAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGa                                                                                                           >  1:246782/1‑149 (MQ=255)                                          ggACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATgc                                                                                                        >  1:35403/1‑149 (MQ=255)                                              aaaaGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAg                                                                                                               <  8:212581/138‑1 (MQ=255)                                              aaaaGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAg                                                                                                               >  7:212581/1‑138 (MQ=255)                                                          gtgGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCg                                                                                                                     >  1:297077/1‑120 (MQ=255)                                                          gtgGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCg                                                                                                                     <  2:297077/120‑1 (MQ=255)                                                          gtgGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACg                                                                                        >  6:180144/1‑149 (MQ=255)                                                          gtgGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACg                                                                                        >  4:3735/1‑149 (MQ=255)                                                              cAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTg                                                                                    <  3:239018/149‑1 (MQ=255)                                                                 tCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCCCCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGAc                                                                                 >  5:45225/1‑149 (MQ=255)                                                                   aGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACgg                                                                               >  8:24799/1‑149 (MQ=255)                                                                        ccTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAATGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACgg                                                                          >  2:68490/1‑149 (MQ=255)                                                                        ccTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACg                                                                           <  3:147805/148‑1 (MQ=255)                                                                            aCGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGtgtt                                                                      <  3:163111/149‑1 (MQ=255)                                                                                  aaCATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGt                                                                <  5:180144/149‑1 (MQ=255)                                                                                  aaCATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGt                                                                <  3:3735/149‑1 (MQ=255)                                                                                    cATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTgg                                                              >  8:217446/1‑149 (MQ=255)                                                                                     aTGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGc                                                             <  4:219078/149‑1 (MQ=255)                                                                                         tCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGata                                                         >  3:238095/1‑149 (MQ=255)                                                                                           aTCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGatat                                                        <  3:58251/148‑1 (MQ=255)                                                                                              cGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGt                                                     <  7:217446/148‑1 (MQ=255)                                                                                                   tCAGAGAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATc                                                       <  4:26589/141‑1 (MQ=255)                                                                                                   tCAGAGAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATc                                                       >  3:26589/1‑141 (MQ=255)                                                                                                   tCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAgcgc                                               >  5:191892/1‑149 (MQ=255)                                                                                                    cAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAgcgc                                               >  2:90289/1‑148 (MQ=255)                                                                                                        aaGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAAcc                                          <  2:35403/149‑1 (MQ=255)                                                                                                                      gggCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAAt                               <  7:24799/149‑1 (MQ=255)                                                                                                                        gCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATgc                             >  7:182183/1‑149 (MQ=255)                                                                                                                        gCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATgc                             <  2:246782/149‑1 (MQ=255)                                                                                                                         cTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATgcg                            >  1:82714/1‑149 (MQ=255)                                                                                                                               aagaagCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGGGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGGGTTGTCGGTGGCGATATCGTCAGTGCGAACCATGCCGGCCATGCGCACCAg                      >  3:184828/1‑149 (MQ=255)                                                                                                                                 gaagCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGt                                                     >  6:213485/1‑116 (MQ=255)                                                                                                                                 gaagCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGt                                                     <  5:213485/116‑1 (MQ=255)                                                                                                                                    cgatgccagccagCCCGGGTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCtt                 >  2:113822/6‑149 (MQ=255)                                                                                                                                     cGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTc                <  1:90289/149‑1 (MQ=255)                                                                                                                                       tgccagccagCCCGGGTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGc              >  3:164862/3‑149 (MQ=255)                                                                                                                                       gTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGc              >  6:312151/1‑149 (MQ=255)                                                                                                                                         ccAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCGCGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCg            <  6:45225/149‑1 (MQ=255)                                                                                                                                          cagccagcCCGGGTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGctc                   <  7:14273/141‑1 (MQ=255)                                                                                                                                          cagccagcCCGGGTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGctc                   >  8:14273/1‑141 (MQ=255)                                                                                                                                          cAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCg            <  6:191892/148‑1 (MQ=255)                                                                                                                                             ccACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGgtgt        <  8:182183/149‑1 (MQ=255)                                                                                                                                                cccGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCa     <  4:238095/149‑1 (MQ=255)                                                                                                                                                   gggCTGGCTGGCATCGGCGAATGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAggg  <  1:68490/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCAATGTCCTTCAGACGCTCGGCCTGGGCAGCGAAGGCGCAGGACAAAAGCAGGGTTGTGGCAATCAGCTGCCTCACGTTGAACATGGTCATCCGTACTCAGAAAGGCCAAA‑‑‑GCGGGCTGATGAAGAAGCGGTCCAGCCACCCGGGCTGGCTGGCATCGGCGAACGAACCGGTACCGGAGTAGGTGATGCGCGCGTCGGCCACACGGGTGGACGGCACGGTGTTGTCGGTGGCGATATCGTCAGCGCGAACCATGCCGGCAATGCGCACCAGCTCTTCGCCGGTGTTCAGGG  >  NC_002947/4977806‑4978097

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