Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	354,288	0	T	G	50.9%	65.7 / 91.7	55	V528G (GTG→GGG)	prpE	propionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (6/21);  new base G (28/0);  total (34/21)
Fisher's exact test for biased strand distribution p-value = 3.52e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-08

CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354145‑354436                                                                                                                                                |                                                                                                                                                      cATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGaa                                                                                                                                                    >  1:519722/1‑147 (MQ=255)  aTCCGTGACGCTGACGGTTATCACTATATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATg                                                                                                                                               <  1:461062/151‑1 (MQ=255)    ccGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATgc                                                                                                                                              <  2:253991/150‑1 (MQ=255)        gACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑g                                                                                                                                         <  1:634549/150‑1 (MQ=255)              gACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGAAGCGCT‑GAAAgggg                                                                                                                                  >  1:7739/1‑150 (MQ=255)               aCGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GGAAggg                                                                                                                                   >  1:919889/1‑149 (MQ=255)                     aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGGAggtgg                                                                                                                            >  1:267657/1‑150 (MQ=255)                     aTCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGAGGGTGGGGGGAAAGAAGCGGC‑GAAAGGGCAggggg                                                                                                                            >  1:584629/1‑150 (MQ=255)                                   gggCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGGTTGt                                                                                                               >  1:92836/1‑149 (MQ=255)                                     gcgcACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCAtt                                                                                                           <  2:241764/151‑1 (MQ=255)                                        cACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCT‑GGAAGGGGAgggggcggggggggtttttatttc                                                                                                         >  2:429281/1‑137 (MQ=255)                                          atGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGTGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGCGGTTCGGGGGGAAGAAGGGCT‑GGAAGGGCGggggggggtgggggtttttat                                                                                                            >  1:398722/2‑135 (MQ=255)                                           tGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGTGGAAGAGGCGCT‑GACAGGGCAGGGGGGGGGGGGGGTTGTTCTTcc                                                                                                         >  2:242003/1‑147 (MQ=255)                                            gACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGTTGGGGGGAGAGAAGGGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGGCATTCCGaaat                                                                                                    >  1:693082/1‑150 (MQ=255)                                             aCGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAaga                                                                                                   <  2:853046/151‑1 (MQ=255)                                                    aaTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGa                                                                                             <  2:565495/150‑1 (MQ=255)                                                     aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACa                                                                                           <  1:8110/151‑1 (MQ=255)                                                     aTTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCCAAAGAGAGCg                                                                                              >  1:867735/1‑148 (MQ=255)                                                        aaCGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTc                                                                                        <  1:799128/151‑1 (MQ=255)                                                           gTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GGAAGGGGAGGGGGCGGGGGCGTTTGGTATTCCGCAAGAGAGGGGCAATCTgg                                                                                     >  2:416630/1‑151 (MQ=255)                                                            ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGGGCT‑GAAAGGGGAGGGGGGGGGGGGGGTTGTTATTTAGAAAGAGGGGGAAAGTCTgg                                                                                     >  2:835984/1‑150 (MQ=255)                                                            ttGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCT‑GGAAGGGCAGGGGGGGGGGGGGGTTGTCATTTCGAAAGAGAGCg                                                                                              >  2:200076/1‑141 (MQ=255)                                                               ccGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGGCGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGGGACAGTCTgg                                                                                     >  1:482322/1‑147 (MQ=255)                                                                 gggCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATc                                                                               <  2:591033/151‑1 (MQ=255)                                                                     aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAGGGGCT‑GGAAAGGGAGGGGGGGGGGGGGGTTGTCCTTTCGAAAGAGGGGGGCAGGCTGGgagagagg                                                                             >  2:28080/1‑145 (MQ=255)                                                                     aTCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGAAGCGCTGGAAAGGGC‑GGGGGGGGGGGGGGTTGTCATTTCGAAAAAGAGGGGCAGACTGGAAGAACg                                                                              >  2:644959/1‑148 (MQ=255)                                                                          cTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGGAAGAAGCGCT‑GGAAGGGCAGGGGGCGGTGGCGGTTGTTATTCCGGAAGAGAGCGACAGTCTGGAAGATCGTGGTgggg                                                                      >  2:797579/1‑151 (MQ=255)                                                                              ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTTGGGGGGAAGATGCGCT‑GAAAGGGGAGGTGGCGGGGGGGTTTGTAGTTCCGTAAGAGGACGACAATGTGGAAGaggggggtgttgcgc                                                                   >  1:688811/1‑136 (MQ=255)                                                                              ggACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCa                                                                  <  1:279238/151‑1 (MQ=255)                                                                               gACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGGGGGTGGGGGGGAAGGTGGCGT‑GGAAAGGCAGGGGGCGGGGGGGGTTGTTGATCCGCAAAAGCGGGACGGTTCTGAAGATCAGGATGTgggtcgc                                                                 >  2:126290/1‑146 (MQ=255)                                                                                 cgcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                               <  1:627216/151‑1 (MQ=255)                                                                                  gcgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCg                                                              <  2:793568/151‑1 (MQ=255)                                                                                   cgcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGc                                                             >  2:811895/1‑151 (MQ=255)                                                                                    gcgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGTGGGGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTc                                                               >  1:127183/1‑148 (MQ=255)                                                                                     cgAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGGGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGATGTGGCGCACTCg                                                              >  2:794021/1‑148 (MQ=255)                                                                                       agaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAaga                                                         <  1:142342/151‑1 (MQ=255)                                                                                        gaTTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCaa                                                           >  1:630606/1‑148 (MQ=255)                                                                                          ttGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAAATGCGCT‑GAAAGGGGAGGGGGGGGGGGGGGTTGTCATTCCGAAAGAGAGCGGAAGTCTGGAAGATCGTGGTGTGGGGCACTCGCaa                                                           >  2:72786/1‑146 (MQ=255)                                                                                            gAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAg                                                     <  2:545782/150‑1 (MQ=255)                                                                                                agagTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGAt                                                <  1:864945/151‑1 (MQ=255)                                                                                                     aTCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGCGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGGAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATgg                                           >  1:369014/1‑151 (MQ=255)                                                                                                         ccAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATgg                                           <  1:562631/147‑1 (MQ=255)                                                                                                          cAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtg                                      <  2:135160/151‑1 (MQ=255)                                                                                                              cATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGGGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCtggtgg                                  >  1:530034/1‑151 (MQ=255)                                                                                                              cATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGGGAAAGATGCGCT‑GAAAGGGCAGGGGGCGGGGGGGGTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGGGGTGTGGCGCACTCGCAAGAGAAGGCGATTATTGCGCtggggg                                  >  2:903623/1‑151 (MQ=255)                                                                                                               aTCCGGGCGTTGCCGAAGTGGCGGTGGGTGGGGGGAAAGAAGCGCT‑GAAAGGGGAGGGGGCGGTGGCGGTTGTTATTTCGTAAGAGAGCGACAGTCTGGGAGATCGTGATGGGGCGCCCTCGCAAGAGAAGGGGGTTATTGCGCtggtgg                                  >  2:567630/1‑150 (MQ=255)                                                                                                                  cGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATCCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACa                               <  2:512257/150‑1 (MQ=255)                                                                                                                       gTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGa                         <  1:372864/151‑1 (MQ=255)                                                                                                                          gCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTg                      <  2:609058/151‑1 (MQ=255)                                                                                                                            cGAAGTGGCGGTGGGTGGGGGGGAAGATGCGCT‑GAAAGGGCAGGGGGGGGGGGCGTTTGTCATTCCGGAAGAGAGGGACAGTCTGGGAGATCGGGATGTGGCGCACTCGCAAGAGAAGGGGATTATGGCGCTGGGGGGCAGCCAGATTgg                     >  2:637236/1‑150 (MQ=255)                                                                                                                                   gCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTg             >  2:675568/1‑151 (MQ=255)                                                                                                                                    cGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTgg            >  2:513067/1‑151 (MQ=255)                                                                                                                                      gtggtTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc          <  1:132850/151‑1 (MQ=255)                                                                                                                                       tggtTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGgcc          >  1:686309/1‑150 (MQ=255)                                                                                                                                              ggTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGgcg  >  2:81244/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                      CATCCGTGACGCTGACGGTTATCACTTTATTCTCGGGCGCACTGACGATGTAATTAACGTTGCCGGGCATCGGCTGGGGACGCGCGAGATTGAAGAGAGTATCTCCAGCCATCCGGGCGTTGCCGAAGTGGCGGTGGTTGGGGTGAAAGATGCGCT‑GAAAGGGCAGGTGGCGGTGGCGTTTGTCATTCCGAAAGAGAGCGACAGTCTGGAAGATCGTGATGTGGCGCACTCGCAAGAGAAGGCGATTATGGCGCTGGTGGACAGCCAGATTGGCAACTTTGGCCGCCCGGCG  >  NC_000913/354145‑354436

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