Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	857,638	C→A	7.4%	A396S (GCG→TCG)	phr ←	deoxyribodipyrimidine photolyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	857,638	0	C	A	7.4%	206.8 / 3.4	66	A396S (GCG→TCG)	phr	deoxyribodipyrimidine photolyase
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base A (2/3);  total (32/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 = 9.98e-01

AAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAG  >  NC_002947/857494‑857786                                                                                                                                                 |                                                                                                                                                      aaaGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCg                                                                                                                                                   <  2:225029/149‑1 (MQ=255)     gTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGATGAAGATACGGAAATAGGGCACCGCGTCGGTGc                                                                                                                                               <  4:106262/149‑1 (MQ=255)        cGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGACGCGGTGccc                                                                                                                                             <  6:214664/148‑2 (MQ=255)         gTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGt                                                                                                                                           <  6:118129/149‑1 (MQ=255)                 aaTGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCa                                                                                                                                   <  2:213285/148‑1 (MQ=255)                   tGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTc                                                                                                                                >  7:42240/1‑149 (MQ=255)                       ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACt                                                                                                                            >  5:132175/1‑149 (MQ=255)                       ttGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACt                                                                                                                            <  2:163874/149‑1 (MQ=255)                         gatggatgGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTg                                                                                                                           <  6:132175/148‑1 (MQ=255)                             gatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑g                                                                                                                                     >  4:135725/1‑134 (MQ=255)                             gatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑g                                                                                                                                     <  3:135725/134‑1 (MQ=255)                             gatggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAAc                                                                                                                      >  1:299771/1‑149 (MQ=255)                              atggCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAAcc                                                                                                                     >  7:55667/1‑149 (MQ=255)                                  cTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGAGACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCg                                                                                                                 >  6:294136/1‑149 (MQ=255)                                     tttCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTc                                                                                                                                                    >  8:46340/1‑112 (MQ=255)                                     tttCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTc                                                                                                                                                    <  7:46340/112‑1 (MQ=255)                                     tttCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCgttg                                                                                                              >  2:47673/1‑149 (MQ=255)                                     tttCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCgtt                                                                                                               <  2:299771/148‑1 (MQ=255)                                          tCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCg                                                                                                                 >  7:197489/1‑141 (MQ=255)                                          tCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCg                                                                                                                 <  8:197489/141‑1 (MQ=255)                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGaa                                                                                                                                       >  3:82209/1‑112 (MQ=255)                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGaa                                                                                                                                       <  4:82209/112‑1 (MQ=255)                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTc                                                                                                                                <  7:19625/118‑1 (MQ=255)                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTc                                                                                                                                >  8:19625/1‑118 (MQ=255)                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAgg                                                                                                  <  3:258809/148‑1 (MQ=255)                                                  cTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAg                                                                                                   <  8:42240/147‑1 (MQ=255)                                                    gTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGACGCGGTGccctatt                                                                                                                                         <  7:148792/108‑6 (MQ=255)                                                    gTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGACGCGGTGccctatt                                                                                                                                         >  8:148792/1‑103 (MQ=255)                                                       gTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCg                                                                                            >  5:303240/1‑149 (MQ=255)                                                        ttCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCg                                                                                                                                                   <  5:65144/94‑1 (MQ=255)                                                        ttCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCg                                                                                                                                                   >  6:65144/1‑94 (MQ=255)                                                           cGGTACCCAGTGGCGAATGAACGGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTa                                                                                        <  6:303240/149‑1 (MQ=255)                                                           cGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTa                                                                                        <  8:55667/149‑1 (MQ=255)                                                            ggTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAACGCACTCCACTGCCAACCCCCGTTGTTGGGGACCTGgccgccaacaa                                                                                        <  1:259652/149‑6 (MQ=255)                                                             gTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTAt                                                                                       >  6:180367/1‑148 (MQ=255)                                                                  cAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTg                                                                                 >  8:120104/1‑149 (MQ=255)                                                                         gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCa                                                                                                                              >  8:15476/1‑97 (MQ=255)                                                                         gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCa                                                                                                                              <  7:15476/97‑1 (MQ=255)                                                                         gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCAt                                                                           >  2:118731/1‑148 (MQ=255)                                                                         gAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCAt                                                                           >  7:111453/1‑148 (MQ=255)                                                                             gAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTgg                                                                                                                                         <  8:281769/83‑1 (MQ=255)                                                                             gAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTgg                                                                                                                                         >  7:281769/1‑83 (MQ=255)                                                                               aCCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGcc                                                                                             >  4:86064/1‑124 (MQ=255)                                                                               aCCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGc                                                                    >  2:326044/1‑149 (MQ=255)                                                                                ccGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCa                                                                                                                              >  5:317099/1‑90 (MQ=255)                                                                                ccGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCa                                                                                                                              <  6:317099/90‑1 (MQ=255)                                                                                   aCCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGTTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCt                                                                <  1:47673/149‑1 (MQ=255)                                                                                     cttgtgGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCg                                                              >  6:121163/3‑149 (MQ=255)                                                                                        gtgGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGccc                                                           >  1:221266/1‑149 (MQ=255)                                                                                           gATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTg                                                        >  8:275869/1‑149 (MQ=255)                                                                                            aTCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTgc                                                       <  7:275869/149‑1 (MQ=255)                                                                                              cGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTgcgc                                                     >  1:253069/1‑149 (MQ=255)                                                                                                aaGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTgcgc                                                     >  3:262048/1‑147 (MQ=255)                                                                                                    gcTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCa                                               <  5:180367/149‑1 (MQ=255)                                                                                                     cTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCaa                                              <  5:121163/149‑1 (MQ=255)                                                                                                          aCTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATcagc                                         <  3:136215/149‑1 (MQ=255)                                                                                                                 aaCCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGGCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCt                                  >  7:293483/1‑149 (MQ=255)                                                                                                                          gAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAgg                          >  3:65342/1‑148 (MQ=255)                                                                                                                          gAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAAGAACCTGCTGatt                          <  2:221266/148‑3 (MQ=255)                                                                                                                               tACGGAAATCGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATg                    <  8:111453/149‑1 (MQ=255)                                                                                                                                      aTAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGa             <  1:326044/149‑1 (MQ=255)                                                                                                                                        ccggCACCGACGCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGa             <  7:322402/145‑1 (MQ=255)                                                                                                                                        ccggCACCGACGCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGa             >  8:322402/3‑147 (MQ=255)                                                                                                                                        aGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGAcc           <  4:262048/149‑1 (MQ=255)                                                                                                                                           gCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATg        <  5:294136/149‑1 (MQ=255)                                                                                                                                            cACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATgc       >  5:18589/1‑149 (MQ=255)                                                                                                                                                gcgTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGCCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCa                                               >  1:29235/1‑105 (MQ=255)                                                                                                                                                gcgTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGCCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCa                                               <  2:29235/105‑1 (MQ=255)                                                                                                                                                 cgTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAg  >  6:50507/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      AAAGGTCCGTCGATTTAATGGGTTGATGGATGGCTTTTTCATCGAGCTCCTGTAGTTCCGGTACCCAGTGGCGAATGAACCGACCCTGTGGATCGAAGCGCTGCGACTGTGAAACCGGGTTGAAGATACGGAAATAGGGCACCGCGTCGGTGCCGGTGGAA‑GCACTCCACTGCCAACCCCCGTTGTTGGCGGCCAGGTCGCCGTCTATCAGGTGGCGCATGAAGTGCCGCTCGCCCTTGCGCCAGTCAATCAGCAGGTTCTTGCTGAGGAACATGGCCACGACCATGCGCAAG  >  NC_002947/857494‑857786

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