Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	962,247	C→G	6.6%	L56L (CTG→CTC)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	962,247	0	C	G	6.6%	197.9 / 4.1	61	L56L (CTG→CTC)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base C (26/31);  new base G (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACT  >  NC_002947/962104‑962389                                                                                                                                                |                                                                                                                                                atggcgatggcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAgcgc                                                                                                                                            <  2:154592/149‑1 (MQ=255)    gcgatggcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTa                                                                                                                                         <  8:263827/149‑1 (MQ=255)          gcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACActgc                                                                                                                                   >  1:325117/1‑149 (MQ=255)            cATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACActgctg                                                                                                                                 >  1:57371/1‑149 (MQ=255)             aTCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACtgctgc                                                                                                                                <  7:158045/149‑1 (MQ=255)                aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTGGACGACGAGcgcg                                                                                                                                           >  2:94298/1‑135 (MQ=255)                aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTGGACGACGAGcgcg                                                                                                                                           <  1:94298/135‑1 (MQ=255)                aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtc                                                                                                                             >  5:238955/1‑149 (MQ=255)                aGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtc                                                                                                                             >  2:92299/1‑149 (MQ=255)                 gCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtcg                                                                                                                            >  7:217602/1‑149 (MQ=255)                  cGGATCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCAtcgt                                                                                                                           >  3:154579/1‑149 (MQ=255)                       cgccgcCCATGGCATAGATGTCGCTGAGGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTcgcgc                                                                                                                      <  6:283931/149‑1 (MQ=255)                               aTGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGc                                                                                                                 <  7:41190/146‑1 (MQ=255)                               aTGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGc                                                                                                                 >  8:41190/1‑146 (MQ=255)                                    aTAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGa                                                                                                        >  5:378325/1‑149 (MQ=255)                                     tAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑tt                                                                                                             <  2:190147/143‑1 (MQ=255)                                     tAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGAATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑tt                                                                                                             >  1:190147/1‑143 (MQ=255)                                              cTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGttt                                                                                              >  6:247568/1‑149 (MQ=255)                                               tGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGAc                                                                                                       >  8:277188/1‑139 (MQ=255)                                               tGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGAc                                                                                                       <  7:277188/139‑1 (MQ=255)                                                 aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTc                                                                                             <  5:21480/147‑1 (MQ=255)                                                 aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCgg                                                                                           <  2:325117/149‑1 (MQ=255)                                                 aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCg                                                                                            <  2:57371/148‑1 (MQ=255)                                                     cGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTcc                                                                                       <  5:203284/149‑1 (MQ=255)                                                       ttGGTGGCGTCAATGCGTCCGAAGTCGTAAGGGTCATCGACGATCGGAATGAAGAAGTCGGTGGTCGAGACCACGCCGTGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAg                                                                                     <  7:251246/149‑1 (MQ=255)                                                         ggtggCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAgg                                                                                    <  6:378325/148‑1 (MQ=255)                                                           tggCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGcc                                                                                 >  8:216044/1‑149 (MQ=255)                                                               gtcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTggg                                                                             <  1:316876/147‑1 (MQ=255)                                                                gcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGt                                                                            >  8:33237/1‑149 (MQ=255)                                                                gcAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGt                                                                            >  7:269530/1‑149 (MQ=255)                                                                  aaTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTTCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGc                                                                          <  4:279224/149‑1 (MQ=255)                                                                   aTGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGcc                                                                         <  7:216044/149‑1 (MQ=255)                                                                     gCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGc                                                                       >  5:106566/1‑149 (MQ=255)                                                                              gTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCaa                                                              <  8:217602/149‑1 (MQ=255)                                                                                 gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCGTTTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAg                                                             <  1:92299/148‑1 (MQ=255)                                                                                 gTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGa                                                           >  5:327131/1‑149 (MQ=255)                                                                                  tAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCAGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGAt                                                          >  4:239482/1‑149 (MQ=255)                                                                                     gggTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCAc                                                       <  7:33237/149‑1 (MQ=255)                                                                                      ggTCATCGACGATCGGCATGAAGAAGTCGGTGGTTGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACg                                                      >  6:271955/1‑149 (MQ=255)                                                                                          aTCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTc                                                    >  6:316880/1‑147 (MQ=255)                                                                                             gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                              <  6:272040/118‑1 (MQ=255)                                                                                             gacgaTCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTgg                                                                              >  5:272040/1‑118 (MQ=255)                                                                                                  tCGGCATGAAGAAGTCGGTGGTTGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCTAGGATCACGTCCAGCACCTTg                                          <  5:271955/149‑1 (MQ=255)                                                                                                     gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       <  6:106566/149‑1 (MQ=255)                                                                                                     gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       <  4:154579/149‑1 (MQ=255)                                                                                                     gCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGGCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       <  6:238955/149‑1 (MQ=255)                                                                                                      cATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTggggg                                      <  6:327131/149‑1 (MQ=255)                                                                                                         gaagaaGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTgggg                                       >  7:65701/1‑145 (MQ=255)                                                                                                             aagTCGGTTGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCtt                               <  4:149623/149‑1 (MQ=255)                                                                                                                       tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACcaca                      <  1:173438/148‑1 (MQ=255)                                                                                                                       tCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACca                        >  2:173438/1‑146 (MQ=255)                                                                                                                            aCCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                          <  7:127572/123‑1 (MQ=255)                                                                                                                            aCCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                          >  8:127572/1‑123 (MQ=255)                                                                                                                              cACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGgcgc              >  5:50732/1‑149 (MQ=255)                                                                                                                              cACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGgcgc              >  6:209017/1‑149 (MQ=255)                                                                                                                               aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAgg                                                                                    >  5:158671/1‑78 (MQ=255)                                                                                                                               aCGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAgg                                                                                    <  6:158671/78‑1 (MQ=255)                                                                                                                                 tacgcgCTGGACGACGAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTc                                                    >  6:64107/3‑108 (MQ=255)                                                                                                                                 tacgcgCTGGACGACGAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTc                                                    <  5:64107/106‑1 (MQ=255)                                                                                                                                   cgcgCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATgg         <  5:247568/149‑1 (MQ=255)                                                                                                                                          gtcgtcCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCACCATGGCTGTACt  <  5:316880/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                ATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGGCTGTACT  >  NC_002947/962104‑962389

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