Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	858,322	C→G	6.1%	A168P (GCA→CCA)	phr ←	deoxyribodipyrimidine photolyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	858,322	0	C	G	6.1%	208.6 / 3.7	65	A168P (GCA→CCA)	phr	deoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base G (2/2);  total (31/34)
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.25e-01

AGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGG  >  NC_002947/858174‑858462                                                                                                                                                       |                                                                                                                                             aGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGc                                                                                                                                              >  1:293685/1‑149 (MQ=255)    aaGCGTGACAGGCGGGACTGTGCCTCAGCGTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAg                                                                                                                                           <  6:369889/149‑1 (MQ=255)        gTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGGAGTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAgggc                                                                                                                                        >  1:122705/1‑149 (MQ=255)        gTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAggga                                                                                                                                        >  6:158200/1‑147 (MQ=255)          gACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAgggcggg                                                                                                                                     >  7:353922/1‑149 (MQ=255)           aCAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAgggcgggc                                                                                                                                    <  1:62884/149‑1 (MQ=255)            cAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCa                                                                                                                                   <  3:175453/149‑1 (MQ=255)            atctcGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTACCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCCTGGCAcggcgggtg                                                                                                                                   <  7:198792/145‑3 (MQ=255)             aGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCTCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAg                                                                                                                                  >  2:271280/1‑149 (MQ=255)              ggCGGGACTGTGCCTCAGCTTCACCGTCTGGCCAGTGGTCGCGCTATGCGGTCGTGG‑TTTTTCCGAACGCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGTGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAgg                                                                                                                                 <  7:286779/149‑1 (MQ=255)               gCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGc                                                                                                                                >  3:295161/1‑149 (MQ=255)                 gggACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTACCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCCTGGCAcggcgggtg                                                                                                                                   >  8:198792/1‑142 (MQ=255)                    aCTGTTCCTCATCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCgg                                                                                                                           <  8:180677/149‑1 (MQ=255)                      tgtgCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGACCCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTg                                                                                                                         >  8:218178/1‑149 (MQ=255)                           cTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGTAGGCTGCGGTGCAGAt                                                                                                                    <  4:345787/149‑1 (MQ=255)                            tCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATg                                                                                                                   >  8:245774/1‑149 (MQ=255)                                      cGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGctg                                                                                                         <  2:293685/149‑3 (MQ=255)                                      cGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAg                                                                                                         <  5:302789/149‑1 (MQ=255)                                           ggCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGACCCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGtctt                                                                                                    <  8:333763/149‑1 (MQ=255)                                            gCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGtcttc                                                                                                   <  4:151237/149‑1 (MQ=255)                                            gCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGtcttc                                                                                                   <  4:298252/149‑1 (MQ=255)                                             ccAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGgattct                                                                                                  >  2:412848/1‑149 (MQ=255)                                                 tGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGaa                                                                                              >  5:322083/1‑149 (MQ=255)                                                     cgcgCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAg                                                                                                             >  6:55169/1‑130 (MQ=255)                                                     cgcgCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAg                                                                                                             <  5:55169/130‑1 (MQ=255)                                                      gcgcAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGc                                                                                         >  7:386018/1‑149 (MQ=255)                                                           aTGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAac                                                                                    <  4:295161/149‑1 (MQ=255)                                                              cGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCGCCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAacac                                                                                  <  4:171210/148‑1 (MQ=255)                                                              cGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACAcc                                                                                 <  1:271280/149‑1 (MQ=255)                                                                   gTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCGCCCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGt                                                                                                       <  1:85173/123‑1 (MQ=255)                                                                   gTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGt                                                                                                       >  2:85173/1‑122 (MQ=255)                                                                   gTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGaa                                                                            <  6:322083/149‑1 (MQ=255)                                                                           ttCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATc                                                                      <  7:218178/148‑1 (MQ=255)                                                                             ccGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGc                                                                    >  7:390369/1‑148 (MQ=255)                                                                                    gatatACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGTGGCTGGCGCTTCACCAG‑CCGTGCCAGGGCGGGCATGCTGCGGTGCAGATGCTCAATGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTTCga                                                            <  4:261214/145‑1 (MQ=255)                                                                                     ctctACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgag                                                           >  8:135295/1‑149 (MQ=255)                                                                                         aCATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGGGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGGAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagtga                                                        <  8:353922/148‑1 (MQ=255)                                                                                          cATGTTGCGTGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCga                                                            >  3:261214/1‑143 (MQ=255)                                                                                           aTGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                          >  7:387026/1‑144 (MQ=255)                                                                                           aTGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                          <  8:387026/144‑1 (MQ=255)                                                                                              ttGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACgg                                                  >  2:66606/1‑149 (MQ=255)                                                                                              ttGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACgg                                                  >  5:163093/1‑149 (MQ=255)                                                                                                gCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTg                                                >  4:245920/1‑149 (MQ=255)                                                                                                    gATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcg                                            >  1:212321/1‑149 (MQ=255)                                                                                                     aTTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcgg                                           >  8:169930/1‑149 (MQ=255)                                                                                                        gggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggac                                         <  2:212321/148‑1 (MQ=255)                                                                                                         ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCgg                                                                                                                           <  3:84710/65‑1 (MQ=255)                                                                                                         ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                          >  2:406131/1‑130 (MQ=255)                                                                                                         ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCgagt                                                          <  1:406131/130‑1 (MQ=255)                                                                                                         ggTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggacg                                        >  2:25547/1‑148 (MQ=255)                                                                                                              cTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                  <  6:163093/149‑1 (MQ=255)                                                                                                              cTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                  <  1:25547/149‑1 (MQ=255)                                                                                                              cTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                  >  8:245892/1‑149 (MQ=255)                                                                                                               tGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                  >  6:56278/1‑148 (MQ=255)                                                                                                                gCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGgaagaa                                  >  1:126216/1‑147 (MQ=255)                                                                                                                  cGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGc                              <  3:245920/149‑1 (MQ=255)                                                                                                                     tGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTgg                           <  7:169930/149‑1 (MQ=255)                                                                                                                                   tGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGt                                                                          >  6:215169/1‑88 (MQ=255)                                                                                                                                   tGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGt                                                                          <  5:215169/88‑1 (MQ=255)                                                                                                                                       gcTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGcc         <  7:245774/149‑1 (MQ=255)                                                                                                                                         ttCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCTCTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTgg                           <  2:276134/129‑1 (MQ=255)                                                                                                                                         ttCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCTCTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTgg                           >  1:276134/1‑129 (MQ=255)                                                                                                                                          ctgcccG‑CCCTGGCACGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggac                                         <  6:268633/111‑1 (MQ=255)                                                                                                                                          ctgcccG‑CCCTGGCACGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCcggac                                         >  5:268633/4‑114 (MQ=255)                                                                                                                                            aCCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGAt    >  7:230003/1‑149 (MQ=255)                                                                                                                                              ccg‑ccgTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATgg  <  7:135295/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                             AGGAAGCGTGACAGGCGGGACTGTGCCTCAGCTTCACCGGCTGGCCAGTGGTCGCGCAATGCGGGCGTGG‑TTTTTCCGAACCCCTCTACATGTTGCGGGATTGGGTCGCTGCCGATGTTCAGCGGGGCCTGGCGCTTCACCCG‑CCGTGCCAGGGCGGGCAGGCTGCGGTGCAGATGCTCAAGGCAGGTCTTCTTGAACTGGCTGAACACCTGGAAGTAATCGCCACTGCGAGTGAGTACGGTGCCCGGACGGAAGAACAGCTGGTCGAGGTGGCTGTGCGCCTGGATGG  >  NC_002947/858174‑858462

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