Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	167,295	T→C	5.2%	intergenic (+31/‑189)	mrcB → / → fhuA	fused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,295	0	T	C	5.2%	286.8 / 4.8	95	intergenic (+31/‑189)	mrcB/fhuA	fused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (50/40);  new base C (3/2);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

ATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATC  >  NC_000913/167155‑167439                                                                                                                                             |                                                                                                                                                 aTCAGTCTTCTCCGCCGCAGCAACAGCCGCAACAGCACCCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGATTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATcag                                                                                                                                          >  1:1182836/1‑147 (MQ=255) aTCAGGCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCAGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCg                                                                                                                                          <  2:519411/149‑1 (MQ=255)    tgTCTTCTCAGCCGCAGCAACAGCCTCAACAGCAACCTGCTCATCAAGAGCAGAAAGACAGAGACGTTGTAGCAGGTTTGATCAAGGATATGTTTGTTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                       <  1:575709/148‑1 (MQ=255)    aGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCGCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTa                                                                                                                                       <  2:37810/149‑1 (MQ=255)      tcttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaa                                                                                                                                     >  1:1100408/1‑149 (MQ=255)       cttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaaa                                                                                                                                    <  1:56146/149‑1 (MQ=255)        ttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAAt                                                                                                                                   <  2:1020210/149‑1 (MQ=255)          ctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCCGCCAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGAGATGTTGGGTAGGAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCTTCCGGTAATATaa                                                                                                                                 >  2:246125/1‑149 (MQ=255)          ctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATaa                                                                                                                                 >  1:932431/1‑149 (MQ=255)           tcAGCCGCAGCAACAGCCGCCACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGGGTAGCCGGTTGGATCATGGATATGTTTGGTAGTAATTAACATCGAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                >  1:664125/1‑149 (MQ=255)           tcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                >  2:1139176/1‑149 (MQ=255)             aGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACAt                                                                                                                              <  2:1139048/149‑1 (MQ=255)              cccGCAGCAACAGCCGCAACAGCAACCTGCTCATCAAGAGCAGAAAGACAGCGACGGTGGAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATc                                                                                                                             <  2:799755/148‑1 (MQ=255)               ccGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATcc                                                                                                                            >  1:303917/1‑149 (MQ=255)                cGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATccc                                                                                                                           >  2:427924/1‑149 (MQ=255)                 gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCa                                                                                                                          >  2:388254/1‑149 (MQ=255)                     caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                      <  2:289256/149‑1 (MQ=255)                     caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                      >  2:345953/1‑149 (MQ=255)                     caACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                      >  1:1122801/1‑149 (MQ=255)                       aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                    <  2:207314/149‑1 (MQ=255)                       aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGACTAGCCGGTAGGATCAAGGATATGTTTGGTAGAAATTAACATCTAAGCGTGAAATACCGGAAGGCGAGTTGCCAACCGGAAAAATAACATCCCAACTAAg                                                                                                                    >  1:209562/1‑149 (MQ=255)                           ccGCAACAGCAAGCTGCTCAGCAAGATCAGAAAGACTGCGACGGTGTAGCCGGTTCGATCAAGGATTTGTTTGGTAGTTATTATCATCTAAGCGTGAATTACGGTATGGCGAGTTGCCATCCGGTAAAATATCATCCCATCTAAGatat                                                                                                                <  2:702490/149‑1 (MQ=255)                           ccGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGatat                                                                                                                >  1:8247/1‑149 (MQ=255)                            cGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGTTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATAtt                                                                                                               >  1:972805/1‑149 (MQ=255)                             gCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTa                                                                                                              >  1:327116/1‑149 (MQ=255)                             gCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTa                                                                                                              >  1:1106539/1‑149 (MQ=255)                                  aGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGACGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACAct                                                                                                         >  1:1021659/1‑149 (MQ=255)                                       ccTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttctt                                                                                                    >  2:291176/1‑149 (MQ=255)                                       ccTGCGCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttctt                                                                                                    <  1:610379/149‑1 (MQ=255)                                         tGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCtttt                                                                                                  <  2:554550/149‑1 (MQ=255)                                          tcTCAGCAAGAGCAGAAAGACATCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGAAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTc                                                                                                 <  1:1036379/148‑1 (MQ=255)                                                cAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTgg                                                                                           >  1:552380/1‑149 (MQ=255)                                                  agagCAGAACGACAGCGACGGTGTCGCCGGGTGGATCAAGGTTATGTGTGGTTGTATTTAACATCTAAGCGTGAAATACCGGATGGCGAGGTGCCATCCGGTAAAATAACATCCCAGCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                         <  1:696255/149‑1 (MQ=255)                                                   gagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttg                                                                                        <  2:399064/149‑1 (MQ=255)                                                     gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCTGgttgtt                                                                                      >  2:646927/1‑149 (MQ=255)                                                     gCAGAAAGACAGAGACGGTTTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCTGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                      <  2:891729/149‑1 (MQ=255)                                                     ccAGAAAGACAGCGACGGGGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                      <  2:970657/148‑1 (MQ=255)                                                      cAGAAAGACAGCTACGGTGTAGCCGGTTGGATGAAGTATATTTATGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                     <  1:661571/149‑1 (MQ=255)                                                      cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                     >  1:1012914/1‑149 (MQ=255)                                                       aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGTTAGTAATTAACATCTAAGCGTTAAATACCTGATATCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACACCTTTCTTTTCATCTGGTTGTtta                                                                                    >  2:789274/1‑149 (MQ=255)                                                       aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                    >  1:464092/1‑149 (MQ=255)                                                        gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgt                                                                                          <  2:432990/142‑1 (MQ=255)                                                        gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgt                                                                                          >  1:432990/1‑142 (MQ=255)                                                        gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACACCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttat                                                                                   <  2:403662/149‑1 (MQ=255)                                                         aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                  >  2:593873/1‑149 (MQ=255)                                                           aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTACCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                                >  1:216562/1‑149 (MQ=255)                                                           aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAAAACCGGATGGCGCGTTGCCATCCCGTAAAATAACATCCCATCAAAGATATTAACCCTTTCTTTTCATCTGGTTCTtaattac                                                                                >  1:707572/1‑148 (MQ=255)                                                                 ctaCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTCAGCGTTAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                          <  2:540420/147‑1 (MQ=255)                                                                 cGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTATCTTTTCATCTGGTTGTTTCTTAACCCTTc                                                                          <  1:688291/149‑1 (MQ=255)                                                                  gACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                       <  1:560941/87‑1 (MQ=38)                                                                  gACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                       >  2:560941/1‑87 (MQ=38)                                                                   aCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAg                                                                        <  1:755086/149‑1 (MQ=255)                                                                      gtgtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCATTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGaa                                                                     >  2:1127734/1‑149 (MQ=255)                                                                        gtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATACCATCTAAGATATTAACCCTTTCTTTTAATCTGGTTGTTTAATAACCCTTCAGGAACg                                                                   >  1:1049610/1‑149 (MQ=255)                                                                           gCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTc                                                                >  2:421031/1‑149 (MQ=255)                                                                             cGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAg                                                              <  2:1106539/149‑1 (MQ=255)                                                                               gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGAt                                                            <  1:617216/149‑1 (MQ=255)                                                                               gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGAt                                                            >  1:16207/1‑149 (MQ=255)                                                                                ttGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGAtt                                                           <  2:303917/149‑1 (MQ=255)                                                                                  ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGc                                                         <  2:361085/149‑1 (MQ=255)                                                                                  ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGc                                                         <  1:297771/149‑1 (MQ=255)                                                                                  ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGc                                                         <  1:287228/149‑1 (MQ=255)                                                                                   gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCg                                                        >  2:618311/1‑149 (MQ=255)                                                                                    aTCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGt                                                       <  1:355708/149‑1 (MQ=255)                                                                                     tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTa                                                      >  1:564494/1‑149 (MQ=255)                                                                                      cAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTAc                                                     >  2:1090884/1‑149 (MQ=255)                                                                                       aaGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGATGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACTCTCAGATTGCGTAcc                                                    >  1:93632/1‑149 (MQ=255)                                                                                         ggATATGTTTTGTAGTAATTAACATATAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGc                                                  <  2:216562/149‑1 (MQ=255)                                                                                         ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGc                                                  >  2:401420/1‑149 (MQ=255)                                                                                            tatGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTg                                               <  2:1100408/149‑1 (MQ=255)                                                                                             atGATTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAAGCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGc                                              <  2:1182836/149‑1 (MQ=255)                                                                                                   ggTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAccc                                        <  1:354480/149‑1 (MQ=255)                                                                                                    gtagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCg                                       >  2:1125999/1‑149 (MQ=255)                                                                                                    gtagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCACATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCg                                       >  2:614499/1‑149 (MQ=255)                                                                                                       gtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACACTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGcaa                                    >  1:1158977/1‑147 (MQ=255)                                                                                                          aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCg                                 >  2:571302/1‑149 (MQ=255)                                                                                                             aaCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                              >  2:57136/1‑149 (MQ=255)                                                                                                              aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                             <  2:689722/149‑1 (MQ=255)                                                                                                              aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTACAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                             >  2:477421/1‑149 (MQ=255)                                                                                                              aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTGTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACACGCCAGCGTTTc                             >  2:863932/1‑149 (MQ=255)                                                                                                                aTCTAAGAGTGAAATACCGGAGGGCGAGTTGCCATCCGGTAAAATAACATCCCATGTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTAATAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGa                           <  2:217868/149‑1 (MQ=255)                                                                                                                 tCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGaa                          <  2:988302/149‑1 (MQ=255)                                                                                                                   tAAGCGTGAAATACCGGATGGCGAGTTGCCCTCCGGTAAAATAACATCCCAACTAAGATATTAACCCTTTCGTTTCATCTGGTTGTTTATTAACCCTTCCGGAACGCTCAGCTTGCGTACCGCTTGCGATCCCGCCAGCGTTTCGaaaa                        >  2:1074121/1‑147 (MQ=255)                                                                                                                   tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAata                        >  1:405820/1‑149 (MQ=255)                                                                                                                        gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                   >  2:79877/1‑149 (MQ=255)                                                                                                                        gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                   <  2:1072315/149‑1 (MQ=255)                                                                                                                         tGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatct                  >  1:935528/1‑149 (MQ=255)                                                                                                                             aTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTTACCCTTTCTTTTCATCTGGTTGTTTATTTACCCTTCAGGAACGCTCAGATTGAGTACCTCTTGCGAACCCGAAAGCGATTCGAATAttatcttata              >  2:431790/1‑148 (MQ=255)                                                                                                                             aTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatcttatc              >  2:942315/1‑149 (MQ=255)                                                                                                                                ccGGATGGGGTTTTGCCATCCGGTAAAATAACATCCCCTTTTAGTTATTAACCCTTTCGTTTGATCTGGTTGTTTATTAACCCTTCAGGAACGCGCAGTTTGCGTACCGCTTGCGTACCCGCCAGCGTTTGGAATATTATCTTATCttt           <  2:972805/149‑1 (MQ=255)                                                                                                                                ccGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTAAGATTGCGTACCGCTTGCGAACCAGCCAGCGTTTCGAATATTATCTTATCttt           >  1:603514/1‑149 (MQ=255)                                                                                                                                 cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGGTTGCGAAGCCGCCAGCGTTTCGAATTTTATCTTATCtttt          >  2:374253/1‑148 (MQ=255)                                                                                                                                 cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTGTTTTACCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATATTATCTTta          >  1:667583/1‑149 (MQ=255)                                                                                                                                         gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATc  <  1:620000/149‑1 (MQ=255)                                                                                                                                         gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATc  <  1:538936/149‑1 (MQ=255)                                                                                                                                            ttGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                             <  1:825909/119‑1 (MQ=255)                                                                                                                                            ttGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                             >  2:825909/1‑119 (MQ=255)                                                                                                                                             |                                                                                                                                                 ATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATC  >  NC_000913/167155‑167439

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