Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,803,225	T→G	6.5%	N45T (AAC→ACC)	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,225	0	T	G	6.5%	205.1 / 4.6	71	N45T (AAC→ACC)	PP_4231	xanthine dehydrogenase accessory factor
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base G (3/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCC  >  NC_002947/4803077‑4803371                                                                                                                                                     |                                                                                                                                                    cGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGtt                                                                                                                                                    <  1:420434/150‑1 (MQ=255)  ggCAGGGGAGCGCGCGGGCGGGGGCAGCTTCCGTGGCCGTAATGCACGAGAGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTc                                                                                                                                                   <  2:373908/150‑1 (MQ=255)      ggCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCg                                                                                                                                               >  1:982322/1‑150 (MQ=255)         gaacgcgcgaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCca                                                                                                                                            <  1:807612/150‑1 (MQ=255)          aacgcgcgaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCcag                                                                                                                                           >  1:471055/1‑150 (MQ=255)             gcgcgaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCcagcag                                                                                                                                        >  2:1291476/1‑150 (MQ=255)                cgaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCTGCTGGCGGTGAACGCcagcggcca                                                                                                                                     <  2:562889/150‑6 (MQ=255)                 gaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCGACGCGTTCGCGAAAGTCGTCTTCGACACAGCGACCGGACAGCGCGCGCAGCCACTGGCCGCTGGCGTTCACAGCCAGCCGCGaa                                                                                                                                    >  2:424939/1‑150 (MQ=255)                 gaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGaa                                                                                                                                    >  2:1067984/1‑150 (MQ=255)                     gcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCg                                                                                                                                <  2:471055/150‑1 (MQ=255)                                ccGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGACACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCTACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACgcg‑gc                                                                                                                    <  1:1105027/150‑1 (MQ=255)                                   tCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGcc                                                                                                                 >  1:569161/1‑150 (MQ=255)                                      ccATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAg                                                                                                              <  2:956484/150‑1 (MQ=255)                                      ccATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCggg                                                                                                              <  2:1277258/150‑3 (MQ=255)                                       cATAACGCACGATGGCAACCGGCCCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGC‑CGCGTGCGCCgggt                                                                                                             <  1:348347/150‑4 (MQ=255)                                             gCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACGGCcagcag                                                                                                                                        <  1:1352405/118‑1 (MQ=255)                                             gCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACGGCcagcag                                                                                                                                        >  2:1352405/1‑118 (MQ=255)                                              cACGATGGCAACCGGCTGCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCCTCTTCGACACAGCCACCCGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGt                                                                                                      <  2:1293936/150‑1 (MQ=255)                                                       aaCCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCTGCTGGCGGTGAACGCcagcggcca                                                                                                                                     >  2:1064368/1‑106 (MQ=255)                                                       aaCCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCTGCTGGCGGTGAACGCcagcggcca                                                                                                                                     <  1:1064368/111‑6 (MQ=255)                                                       aaCCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCg                                                                                                                                >  2:430247/1‑116 (MQ=255)                                                       aaCCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCg                                                                                                                                <  1:430247/116‑1 (MQ=255)                                                          cGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACtgt                                                                                          <  2:982322/150‑1 (MQ=255)                                                                cgcgAATTCCCCGAGGGCCACGCGTTCGAGAAAGCCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACaa                                                                                    <  1:424939/147‑1 (MQ=255)                                                                cGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAg                                                                                                >  1:7768/1‑138 (MQ=255)                                                                cGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAg                                                                                                <  2:7768/138‑1 (MQ=255)                                                                  ggAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAAcc                                                                                  <  1:1133236/150‑1 (MQ=255)                                                                      ttCCCGGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATAc                                                                              >  2:821093/1‑150 (MQ=255)                                                                         cccGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCg                                                                           >  2:494331/1‑150 (MQ=255)                                                                         cccGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCg                                                                           >  1:170997/1‑150 (MQ=255)                                                                          ccGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGc                                                                          >  1:802762/1‑150 (MQ=255)                                                                          ccGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGc                                                                          >  1:932831/1‑150 (MQ=255)                                                                                  cACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACCAACACTCCAGTTACTGGCCGCTTGCGGTTACTGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGtt                                                                  >  1:739761/1‑150 (MQ=255)                                                                                   aCGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTg                                                                 <  2:285408/150‑1 (MQ=255)                                                                                    cgcgTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGc                                                                >  2:893996/1‑150 (MQ=255)                                                                                      cgTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAg                                                              <  2:932175/150‑1 (MQ=255)                                                                                        ttCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAg                                                                                                        >  1:1011263/1‑106 (MQ=255)                                                                                        ttCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAg                                                                                                        <  2:1011263/106‑1 (MQ=255)                                                                                        ttCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGa                                                            <  2:1097528/150‑1 (MQ=255)                                                                                           gagaAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCa                                                         >  1:142716/1‑150 (MQ=255)                                                                                           gagaAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCa                                                         <  2:86295/150‑1 (MQ=255)                                                                                                   cgtcTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAg                                                 <  2:723621/150‑1 (MQ=255)                                                                                                    gtcTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGc                                                <  1:821093/150‑1 (MQ=255)                                                                                                       ttCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTg                                             >  1:325681/1‑150 (MQ=255)                                                                                                            cacaGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGtt                                                                                                                                                    >  2:104331/1‑43 (MQ=255)                                                                                                            cacaGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGtt                                                                                                                                                    <  1:104331/43‑1 (MQ=255)                                                                                                            cacaGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGCGCCTTAGTTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGa                                        >  1:88980/1‑150 (MQ=255)                                                                                                             ccaGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGAt                                       <  1:893254/149‑1 (MQ=255)                                                                                                              caGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCTGTAGGTGGAGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATc                                      <  2:1326959/150‑1 (MQ=255)                                                                                                               aGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCa                                     >  2:1290240/1‑150 (MQ=255)                                                                                                                  cACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCt                                  <  2:569161/150‑1 (MQ=255)                                                                                                                          aGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTg                                 <  1:192720/143‑1 (MQ=255)                                                                                                                          aGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTg                                 >  2:192720/1‑143 (MQ=255)                                                                                                                           gCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCga                         <  1:1290240/150‑1 (MQ=255)                                                                                                                            cGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCgag                        <  1:405540/150‑1 (MQ=255)                                                                                                                            cGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACACGGCGAACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCgag                        >  1:1106084/1‑150 (MQ=255)                                                                                                                                   aGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGAc                                                           <  1:95953/108‑1 (MQ=255)                                                                                                                                   aGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGAc                                                           >  2:95953/1‑108 (MQ=255)                                                                                                                                   aGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGa                                            <  1:120456/123‑1 (MQ=255)                                                                                                                                   aGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGa                                            >  2:120456/1‑123 (MQ=255)                                                                                                                                      cACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGAc                                           >  2:479425/1‑121 (MQ=255)                                                                                                                                      cACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGAc                                           <  1:479425/121‑1 (MQ=255)                                                                                                                                       aCTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGc             >  1:231904/1‑150 (MQ=255)                                                                                                                                        cTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCt            <  1:1067984/150‑1 (MQ=255)                                                                                                                                        cTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCt            <  1:1291476/150‑1 (MQ=255)                                                                                                                                        cTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCt            >  1:158276/1‑150 (MQ=255)                                                                                                                                             tgCTGGCGGTGAACGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACt       >  1:608322/2‑150 (MQ=255)                                                                                                                                              gCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTc      >  2:352560/1‑150 (MQ=255)                                                                                                                                                  gCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCgag                        >  2:78580/1‑128 (MQ=255)                                                                                                                                                  gCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCgag                        <  1:78580/128‑1 (MQ=255)                                                                                                                                                  gCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGAtgct                   <  1:893037/133‑1 (MQ=255)                                                                                                                                                  gCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGAtgct                   >  2:893037/1‑133 (MQ=255)                                                                                                                                                  gCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGcc  <  2:1368653/150‑1 (MQ=255)                                                                                                                                                  gCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGcc  <  2:1201412/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCC  >  NC_002947/4803077‑4803371

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