Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,336,499	A→G	8.9%	intergenic (‑43/+37)	PP_2052 ← / ← PP_5500	bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,336,499	0	A	G	8.9%	151.0 / 3.9	56	intergenic (‑43/+37)	PP_2052/PP_5500	bifunctional sugar‑phosphatase/mannitol‑1‑phosphate 5‑dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base G (3/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGAAGCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTGGCAACCGTTGCGGC  >  NC_002947/2336352‑2336637                                                                                                                                                     |                                                                                                                                           agaAGCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAAc                                                                                                                                           <  7:164963/149‑1 (MQ=255)  gaAGCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACa                                                                                                                                          <  4:42137/149‑1 (MQ=255)    aGCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAgg                                                                                                                                        <  1:163822/149‑1 (MQ=255)         gcttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCga                                                                                                                                     <  1:341023/147‑1 (MQ=255)          cttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGTTGCAGGCgagat                                                                                                                                  <  6:331161/149‑2 (MQ=255)          cttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGaa                                                                                                                                  >  1:175908/1‑149 (MQ=255)           ttgAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAAc                                                                                                                                 >  3:287227/1‑149 (MQ=255)             gAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACgg                                                                                                                               >  6:281400/1‑149 (MQ=255)               gtgtCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGgt                                                                                                                              <  7:108345/148‑1 (MQ=255)               gtgtCCTGAAGCGCACGCGCTCGGTGTCGAAGAGTGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACAGTGAAAGAAAAGCATTGGCTTTTGCCCCTTGTCGTTTCGTCA‑GGCATCTCGCCTGCAACATGCGAGAACGgtg                                                                                                                             >  2:260299/1‑149 (MQ=255)                    cTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGa                                                                                                                        <  1:178181/149‑1 (MQ=255)                     tGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAg                                                                                                                       >  3:56070/1‑149 (MQ=255)                       aaGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAg                                                                                                                       <  4:287227/147‑1 (MQ=255)                          cgcAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGTTGCAGGCgagatgc                                                                                                                                >  1:285900/1‑132 (MQ=255)                          cgcAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGTTGCAGGCgagatgc                                                                                                                                <  2:285900/135‑4 (MQ=255)                               gcgcTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTggg                                                                                                             <  6:156184/149‑1 (MQ=255)                               gcgcTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTggg                                                                                                             <  5:321967/149‑1 (MQ=255)                               gcgcTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGGCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTggg                                                                                                             <  8:87689/149‑1 (MQ=255)                                      gtgtCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTc                                                                                                      <  5:281400/149‑1 (MQ=255)                                           gAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTc                                                                                                 >  7:206530/1‑149 (MQ=255)                                               agGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCgggg                                                                                             <  8:206530/149‑1 (MQ=255)                                                   tGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATCGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGGAGTGGGGGCAGTCAGGGGGGGGGcg                                                                                          >  4:287915/1‑147 (MQ=255)                                                       atccatATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGACAACCGTGTGTGAGATGCAGTGGGCGCAGTCAGGTCGGGGGCACACGt                                                                                     >  1:311364/1‑149 (MQ=255)                                                        tccatATCGAAGATGGCGCTGAAAATCTGTATACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGtt                                                                                    >  1:28274/1‑149 (MQ=255)                                                          catATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTc                                                                                   >  7:56243/1‑148 (MQ=255)                                                          catATCGAAGATGGCGCGGAAAATCTGATTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTc                                                                                   <  8:218495/148‑1 (MQ=255)                                                           atatCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGa                                                                                 >  3:327366/1‑149 (MQ=255)                                                            tatCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGa                                                                                 >  1:234330/1‑148 (MQ=255)                                                             atCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGa                                                                                 >  5:46323/1‑147 (MQ=255)                                                              tCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAg                                                                              <  7:227882/149‑1 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCTTGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAgg                                                                             <  1:260299/148‑1 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCgggg                                                                                             <  7:7501/132‑1 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCgggg                                                                                             >  8:7501/1‑132 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAggg                                                                            <  3:258430/149‑1 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAggg                                                                            <  8:56243/149‑1 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAggg                                                                            >  7:288354/1‑149 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAggg                                                                            <  2:28274/149‑1 (MQ=255)                                                                gAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAgg                                                                             <  2:107336/148‑1 (MQ=255)                                                                          cTGAACATCTGTTGACCGTGAAAGAAAAGCATTCGCTGCTGCCCCTTGGCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTCCAGTGGGGGCAGTCAGGCCGGGGGCACACGTTCGACAGGGGGACGCTCAg                                                                  <  2:263473/149‑1 (MQ=255)                                                                                     tttACCGTGCACGAAAAGCATTGCCCGTTGCCGCTTGGCGTTGCGTCC‑GTCATCTCGCCTGCAACCGGCGAGAACGGTGTGTGCGCTGCAGTGGCGCCAGTCAGGTCGGGGGCCCACGTTCGACAGGGGGACGCTCAGGAGTGACGCag                                                       <  6:317456/149‑1 (MQ=255)                                                                                     gttACCGTGAAACCCCAGCCGTGCCTGTTGCCCCTTGGCGGTGCGCCC‑GGCAGCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCGGCCGGGGGGGCAGTCCGGGCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCag                                                       <  2:29942/148‑1 (MQ=255)                                                                                          cGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGt                                                   >  1:277725/1‑148 (MQ=255)                                                                                             gaaagaaaAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGCGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAgg                                               >  7:8792/1‑149 (MQ=255)                                                                                               aagaaaAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAggcg                                             >  6:70446/1‑149 (MQ=255)                                                                                                agaaaAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGa                                                                                 >  1:264068/1‑112 (MQ=255)                                                                                                agaaaAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGa                                                                                 <  2:264068/112‑1 (MQ=255)                                                                                                  aaaaGCATTGGCTGTTGCCCGTTGTCGTTGCGTCA‑GGCATCTCGCCTGGAGCAGGGGAGAAGGGTGTGTGTGCTGCAGTGGGGGCTGGCAGGTCGGGGGGACAGGTTGGGGAGGGGGAGGGTCTGGAGGGGGGGAGTGGTAGggcgggg                                          >  6:73274/1‑149 (MQ=255)                                                                                                    aaGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAggcggggcg                                        >  6:132275/1‑149 (MQ=255)                                                                                                            gCTGTTGCCCCTTGTCGTTGCGTCAGGGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCggg                                 >  1:177741/1‑149 (MQ=255)                                                                                                                      cTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGGAAGGCGGGGCGGAGCGGGGCGTGCagag                       >  3:78748/1‑148 (MQ=255)                                                                                                                         gTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGaggagg                   >  5:230989/1‑149 (MQ=255)                                                                                                                                       cgttctcGCCTGTTGCAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTGGCACCCGTTg      >  3:288304/4‑149 (MQ=255)                                                                                                                                           ctcGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAggg                                                                            >  1:316736/1‑75 (MQ=255)                                                                                                                                           ctcGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAggg                                                                            <  2:316736/75‑1 (MQ=255)                                                                                                                                            tcGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTGGCTACCGTTgcggc  >  8:221613/1‑148 (MQ=255)                                                                                                                                              gCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGCGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTGGCAACCGTTgcgt   <  8:8792/145‑2 (MQ=255)                                                                                                                                                     |                                                                                                                                           AGAAGCTTGCTTGAGTGTCCTGAAGCGCAGGCGCTCGGTGTCGAAGAGGGTGCCATCCATATCGAAGATGGCGCTGAAAATCTGTTTACCGTGAAAGAAAAGCATTGGCTGTTGCCCCTTGTCGTTGCGTCA‑GGCATCTCGCCTGCAACAGGCGAGAACGGTGTGTGAGCTGCAGTGGGGGCAGTCAGGTCGGGGGCACACGTTCGACAGGGGGACGCTCAGGAGTGACGCAGAGGTAAGGCGGGGCGGAGCGGGGCGTGCAGAGGAGGTTGGCAACCGTTGCGGC  >  NC_002947/2336352‑2336637

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