Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	58,435	G→C	10.2%	intergenic (‑346/+342)	PP_5424 ← / ← PP_0050	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	58,435	0	G	C	10.2%	173.0 / 10.1	59	intergenic (‑346/+342)	PP_5424/PP_0050	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base C (4/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

ACCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAG  >  NC_002947/58289‑58579                                                                                                                                                   |                                                                                                                                                  tccTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCt                                                                                                                                                 <  7:149154/148‑1 (MQ=255)  ccTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTg                                                                                                                                                <  3:265362/149‑1 (MQ=255)   cTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCt                                                                                                                                                 >  8:149154/1‑147 (MQ=255)            gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTg                                                                                                                                       >  5:13761/1‑148 (MQ=255)            gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTg                                                                                                                                       <  6:13761/148‑1 (MQ=255)              cggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAag                                                                                                                                    <  7:16258/149‑1 (MQ=255)                      gTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcg                                                                                                                            >  1:168530/1‑149 (MQ=255)                       tCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcgg                                                                                                                           <  2:150470/149‑1 (MQ=255)                             aTCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCa                                                                                                                     >  3:193531/1‑149 (MQ=255)                                       gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCg                                                                                                           >  2:167623/1‑149 (MQ=255)                                         ttgaCACGCGTGCCGATTTGGATAGAATTGCAATAGAGGGCTGCAACTGCGCATTTTTCGCGCCTGGTTACCACAATGACCACCAGCAATATTGGCAGAAGAGGAAGCTGATTCCTATGAAGAAGTGACTGCGGGCATAGCCAGACGGa                                                                                                         >  6:154172/1‑149 (MQ=255)                                            acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                     <  6:235327/118‑1 (MQ=255)                                            acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                     >  5:235327/1‑118 (MQ=255)                                            acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTc                                                                                                      >  1:27933/1‑149 (MQ=255)                                              acGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTg                                                                                                    <  8:152214/149‑1 (MQ=255)                                                gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                  >  6:273022/1‑149 (MQ=255)                                                    gCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCATTTTTGGCAGCAGAGGAAGCTGGTTGCTCTGAAGTAGTGCCGGCGGGCATTGCCAGGCTGATTCTGAAcc‑c                                                                                              >  2:265937/1‑148 (MQ=255)                                                            ggATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGc                                                                                      <  2:27933/148‑1 (MQ=255)                                                             gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGc                                                                                      >  5:72980/1‑147 (MQ=255)                                                             gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCa                                                                                    <  5:144089/149‑1 (MQ=255)                                                                 gCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGAc                                                                                 <  5:273022/148‑1 (MQ=255)                                                                  ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGAc                                                                                 >  8:227244/1‑147 (MQ=255)                                                                   cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                     >  6:250333/1‑93 (MQ=255)                                                                   cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                     <  5:250333/95‑3 (MQ=255)                                                                    ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                              >  6:274644/1‑148 (MQ=255)                                                                    ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                              >  8:87790/1‑148 (MQ=255)                                                                      gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAgg                                                                             <  2:168530/147‑1 (MQ=255)                                                                      gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                              <  7:87790/146‑1 (MQ=255)                                                                             aGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAg                                                                     >  3:104735/1‑148 (MQ=255)                                                                             aGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAg                                                                     <  5:22532/148‑1 (MQ=255)                                                                                  tGCAACTGCGCATTCGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                       <  1:54880/110‑1 (MQ=255)                                                                                  tGCAACTGCGCATTCGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                       >  2:54880/1‑110 (MQ=255)                                                                                      aCTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCg                                                           <  4:104735/149‑1 (MQ=255)                                                                                        tGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                              <  8:281484/128‑1 (MQ=255)                                                                                        tGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                              >  7:281484/1‑128 (MQ=255)                                                                                           gcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                               <  5:143613/61‑1 (MQ=255)                                                                                           gcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                               >  6:143613/1‑61 (MQ=255)                                                                                             aTTTGTCGCGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGCGCATTGCCAGTCGGATTCTGAACC‑AAAGGTCGCCACACGAGGAGCTAAAGGATAGCTCCGACTActc                                                    >  2:90571/1‑149 (MQ=255)                                                                                              tttGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                   <  4:193531/149‑1 (MQ=255)                                                                                                  tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                               >  4:336597/1‑149 (MQ=255)                                                                                                  tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                               <  6:72980/149‑1 (MQ=255)                                                                                                           ggTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTAcc                                       <  5:274644/148‑1 (MQ=255)                                                                                                             tACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGa                                     >  5:129162/1‑148 (MQ=255)                                                                                                                    aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATccc                                            <  3:340089/134‑1 (MQ=255)                                                                                                                    aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATccc                                            >  4:340089/1‑134 (MQ=255)                                                                                                                           cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcgg                      <  3:336597/149‑1 (MQ=255)                                                                                                                            ccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCCCGCgtcggt                     >  1:263307/1‑149 (MQ=255)                                                                                                                             cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                             <  1:322056/140‑1 (MQ=255)                                                                                                                             cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                             >  2:322056/1‑140 (MQ=255)                                                                                                                                      tacttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGc           >  8:173846/6‑149 (MQ=255)                                                                                                                                      ggCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGc           >  8:66877/1‑149 (MQ=255)                                                                                                                                        cttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGc                                                            <  3:221898/95‑1 (MQ=255)                                                                                                                                        cttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGc                                                            >  4:221898/4‑98 (MQ=255)                                                                                                                                         ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGCCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGAGACTACTCCTCATACCTTACCGAAAAGTCAAGCGTCGGTAGTAAACGGCCCt        >  5:97768/3‑149 (MQ=255)                                                                                                                                         agcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCTGAAGGTAGCCAGACTAGGAGCTAAAGGATAGCTGCGACTACAACGAATCCCTTAACGCAAATTAACGCGGCGGTCGTTAACGGccc         >  1:107067/1‑149 (MQ=255)                                                                                                                                         agcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt        <  6:129162/149‑1 (MQ=255)                                                                                                                                           cagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCttt      >  4:62937/1‑149 (MQ=255)                                                                                                                                              agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc     >  1:235294/1‑147 (MQ=255)                                                                                                                                               ggAAGCTGTTTACTCTGAATTATTGACTGCGGGCATTGCCAGTCGGATTCTGAACC‑CAATGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTATCCAg  <  5:154172/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ACCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAG  >  NC_002947/58289‑58579

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