Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	962,240	C→G	6.7%	E59Q (GAG→CAG)	selD ←	selenide, water dikinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	962,240	0	C	G	6.7%	194.5 / 3.8	60	E59Q (GAG→CAG)	selD	selenide, water dikinase
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (2/2);  total (28/32)
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.63e-01

ACCGAGGATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGG  >  NC_002947/962097‑962382                                                                                                                                                |                                                                                                                                                cccgctgatggccatcgcCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTcgtcgt                                                                                                                                            >  8:158045/7‑149 (MQ=255) aCCGAGGATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTcgtcgt                                                                                                                                            >  6:21480/1‑149 (MQ=255)        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)                                                        aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCtcg                                                                                                                                                <  1:229753/90‑1 (MQ=255)                                                        aTGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCtcg                                                                                                                                                >  2:229753/1‑90 (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                                   >  8:127572/1‑123 (MQ=255)                                                                                                                                   aCCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTg                                   <  7:127572/123‑1 (MQ=255)                                                                                                                                     cACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGgcgc       >  6:209017/1‑149 (MQ=255)                                                                                                                                     cACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGgcgc       >  5:50732/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)                                                                                                                                                |                                                                                                                                                ACCGAGGATGGCGATGGCCATCAGCGGGTCGCCGCCCATGGCATAGATGTCGCTGATGGCGTTGGTGGCGGCAATGCGGCCGAAGTCGTAAGGGTCATCGACGATCGGCATGAAGAAGTCGGTGGTCGAGACCACGCCGCGCTCGTCGTCCAGCGCGTACACTGCTGCATCGTCGCGCGATGCG‑TTGCCGACCCACAGTTTCGGGTCCAGGGCCTGGGTGCCGCTTTCGGCAAGGATCACGTCCAGCACCTTGGGGGAGATCTTGCAACCACAGCCGGCGCCATGG  >  NC_002947/962097‑962382

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