Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,809,454	T→G	V305V (GTT→GTG)	waaA →	KDO transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,809,454	0	T	G	98.4%	214.8 / NA	63	V305V (GTT→GTG)	waaA	KDO transferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (35/27);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAC  >  NC_000913/3809213‑3809696                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   gccgccACTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGAGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGAGGGCGATACg                                                                                                                                                                                                                                           >  2:172815/1‑251 (MQ=255)    gccAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGGGATACGATg                                                                                                                                                                                                                                        >  2:181337/1‑251 (MQ=255)            gCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGtt                                                                                                                                                                                                                                <  2:268791/251‑1 (MQ=255)                          ggTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTacg                                                                                                                                                                                                                  >  1:179341/1‑249 (MQ=255)                           gTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTacgg                                                                                                                                                                                                                 <  1:305821/251‑4 (MQ=255)                            tATGGATTGCCACCAGCACTCACGAAGGCGCAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCGCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGCTACGATGGGCGGGTTGATGTTACTGTACGGc                                                                                                                                                                                                                >  2:80493/1‑251 (MQ=255)                             aTGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCg                                                                                                                                                                                                                                    <  1:323928/230‑1 (MQ=255)                             aTGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCg                                                                                                                                                                                                                                    >  2:323928/1‑230 (MQ=255)                               ggATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCAtt                                                                                                                                                                                                             <  1:424813/251‑1 (MQ=255)                                   tGCCACCAGCAATCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGGGGGCGATACGATGGGCGTGTTGATGTTACTGTACGGCCTTGCCg                                                                                                                                                                                                         >  2:393255/1‑251 (MQ=255)                                     ccaccaGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGAt                                                                                                                                                                                                       >  1:159359/1‑251 (MQ=255)                                           gCACTCACGAAGGCGAAGATAGTGTGGTGATCGCCGCACATCAGGCATTGCTACAGCAACTCCCGAATTTATTGCTCATCCTGGTACCCCGGCAACCAGAACGCTGCCTGGATGCGATTAAACGTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCTCCAGCACGCAGGTGGTGGTGGGCGATACGATGGGCGAGTTGATATTACTGTACGTCATTGCCGATCTCGcc                                                                                                                                                                                                 <  2:152495/251‑1 (MQ=255)                                            cACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCt                                                                                                                                                                                                >  1:314245/1‑251 (MQ=255)                                                 aCGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTTGTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGCACGGCATTGCCGATCTCGCCTttgtt                                                                                                                                                                                           <  1:242195/251‑1 (MQ=255)                                                             agagTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTAGTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTGGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCGTTGTTGGCGGTTCACTg                                                                                                                                                                               <  2:179341/251‑1 (MQ=255)                                                                 tgtgGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTg                                                                                                                                                                           >  2:125539/1‑251 (MQ=255)                                                                     gtgATCGCCGCACATCAGGCATTGTGACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCTCCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACg                                                                                                                                                                       <  2:256107/251‑1 (MQ=255)                                                                       gATCGCCGCACATCAGGCATTGTTCCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACgtg                                                                                                                                                                     >  2:296621/1‑251 (MQ=255)                                                                            ccgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTggg                                                                                                                                                                >  1:258740/1‑251 (MQ=255)                                                                             cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                                                                                               >  2:221749/1‑251 (MQ=255)                                                                             cgcACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTCGTAACCCGTCATCCAGAACGCTTCCCGGCTGCGATTCATCTTGTCAGCCAGGCTGGACTAAGCTATATCACCAGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGGTGAACGTGGTGGGc                                                                                                                                                               >  1:415969/1‑251 (MQ=255)                                                                                        cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                                                                                               >  2:131547/1‑240 (MQ=255)                                                                                        cATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGc                                                                                                                                                               <  1:131547/240‑1 (MQ=255)                                                                                             ttACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGc                                                                                                                                               >  1:267692/1‑251 (MQ=255)                                                                                               aCAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTg                                                                                                                                             >  1:259727/1‑251 (MQ=255)                                                                                                   cAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGc                                                                                                                                         >  2:252935/1‑251 (MQ=255)                                                                                                       tCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGcaca                                                                                                                                      <  2:159359/250‑1 (MQ=255)                                                                                                        cccGAATTTATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCATCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCTCTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCacac                                                                                                                                     <  2:144006/250‑1 (MQ=255)                                                                                                               ttATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTcc                                                                                                                             >  1:68380/1‑251 (MQ=255)                                                                                                                tATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCg                                                                                                                            >  1:219828/1‑251 (MQ=255)                                                                                                                tATTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAAGCTTGTCCGCCAGGCTGGACTAAGCTATATCTCACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCTCACGCTATTcc                                                                                                                             <  2:190323/250‑1 (MQ=255)                                                                                                                 aTTGCTCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCgg                                                                                                                           >  1:375838/1‑251 (MQ=255)                                                                                                                      tCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCt                                                                                                                                              <  1:373584/227‑1 (MQ=255)                                                                                                                      tCATCCTGGTACCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCt                                                                                                                                              >  2:373584/1‑227 (MQ=255)                                                                                                                            gggTACCCCGTCATCCAGAAAGCTTCCCGGATTCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCCCTCCAGGGGAAGGCCCCTCGACCAGCACGCAGGTTGTGGTGGGCGATTCGATGGGCGAGTTGATGTTACTGGTCGGCATGTCCGATGCCGTCTTTGTTGGCGATTCTCTGGTTGAGCGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACATGCTATTCCGGGATTGATgggg                                                                                                                <  2:415969/250‑1 (MQ=255)                                                                                                                                aCCCCGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGc                                                                                                            >  1:277970/1‑251 (MQ=255)                                                                                                                                  cccGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCAt                                                                                                          <  2:348174/251‑1 (MQ=255)                                                                                                                                   ccGTCATCCAGAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATa                                                                                                         <  2:430504/251‑1 (MQ=255)                                                                                                                                           cagAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAAc                                                                                                 <  2:258740/249‑1 (MQ=255)                                                                                                                                             gAACGCTTCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACtt                                                                                               >  1:131849/1‑251 (MQ=255)                                                                                                                                                    tCCCGGATGCGATTAATCTTGTCCGCCAGGCTGGCCTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGGGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGATGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATCCTTTTAACTTTAAAGa                                                                                         <  2:40555/250‑1 (MQ=255)                                                                                                                                                               aTTAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGcgcgcg                                                                             <  2:122361/251‑1 (MQ=255)                                                                                                                                                                ttAATCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg                                                                            >  2:8376/1‑251 (MQ=255)                                                                                                                                                                     cTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCt                                                                          <  1:241334/248‑1 (MQ=255)                                                                                                                                                                     cTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCt                                                                          >  2:241334/1‑248 (MQ=255)                                                                                                                                                                               aGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGCAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCAAACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGc                                                             >  1:55225/1‑251 (MQ=255)                                                                                                                                                                                     gACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTg                                                       >  1:437273/1‑251 (MQ=255)                                                                                                                                                                                          aGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg                                                                            >  2:426253/1‑225 (MQ=255)                                                                                                                                                                                          aGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCgg                                                                            <  1:426253/225‑1 (MQ=255)                                                                                                                                                                                            cTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGAttaccg                                                >  1:67943/1‑251 (MQ=255)                                                                                                                                                                                                atCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGAttaccgttac                                            >  1:235467/1‑251 (MQ=255)                                                                                                                                                                                                  cacaCGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGGCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGAttaccgttacc                                           <  1:169485/250‑1 (MQ=255)                                                                                                                                                                                                     acGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATg                                       <  1:167191/251‑1 (MQ=255)                                                                                                                                                                                                      cgcgctTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGCTGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGACCGTGGTGGGCATAATCCGCTGGAAGCTGCCGAACACGCTATTCCGGTATTGGTGGGGCCGCATACTTTTACCTTTAAAGCAATTTGCGCGCGGCTGGCGCAGGAAAGCGGGCTGGTTACCGTTACCGaggc                                      >  2:52437/5‑248 (MQ=255)                                                                                                                                                                                                           ttCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcgt                                 >  2:414311/1‑249 (MQ=255)                                                                                                                                                                                                            tCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcg                                  <  1:177895/249‑2 (MQ=255)                                                                                                                                                                                                            tCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGACGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAc                                   <  1:331572/248‑1 (MQ=255)                                                                                                                                                                                                               ggggAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcgtcgct                             <  1:289548/251‑7 (MQ=255)                                                                                                                                                                                                               ggggAAGTCCCCTCCACCAGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCAcgtcgct                             >  1:113130/1‑245 (MQ=255)                                                                                                                                                                                                                              cccaGCACGCAGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTc              <  1:181337/250‑1 (MQ=255)                                                                                                                                                                                                                                      gCAGGTTGTGGTGGGAGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGGACGTGGTGGGGATAATCCGCTGGACGCTGCCGCACACGCTATTCCGGCATTGAGGGGGCCGCAGACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGGGCAGGCAAGCGGGCTGATTAACGTTACCGATGCAACGTCGCTTGACAAAGAGGTTTCCTCTTTa       >  2:358669/1‑250 (MQ=255)                                                                                                                                                                                                                                        aGGTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTc    >  1:17515/1‑251 (MQ=255)                                                                                                                                                                                                                                          gTTGTGGTGGGCGATACGATGGGCGAGTTGATGTTACTGTACGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACGTCGCTTGCAAAAGAGGTTTCCTCTTTACTCAc  <  2:131849/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   GCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTACAGCAATTCCCGAATTTATTGCTCATCCTGGTACCCCGTCATCCGGAACGCTTCCCGGATGCGATTAACCTTGTCCGCCAGGCTGGACTAAGCTATATCACACGCTCTTCAGGGGAAGTCCCCTCCACCAGCACGCAGGTTGTGGTTGGCGATACGATGGGCGAGTTGATGTTACTGTATGGCATTGCCGATCTCGCCTTTGTTGGCGGTTCACTGGTTGAACGTGGTGGGCATAATCCGCTGGAAGCTGCCGCACACGCTATTCCGGTATTGATGGGGCCGCATACTTTTAACTTTAAAGACATTTGCGCGCGGCTGGAGCAGGCAAGCGGGCTGATTACCGTTACCGATGCCACTACGCTTGCAAAAGAGGTTTCCTCTTTACTCAC  >  NC_000913/3809213‑3809696

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑