Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,803,229	C→A	5.6%	V44L (GTG→TTG) ‡	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,229	0	C	A	5.6%	233.7 / 2.1	71	V44L (GTG→TTG) ‡	PP_4231	xanthine dehydrogenase accessory factor
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base A (2/2);  total (33/38)
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.73e-01

GGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCT  >  NC_002947/4803082‑4803378                                                                                                                                                    |                                                                                                                                                       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                                                                                                                                               >  2:1352405/1‑118 (MQ=255)                                        gCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACGGCcagcag                                                                                                                                               <  1:1352405/118‑1 (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)                                                                                                       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                                                   >  2:120456/1‑123 (MQ=255)                                                                                                                              aGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGa                                                   <  1:120456/123‑1 (MQ=255)                                                                                                                                 cACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGAc                                                  <  1:479425/121‑1 (MQ=255)                                                                                                                                 cACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGAc                                                  >  2:479425/1‑121 (MQ=255)                                                                                                                                  aCTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGc                    >  1:231904/1‑150 (MQ=255)                                                                                                                                   cTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCt                   <  1:1291476/150‑1 (MQ=255)                                                                                                                                   cTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCt                   >  1:158276/1‑150 (MQ=255)                                                                                                                                   cTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCt                   <  1:1067984/150‑1 (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)                                                                                                                                                 tCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCa     <  2:170997/150‑1 (MQ=255)                                                                                                                                                    ccGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACGGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCt  >  2:558429/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAACCCGGCGCACGCG‑GCGCCGAGCCGTAGGTGGCGAGCACTGTGCACAACCATACCCGCTCGCCGTTGCAGGACCACTGCAAAGCCTGACGGATCACCTGTAGATCGAGATGCTGCACGCTTTACTCCGCCTTCAGCT  >  NC_002947/4803082‑4803378

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