Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,612,649	T→A	5.2%	F149Y (TTT→TAT)	ung →	uracil‑DNA glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,612,649	0	T	A	5.2%	237.7 / 2.5	76	F149Y (TTT→TAT)	ung	uracil‑DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base T (34/38);  new base A (2/2);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.20e-01

AAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGG  >  NC_002947/1612501‑1612795                                                                                                                                                      |                                                                                                                                                    aaGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGtt                                                                                                                                                     >  8:285004/1‑149 (MQ=255) aaGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGtt                                                                                                                                                     <  3:200892/149‑1 (MQ=255)   ggAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTa                                                                                                                                                   >  2:192878/1‑149 (MQ=255)     aGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTAcc                                                                                                                                                 <  8:312970/149‑1 (MQ=255)              gcgACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑a                                                                                                                                         <  1:187496/147‑1 (MQ=255)              gcgACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑a                                                                                                                                         >  2:187496/1‑147 (MQ=255)                     gAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTATACCGACAGG‑ATTATCCAgg                                                                                                                                >  5:151179/1‑149 (MQ=255)                     gAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTATACCGACAGG‑ATTATCCAgg                                                                                                                                >  5:145494/1‑149 (MQ=255)                        cATGCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGAGGAACAAGTCGTTCTGTTGCTCAACAG‑GACAATGAGCGTCGAGCGCGCCAATGCCGCGTAGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAggtgg                                                                                                                             <  4:339328/149‑1 (MQ=255)                         aTCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAggtggt                                                                                                                            <  4:240236/149‑1 (MQ=255)                          tCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAggtggt                                                                                                                            >  2:336391/1‑148 (MQ=255)                          tCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAggtg                                                                                                                              <  6:62936/146‑1 (MQ=255)                          tCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAggtg                                                                                                                              >  5:62936/1‑146 (MQ=255)                           cccGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGGTGGGAATTGTTTACCGACAGG‑ATTATCCAggtggt                                                                                                                            >  4:205164/1‑147 (MQ=255)                           cccGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCa                                                                                                                          >  4:93328/1‑149 (MQ=255)                             cGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTAcc                                                                                                                                                 <  1:125652/125‑1 (MQ=255)                             cGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTAcc                                                                                                                                                 >  2:125652/1‑125 (MQ=255)                              gcaccccAGCCACGGCTACCTCCAGATCTCGGCGGAACAAGGCGTTCAGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCTATGCCGCGTCGCATGGCAAGAAAGCCTGGGAATTGTTTACCGACAGG‑ATTATCCACGTGGTCAGTg                                                                                                                       <  6:305632/146‑1 (MQ=255)                              gATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GCCCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAg                                                                                                                                 >  5:56287/1‑139 (MQ=255)                              gATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAg                                                                                                                                 <  6:56287/139‑1 (MQ=255)                              gATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTg                                                                                                                       <  1:336391/149‑1 (MQ=255)                               aTCCCCAGCCACGGCTACCTCCACAGCTGGGCGGAACAAGGCGTACTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATCCCAAGAAAGGCTCGGAATTGTTTACCGACAGC‑ATTATCCAGGAGGTCAGTg                                                                                                                       >  6:324735/1‑148 (MQ=255)                                   ccAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACACGGACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGc                                                                                                                    >  1:103944/1‑148 (MQ=255)                                    cAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAAt                                                                                                                 <  3:322114/149‑1 (MQ=255)                                     aGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGGATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATg                                                                                                                >  7:88016/1‑149 (MQ=255)                                     aGCCACGGCTAACTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATg                                                                                                                <  7:248539/149‑1 (MQ=255)                                             cTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACg                                                                                                        >  1:55698/1‑149 (MQ=255)                                                 cTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACgtggt                                                                                                    <  2:190566/149‑1 (MQ=255)                                                    cAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTa                                                                                                                                      <  2:258080/112‑1 (MQ=255)                                                    cAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTa                                                                                                                                      >  1:258080/1‑112 (MQ=255)                                                        gCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTg                                                                                             <  1:192878/149‑1 (MQ=255)                                                        gCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTATACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTg                                                                                             <  6:145494/149‑1 (MQ=255)                                                        gCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTATACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTg                                                                                             <  6:151179/149‑1 (MQ=255)                                                            ggCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTAt                                                                                         >  4:19495/1‑149 (MQ=255)                                                             gCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATg                                                                                        >  4:168959/1‑149 (MQ=255)                                                                gAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATgggg                                                                                     <  2:103944/149‑1 (MQ=255)                                                                 aaCAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATggcgg                                                                                    >  1:142694/1‑149 (MQ=255)                                                                   cAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGTGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTAt                                                                                         <  7:60318/142‑1 (MQ=255)                                                                   cAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCTACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGtt                                                                                                 <  2:293597/134‑1 (MQ=255)                                                                   cAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGgc                                                                                  <  7:285004/149‑1 (MQ=255)                                                                   cAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCCACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGtt                                                                                                 >  1:293597/1‑134 (MQ=255)                                                                   cAAGGCGATCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTAt                                                                                         >  8:60318/1‑142 (MQ=255)                                                                    aaGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATGCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGgcg                                                                                 >  8:326882/1‑149 (MQ=255)                                                                     aGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGgcgc                                                                                <  7:291356/149‑1 (MQ=255)                                                                      ggCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCa                                                                               <  3:19495/149‑1 (MQ=255)                                                                         gTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGc                                                                            <  3:93328/149‑1 (MQ=255)                                                                              gTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGcaca                                                                         >  8:39262/1‑147 (MQ=255)                                                                               ttGCTCAACAC‑GACCATGACCGTCGAGCGCGCCACTGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGcaca                                                                         <  8:313719/146‑1 (MQ=255)                                                                               ttGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAGCGTGGTGTTCCTGCTATGGGGGGCGCATGcaca                                                                         <  2:59602/146‑1 (MQ=255)                                                                               ttGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACagag                                                                      <  7:39262/149‑1 (MQ=255)                                                                               ttGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACaga                                                                       <  5:324735/148‑1 (MQ=255)                                                                               ttGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACag                                                                        <  2:55698/147‑1 (MQ=255)                                                                                 gCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGca                                                                           >  2:45918/1‑142 (MQ=255)                                                                                 gCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGca                                                                           <  1:45918/142‑1 (MQ=255)                                                                                   tCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGt                                                                                                                        >  2:209065/1‑95 (MQ=255)                                                                                   tCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGt                                                                                                                        <  1:209065/95‑1 (MQ=255)                                                                                      acac‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACgtg                                                                                                      <  7:245578/110‑1 (MQ=255)                                                                                      acac‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACgtg                                                                                                      >  8:245578/1‑110 (MQ=255)                                                                                              cATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGa                                                        >  1:204798/1‑149 (MQ=255)                                                                                                 gACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGGCAGTGAGCAATGCGCGAACGTGGTGTTCCTGCTGTGGGGGGCGCATGCACAGCGCAAGCAGAGGCTGATc                                                      >  8:342653/1‑148 (MQ=255)                                                                                                   ccGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTCGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGa                                                    <  4:1154/148‑1 (MQ=255)                                                                                                    cGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGaa                                                             <  6:4114/138‑1 (MQ=255)                                                                                                    cGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGaa                                                             >  5:4114/1‑138 (MQ=255)                                                                                                          gcgcgcCAATGCCGCGTCGCATGCCAAGAAAGGGTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCa                                            <  3:205164/149‑1 (MQ=255)                                                                                                            gcgcCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGagcaag                                                                  >  5:217779/1‑125 (MQ=255)                                                                                                            gcgcCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGagcaag                                                                  <  6:217779/125‑1 (MQ=255)                                                                                                              gcCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCa                                        >  3:218596/1‑149 (MQ=255)                                                                                                                      cgcgTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGttt                                <  2:366610/149‑1 (MQ=255)                                                                                                                      cgcgTCGCATGCCAAGAAAGGCTGGGAATTGTTAACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGttt                                >  6:100788/1‑149 (MQ=255)                                                                                                                          tCGCATGCCAAGAAAGGCTGGGGATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGaa                            <  8:88016/149‑1 (MQ=255)                                                                                                                              aTGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCg                        <  4:218596/149‑1 (MQ=255)                                                                                                                                gCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGt                      >  2:287570/1‑149 (MQ=255)                                                                                                                                   aaGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCg                        <  3:182745/144‑1 (MQ=255)                                                                                                                                   aaGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCg                        >  4:182745/1‑144 (MQ=255)                                                                                                                                     gAAAGGCTGGGAATTGTTTACCGACAGGAATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCAt                  >  1:220449/1‑149 (MQ=255)                                                                                                                                           cTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTcgc           <  2:142694/149‑1 (MQ=255)                                                                                                                                              ggAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTcgccgc        >  7:26839/1‑149 (MQ=255)                                                                                                                                                    gTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCgg  >  8:253007/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    AAGGAGCTGCAGCGCGACCTGAACATCCCGATCCCCAGCCACGGCTACCTCCAGAGCTGGGCGGAACAAGGCGTTCTGTTGCTCAACAC‑GACCATGACCGTCGAGCGCGCCAATGCCGCGTCGCATGCCAAGAAAGGCTGGGAATTGTTTACCGACAGG‑ATTATCCAGGTGGTCAGTGAGCAATGCCCGAACGTGGTGTTCCTGCTATGGGGGGCGCATGCACAGAGCAAGCAGAAGCTGATCGACGGTACCAAGCATCTGGTTTTGAAGTCGGTCCATCCTTCGCCGCTGTCGG  >  NC_002947/1612501‑1612795

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑