Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,803,206	Δ1 bp	9.1%	coding (108/972 nt)	PP_4231 ←	xanthine dehydrogenase accessory factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,803,206	0	A	.	9.1%	193.0 / 14.6	66	coding (108/972 nt)	PP_4231	xanthine dehydrogenase accessory factor
Reads supporting (aligned to +/- strand):  ref base A (27/32);  new base . (2/4);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

GCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCTCGGCGAGCTGCTGTACGTT  >  NC_002947/4803062‑4803352                                                                                                                                                   |                                                                                                                                                     gCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCAcgcg                                                                                                                                                  <  7:263207/149‑1 (MQ=255)       gCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCg                                                                                                                                           <  8:220123/149‑1 (MQ=255)         cATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCggg                                                                                                                                         <  2:80429/149‑3 (MQ=255)          aTAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCgggt                                                                                                                                        >  4:206188/1‑146 (MQ=255)            aaCGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCgggttc                                                                                                                                      <  3:206188/149‑6 (MQ=255)             aCGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCgggttcg                                                                                                                                     <  7:283905/149‑7 (MQ=255)               gCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAg                                                                                                                                   >  8:133442/1‑149 (MQ=255)               gCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAg                                                                                                                                   >  8:144815/1‑149 (MQ=255)               gCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAg                                                                                                                                   <  4:188506/149‑1 (MQ=255)                cACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAgg                                                                                                                                  <  2:252079/149‑1 (MQ=255)                 aCGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGt                                                                                                                                 >  3:36311/1‑149 (MQ=255)                    aTGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGc                                                                                                                              <  8:333858/149‑1 (MQ=255)                          aCCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCgccgcg                                                                                                                        <  4:36311/149‑6 (MQ=255)                             ggCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACtgt                                                                                                                     <  2:147738/149‑1 (MQ=255)                                  cgtgAATTCCCCGAGGGCCACGCGCTCGAGAAAGTCGTCTTCCCCAGACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACtgtg                                                                                                                    <  1:289728/144‑1 (MQ=255)                                  cGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGcac                                                                                                                 <  7:133442/148‑1 (MQ=255)                                    ggAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAgg                                                                                                                                  >  3:140364/1‑129 (MQ=255)                                    ggAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAgg                                                                                                                                  <  4:140364/129‑1 (MQ=255)                                        ttCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATa                                                                                                          >  5:128348/1‑149 (MQ=255)                                        ttCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑CCGCCGAGCCGTAGGTCGCGAGCACTGTGCACAACCATa                                                                                                          >  7:196839/1‑149 (MQ=255)                                          ccccGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGCCAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGCGCCGTAGGTGGCGAGCACTGTGCACAGCCaga                                                                                                          >  6:340158/1‑145 (MQ=255)                                          ccccGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGc                                                                                                                                        <  5:121564/117‑1 (MQ=255)                                          ccccGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATa                                                                                                          <  5:340158/147‑1 (MQ=255)                                           cccGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑ccc                                                                                                      >  7:232471/1‑149 (MQ=255)                                              gAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGa                                                                                                                            >  4:249094/1‑125 (MQ=255)                                                       gcgTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGtt                                                                                            >  1:65474/1‑147 (MQ=255)                                                       gcgTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGtt                                                                                            <  2:65474/147‑1 (MQ=255)                                                         gTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAg                                                                                        <  8:281680/149‑1 (MQ=255)                                                          ttCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCgggttc                                                                                                                                      <  4:284744/103‑6 (MQ=255)                                                          ttCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCgggttc                                                                                                                                      >  3:284744/1‑98 (MQ=255)                                                              agaAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCa                                                                                   >  1:195589/1‑149 (MQ=255)                                                              agaAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCa                                                                                   <  5:318247/149‑1 (MQ=255)                                                                  aGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTg                                                                                >  2:269785/1‑148 (MQ=255)                                                                   gtcgtcTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCa                                                                              <  2:195589/149‑1 (MQ=255)                                                                     cgtcTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTc                                                                                                  >  5:164144/1‑127 (MQ=255)                                                                     cgtcTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTc                                                                                                  <  6:164144/127‑1 (MQ=255)                                                                         tt‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAg                                                                                        <  6:90750/133‑1 (MQ=255)                                                                         tt‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAg                                                                                        >  5:90750/1‑133 (MQ=255)                                                                             cGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCt                                                                        <  8:232471/147‑1 (MQ=255)                                                                             cGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGa                                                                      >  6:111397/1‑149 (MQ=255)                                                                              gACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGa                                                                      <  1:52460/149‑1 (MQ=255)                                                                              gACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGa                                                                      >  2:52460/1‑148 (MQ=255)                                                                               acacAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGc                                                                               >  5:27565/1‑138 (MQ=255)                                                                               acacAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGc                                                                               <  6:27565/138‑1 (MQ=255)                                                                                   aGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGAt                                                                 >  3:266967/1‑148 (MQ=255)                                                                                     ccACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCAc                                                              <  1:269785/149‑1 (MQ=255)                                                                                        ccGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTg                                                           <  4:266967/149‑1 (MQ=255)                                                                                                cGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGc                                                                                                                   >  7:22597/1‑86 (MQ=255)                                                                                                cGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGc                                                                                                                   <  8:22597/86‑1 (MQ=255)                                                                                                          cACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCAc                                         >  1:306046/1‑149 (MQ=255)                                                                                                           aCTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCgag                                                  <  6:240615/139‑1 (MQ=255)                                                                                                           aCTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCgag                                                  >  5:240615/1‑139 (MQ=255)                                                                                                           aCTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCa                                          <  1:276739/147‑1 (MQ=255)                                                                                                           aCTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCa                                          >  2:276739/1‑147 (MQ=255)                                                                                                                ctgCTGGCGGTGAACGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCtt                                     <  6:22318/145‑1 (MQ=255)                                                                                                                ctgCTGGCGGTGAACGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCtt                                     >  5:22318/3‑147 (MQ=255)                                                                                                                ccGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTa                                   <  2:306046/149‑1 (MQ=255)                                                                                                                 tgCTGGCGGTGAACGCCAGCGGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTg                                                           <  8:86835/123‑1 (MQ=255)                                                                                                                 tgCTGGCGGTGAACGCCAGCGGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTg                                                           >  7:86835/2‑124 (MQ=255)                                                                                                                            aCCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCa                        <  6:128348/148‑1 (MQ=255)                                                                                                                                ccAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGc                                                                          <  8:177850/94‑1 (MQ=255)                                                                                                                                ccAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGc                                                                          >  7:177850/1‑94 (MQ=255)                                                                                                                                 cagcagCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTATTCCGCCTTCAGCTCgg                  >  7:25617/1‑149 (MQ=255)                                                                                                                                           cccGGCGCTCGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCTCGCTGAgctgaag        <  8:196839/149‑4 (MQ=255)                                                                                                                                           cccGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGc                                                                                          >  8:159370/1‑67 (MQ=255)                                                                                                                                           cccGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGc                                                                                          <  7:159370/67‑1 (MQ=255)                                                                                                                                                   aCGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCTCGGCGAGCTGCTGTACGtt  >  8:125544/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTT‑‑CGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATA‑CCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCTCGGCGAGCTGCTGTACGTT  >  NC_002947/4803062‑4803352

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