Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	616,452	G→T	8.1%	noncoding (70/142 nt)	PP_mr11 ←	FMN regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	616,452	0	G	T	8.1%	170.6 / 2.2	56	noncoding (70/142 nt)	PP_mr11	FMN regulatory region
Reads supporting (aligned to +/- strand):  ref base G (27/24);  new base T (2/3);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGCAGCAG  >  NC_002947/616304‑616593                                                                                                                                                      |                                                                                                                                               aaCGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATg                                                                                                                                                <  1:316214/149‑1 (MQ=255) aaCGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATg                                                                                                                                                <  6:405338/149‑1 (MQ=255)    gCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTc                                                                                                                                             <  7:314608/149‑1 (MQ=255)           gTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACCATGTCggggg                                                                                                                                        <  2:276544/147‑2 (MQ=255)           gTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGgcgc                                                                                                                                      <  4:223420/149‑1 (MQ=255)                   tCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGc                                                                                                                              <  5:41999/149‑1 (MQ=255)                       tGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCAt                                                                                                                          <  8:42305/149‑1 (MQ=255)                           cGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCtt‑t                                                                                                                     <  6:360878/149‑1 (MQ=255)                                   ggCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTAc                                                                                                             <  7:221887/149‑1 (MQ=255)                                    gCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTAcc                                                                                                            >  3:356277/1‑149 (MQ=255)                                      gAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTAccgc                                                                                                          >  6:121601/1‑149 (MQ=255)                                       aaCAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAGCATCTGCTGACCCCCGACCATGTCgggggtca                                                                                                                                     <  1:247912/123‑5 (MQ=255)                                       aaCAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACCATGTCgggggtca                                                                                                                                     >  2:247912/1‑118 (MQ=255)                                            aaTGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTgg                                                                                                    >  3:42192/1‑149 (MQ=255)                                              tGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTgggg                                                                                                  <  4:136880/149‑1 (MQ=255)                                               gcgcgcgcAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTggg                                                                                                   >  7:13172/1‑147 (MQ=255)                                                   gcgcAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTa                                                                                             >  2:247896/1‑149 (MQ=255)                                                     gcAAGGGCGCCCGGATGGCCACTTCATGACGCCTTCTTTCATCCTGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATTGTCGGGTCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTaca                                                                                           <  5:139766/149‑1 (MQ=255)                                                      cAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTacac                                                                                          <  7:43506/149‑1 (MQ=255)                                                         gggCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACAcccc                                                                                       >  8:169148/1‑149 (MQ=255)                                                             gcCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGccc                                                                                   <  4:42192/149‑1 (MQ=255)                                                                cGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGa                                                                                <  7:91027/149‑1 (MQ=255)                                                                 ggATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGaa                                                                               <  7:354688/149‑1 (MQ=255)                                                                        cccTGCATGACGCCTTCTTTCATCCGGACTATGATCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgaccgac                                                                        >  5:149151/1‑149 (MQ=255)                                                                        cccTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAgaccgac                                                                        >  8:246043/1‑149 (MQ=255)                                                                        cccTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCAT‑TCGATTTACAGCCGGTGGGGACTTAAACACCGCCCTGAAgaccgac                                                                        >  8:33747/1‑149 (MQ=255)                                                                              aTGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCa                                                                  >  1:179577/1‑149 (MQ=255)                                                                                 tcGCCTTATTTCATCCTGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATa                                                                <  7:33747/147‑1 (MQ=255)                                                                                 aCGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATa                                                                <  8:160477/148‑1 (MQ=255)                                                                                 aCGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATa                                                                >  1:157755/1‑148 (MQ=255)                                                                                   gCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATATGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGAGTTACCGCCGGTGGGGACTTACACGCCGCCCTGAAGCCGGACCCAGCATACAg                                                             >  2:384836/1‑149 (MQ=255)                                                                                     cttcttTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGca                                                           >  1:248881/1‑149 (MQ=255)                                                                                        cttTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGcac                                                          <  2:179577/147‑1 (MQ=255)                                                                                          ttCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTg                                                      >  7:23585/1‑149 (MQ=255)                                                                                             aTCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGCCCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCCTACAGCACAGTGCCa                                                   >  2:393572/1‑149 (MQ=255)                                                                                               ccGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCAcc                                                 >  6:225664/1‑149 (MQ=255)                                                                                                  gACTATGATCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCa                                              <  6:149151/149‑1 (MQ=255)                                                                                                            gTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAAccc                                    <  7:169148/149‑1 (MQ=255)                                                                                                             tCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTggg                                                                                                   >  2:357565/1‑85 (MQ=255)                                                                                                             tCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTggg                                                                                                   <  1:357565/85‑1 (MQ=255)                                                                                                              cGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCcgc                                  <  4:356277/149‑1 (MQ=255)                                                                                                                 ctctGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGt                               >  1:228554/1‑149 (MQ=255)                                                                                                                      gAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctac                          <  8:13172/149‑1 (MQ=255)                                                                                                                      gAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTcta                           >  2:78647/1‑148 (MQ=255)                                                                                                                      gAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTcta                           <  1:78647/148‑1 (MQ=255)                                                                                                                       aGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctact                         <  5:121601/149‑1 (MQ=255)                                                                                                                        gTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctacta                        >  2:37991/1‑149 (MQ=255)                                                                                                                            cACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACCGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCtt                    >  5:208036/1‑149 (MQ=255)                                                                                                                             aCCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACAccc                                                                                        >  6:111582/1‑80 (MQ=255)                                                                                                                             aCCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACAccc                                                                                        <  5:111582/80‑1 (MQ=255)                                                                                                                              ccAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACGCCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCg                  >  4:389455/1‑149 (MQ=255)                                                                                                                                aG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGc                >  7:179216/1‑149 (MQ=255)                                                                                                                                             ccccGACATTGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGTCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcag     >  6:346298/1‑148 (MQ=255)                                                                                                                                             ccccGACATGGTCGGGGCGCTCGCGGGCTCATCTTGTCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctac                          >  1:265942/1‑128 (MQ=255)                                                                                                                                             ccccGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTctac                          <  2:265942/127‑1 (MQ=255)                                                                                                                                              cccGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTgcagca   <  2:228554/149‑1 (MQ=255)                                                                                                                                               ccGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGcagcag  >  5:369131/1‑149 (MQ=255)                                                                                                                                               ccGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGAATTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGcagcag  >  4:386304/1‑149 (MQ=255)                                                                                                                                                      ggTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGt                                                                                                      >  5:24234/1‑42 (MQ=255)                                                                                                                                                      ggTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGt                                                                                                      <  6:24234/42‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               AACGCTCCTCGTGAATGATCGATGAACGTCTCGGGGCGAACAAAATGCGCGCGCAAGGGCGCCCGGATGGCCCCTGCATGACGCCTTCTTTCATCCGGACTATGACCGTCGGCTCTGGAGTTGCACCAG‑ATCTGCTGACCCCCGACATGGTCGGGGCGCTCGCGGGCTCATCTT‑TCGATTTACCGCCGGTGGGGACTTACACCCCGCCCTGAAGACCGACCCAGCATACAGCACAGTGCCACCCCAGTGGCAACCCGCGGTTCTACTACCTTCGGCGAGACTTGCAGCAG  >  NC_002947/616304‑616593

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