Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,524,463	Δ1 bp	6.9%	intergenic (+95/+18)	yncE → / ← ansP	PQQ‑like domain‑containing protein YncE/L‑asparagine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,524,463	0	A	.	6.9%	235.1 / 10.2	72	intergenic (+95/+18)	yncE/ansP	PQQ‑like domain‑containing protein YncE/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base A (38/29);  new base . (2/3);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGA  >  NC_000913/1524314‑1524611                                                                                                                                                       |                                                                                                                                                     cACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTa                                                                                                                                                      >  2:699416/1‑150 (MQ=255)  aCTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGAGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                                                                                     >  1:119875/1‑150 (MQ=255)     aaaCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                                                                                  >  1:1430501/1‑150 (MQ=255)         agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCa                                                                                                                                              <  1:445890/150‑1 (MQ=255)         agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCa                                                                                                                                              <  1:1151021/150‑1 (MQ=255)         agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCa                                                                                                                                              <  2:793637/150‑1 (MQ=255)          gcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCaa                                                                                                                                             >  1:163275/1‑150 (MQ=255)            agGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCaaat                                                                                                                                           <  1:924144/150‑3 (MQ=255)            agGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAAtt                                                                                                                                           <  2:572974/150‑1 (MQ=255)             ggAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAAttt                                                                                                                                          >  2:6504/1‑150 (MQ=255)                 gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCa                                                                                                                                      <  2:1003191/150‑1 (MQ=255)                 gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCa                                                                                                                                      <  2:1314545/150‑1 (MQ=255)                                   gtgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAg                                                                                                                    <  1:1328751/150‑1 (MQ=255)                                     gATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGc                                                                                                                  >  2:1138764/1‑150 (MQ=255)                                            aTTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                                                                           <  1:1097558/150‑1 (MQ=255)                                               gcgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                        >  2:639911/1‑150 (MQ=255)                                               gcgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                        >  2:711299/1‑150 (MQ=255)                                                   tGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTc                                                                                                    >  2:1309864/1‑150 (MQ=255)                                                    gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                                   <  1:893565/150‑1 (MQ=255)                                                      aaTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                                 >  2:272564/1‑150 (MQ=255)                                                      aaTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                                 >  1:657339/1‑150 (MQ=255)                                                      aaTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                                 >  1:1124588/1‑150 (MQ=255)                                                              tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                                                                         <  1:761438/150‑1 (MQ=255)                                                                 cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAgg                                                                                                                   >  2:592088/1‑121 (MQ=255)                                                                 cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAgg                                                                                                                   <  1:592088/121‑1 (MQ=255)                                                                 cGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                                                                                      >  2:1410064/1‑150 (MQ=255)                                                                   aGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCt                                                                                    >  1:90949/1‑150 (MQ=255)                                                                        gAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                                                                               >  2:1213998/1‑150 (MQ=255)                                                                         aaaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcg                                                                              <  1:185090/150‑1 (MQ=255)                                                                         aaaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcg                                                                              >  2:643594/1‑150 (MQ=255)                                                                          aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcga                                                                             >  1:403483/1‑150 (MQ=255)                                                                              ggTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCgacgac                                                                         >  1:1419585/1‑150 (MQ=255)                                                                               gTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACt                                                                        >  1:431361/1‑150 (MQ=255)                                                                                  agaCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgc                                                                     <  2:290679/150‑1 (MQ=255)                                                                                  agaCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgc                                                                     <  2:946006/150‑1 (MQ=255)                                                                                         cctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc                                                              >  2:1115695/1‑150 (MQ=255)                                                                                         cctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc                                                              <  1:1101755/150‑1 (MQ=255)                                                                                          ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt                                                             >  2:1329120/1‑150 (MQ=255)                                                                                            gcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgt                                                           >  1:285089/1‑150 (MQ=255)                                                                                             cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg                                                          <  2:1419585/150‑1 (MQ=255)                                                                                              ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                                                                               <  2:545656/128‑1 (MQ=255)                                                                                              ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                                                                               >  1:545656/1‑128 (MQ=255)                                                                                              ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa                                                         >  1:821495/1‑150 (MQ=255)                                                                                              ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa                                                         >  2:22360/1‑150 (MQ=255)                                                                                               ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGaa                                                        >  2:887665/1‑150 (MQ=255)                                                                                                   aTCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAttta                                                    >  1:508957/1‑149 (MQ=255)                                                                                                    tCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                                   >  1:396168/1‑150 (MQ=255)                                                                                                         gCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAAc                                              <  2:513435/150‑1 (MQ=255)                                                                                                           gTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAAcgc                                            >  1:473766/1‑150 (MQ=255)                                                                                                               ccGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGttt                                        >  2:633295/1‑150 (MQ=255)                                                                                                                cGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTg                                       <  2:505543/150‑1 (MQ=255)                                                                                                                 ggCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTgc                                      <  1:1410064/150‑1 (MQ=255)                                                                                                                   cACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTgcgc                                    <  2:1124588/150‑1 (MQ=255)                                                                                                                     cTTTCTTCAGTAATCTTAATTTGGACGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTACTGTTTTTCTTCATCTTCCTCTACGACTGGCGCAGTGCTGTTAATTTCAGCAACGCGTTTGCGcac                                  <  2:657339/150‑1 (MQ=255)                                                                                                                         cGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCa                              >  1:334151/1‑150 (MQ=255)                                                                                                                               tAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAg                        <  1:1080952/150‑1 (MQ=255)                                                                                                                                 aTCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc                      <  1:887665/150‑1 (MQ=255)                                                                                                                                 aTCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc                      <  1:272564/150‑1 (MQ=255)                                                                                                                                    tGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCACAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCtat                   >  2:1045127/1‑150 (MQ=255)                                                                                                                                       aTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAc                >  2:810989/1‑150 (MQ=255)                                                                                                                                       aTTTGCCCGGA‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAATTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAc                >  2:315483/1‑150 (MQ=255)                                                                                                                                        tttGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAcc               <  1:1213998/150‑1 (MQ=255)                                                                                                                                        tttGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAcc               >  2:46714/1‑150 (MQ=255)                                                                                                                                        attGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAcc               >  1:1161316/2‑150 (MQ=255)                                                                                                                                          tGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAg             >  1:1281937/1‑150 (MQ=255)                                                                                                                                           gCCCGTC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc            <  1:881068/150‑1 (MQ=255)                                                                                                                                           gCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc            >  1:1238349/1‑150 (MQ=255)                                                                                                                                            cccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa           <  2:163275/150‑1 (MQ=255)                                                                                                                                            cccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa           >  2:755611/1‑150 (MQ=255)                                                                                                                                            cccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa           >  2:908581/1‑150 (MQ=255)                                                                                                                                                    cGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATAcc    <  2:285089/150‑1 (MQ=255)                                                                                                                                                     gTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCg   <  2:821495/150‑1 (MQ=255)                                                                                                                                                      tAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGa  <  2:396168/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGA  >  NC_000913/1524314‑1524611

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