Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,040,433	C→A	A319D (GCC→GAC)	msrP →	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,433	0	C	A	96.8%	78.3 / NA	33	A319D (GCC→GAC)	msrP	protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (18/12);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCAC  >  NC_000913/2040191‑2040678                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      aCTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTgg                                                                                                                                                                                                                                               <  2:79384/251‑1 (MQ=255)          gacGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTTCTGTGTAATGGTTACGCCGACCAGGTGGCATCGCTGt                                                                                                                                                                                                                                      <  2:195885/249‑1 (MQ=255)             gTGGATATATGGCCTTAAAGGTATTAAATCGCTAGACAGCATTCCGCTGCCCCGCGAGCGTCCGCCAACCACCTGGTATCTGGGAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACATCTCCACCCATCGCTACTGTTTATTCGTTACGCCGCCCAGGTGGCATCGCTGTATc                                                                                                                                                                                                                                   >  2:280696/1‑251 (MQ=255)                    tatGGCTTTAAAGGGATTAAATCGATAGTCAGTATTACGCTGACACGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGACTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGCACGATTTATTGGTTCAGGCGGCCTCCTCGATGTACAGCGCCAGACAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCt                                                                                                                                                                                                                            >  2:183067/1‑251 (MQ=255)                      tGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTgg                                                                                                                                                                                                                          >  1:28098/1‑251 (MQ=255)                             aaaGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGc                                                                                                                                                                                                                    <  1:103185/250‑1 (MQ=255)                              aaGGGATTAAATCGATAGTCAGTATTAATCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTAGCATCGCTGCATCGCGGCCTGGATCTGCgg                                                                                                                                                                                                                  >  2:30708/1‑251 (MQ=255)                                  gATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCggg                                                                                                                                                                                                                 >  1:127673/1‑248 (MQ=255)                                           gATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGACTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCCCAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTaa                                                                                                                                                                                                      >  2:45173/1‑250 (MQ=255)                                                     aTTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTTGATTTGCGGGAGAATTTCTAAATGCGTCTGAc                                                                                                                                                                                           <  2:197559/251‑1 (MQ=255)                                                        aaGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGc                                                                                                                                                                                        >  1:220937/1‑251 (MQ=255)                                                         aGCTTACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTCTACGCCAACGGTAATCCGTATGTTGATCACCCGCGCGGGTGACAGGCTACCGAACGATGTATTGGTTCAGGCGGCAGCCTCGATGGACAGCGCCAGCGTACGTTACTGTTTAATGGTTACGCCGACCAGGTGGCATCTCTGTTTCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCa                                                                                                                                                                                       <  1:31895/251‑1 (MQ=255)                                                             gACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCaaa                                                                                                                                                                                     <  1:39279/249‑1 (MQ=255)                                                                ccGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGAGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTCTCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAgg                                                                                                                                                                                >  1:317488/1‑251 (MQ=255)                                                                           cGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTg                                                                                                                                                                     >  1:216299/1‑251 (MQ=255)                                                                                                             attcGGTTTTTACGCCAACGTTAATCCGTATGGTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACGGTTTAATGGTTACGCCGTCCAGGTGGCTTCGCTGTATCGTGGCCTGGATTTGCGTGAGAATTTCTAAATTCGTCTTACAGCAAAACAGGTGACATGGCTGAAAGTTTTCCTGCATCTTGCCGGGTTGTTGCCGt                                                                                                                                   <  2:123897/248‑1 (MQ=255)                                                                                                                     tttACGCCAACGTTAATCCGTATGTTGATCACCAGCGCTGGTCACAGGCTCCCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGTTTGTTGCCGTTTCTctgg                                                                                                                           >  2:215991/1‑251 (MQ=255)                                                                                                                             aaCGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGATGCCGTTTCTCTGGCTGGTCTg                                                                                                                   >  2:317028/1‑251 (MQ=255)                                                                                                                                       cGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGatcaat                                                                                                         <  2:79230/251‑1 (MQ=255)                                                                                                                                                                    gCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGaa                                                                            <  1:273875/251‑1 (MQ=255)                                                                                                                                                                               aTTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTGTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCTTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGc                                                                 <  2:28098/251‑1 (MQ=255)                                                                                                                                                                                tttATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTAATGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCATGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGc                                                                 >  2:290707/1‑250 (MQ=255)                                                                                                                                                                                           aGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTc                                                     >  1:190160/1‑251 (MQ=255)                                                                                                                                                                                                gcATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCGGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTGGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcact                                                <  1:280696/251‑1 (MQ=255)                                                                                                                                                                                                 cATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactg                                               >  1:125409/1‑251 (MQ=255)                                                                                                                                                                                                     cTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGCTTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcact                                           >  1:324701/1‑251 (MQ=255)                                                                                                                                                                                                                 cgcCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTTGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCt                               >  1:62706/1‑251 (MQ=255)                                                                                                                                                                                                                    cagccaACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGc                            >  2:130029/1‑251 (MQ=255)                                                                                                                                                                                                                    cagccaACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGGCGCACTGCACTGAAATTTTTGCTGGc                            <  2:216299/251‑1 (MQ=255)                                                                                                                                                                                                                     agccaACGCTACTGTTTAATGGCTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCg                           >  1:38598/1‑251 (MQ=255)                                                                                                                                                                                                                        caACGCTACTGGTTAATGGTGACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGTCATGGCTGAAAGTGTGCCTGCACCTTGCCGTATTGTTGCCGTTTCTCTGGCTGTTCTTGGCGATCAATCACGGTAGACTGGGTGCCGATCCGGTGAAATATAATCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAcc                        <  1:134444/251‑1 (MQ=255)                                                                                                                                                                                                                                       aTGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGACGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCct         >  1:319480/1‑251 (MQ=255)                                                                                                                                                                                                                                         ggTTACGCCGGCCAGGTGGCAGCGCGTTATCGTGGCCTGTATTTGCGGGAGAATTTCTAAATGCGGCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTGGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGGTGACTGGGTGCCGATCCGGTGAAAGATATTCACCATTTTACTGGTAGCACTGCACTGAAATTTTGGCTGGCTTCCTTGTTAATCGCCCctc        <  2:204626/250‑1 (MQ=255)                                                                                                                                                                                                                                              cgccgaCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCAc  <  1:130029/251‑1 (MQ=255)                                                                                                                                                                                                                                              cgccgaCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAAGCACCCCTCTGGCAc  <  2:220937/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      ACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCAC  >  NC_000913/2040191‑2040678

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