Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	353,100	G→A	G132D (GGC→GAC)	prpE →	propionyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	353,100	0	G	A	98.4%	218.8 / ‑3.9	65	G132D (GGC→GAC)	prpE	propionyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (34/29);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01

CGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGGCGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTC  >  NC_000913/352853‑353346                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       cGAACAAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGacgc                                                                                                                                                                                                                                                     >  2:26582/1‑251 (MQ=255)  gAACCAACTTGTGCCACAACGCCATCGACCGCTGGCGGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTGTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAAGGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCTTATTACTCTGCTGGCCTGCGCGCGCATTGACGCt                                                                                                                                                                                                                                                    <  1:275216/251‑1 (MQ=255)   aaCCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCt                                                                                                                                                                                                                                                    >  2:297569/1‑250 (MQ=255)     ccAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTAtt                                                                                                                                                                                                                                                 >  1:87736/1‑251 (MQ=255)         cTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACt                                                                                                                                                                                                                                             >  1:4225/1‑251 (MQ=255)              gCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCggtg                                                                                                                                                                                                                                        >  1:93767/1‑251 (MQ=255)                  caACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGgtgt                                                                                                                                                                                                                                    <  1:100730/251‑1 (MQ=255)                           cGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCGATTCACTCGGTGGTGTTTGGTGGAt                                                                                                                                                                                                                           <  2:413172/251‑1 (MQ=255)                                       ggAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGcaca                                                                                                                                                                                                               <  1:317456/251‑1 (MQ=255)                                            aaCAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTg                                                                                                                                                                                                          >  2:365881/1‑251 (MQ=255)                                              cAGCCAGAGGCGCTGGCGCCGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTggc                                                                                                                                                                                                        >  2:75162/1‑251 (MQ=255)                                                    gagGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGcgcg                                                                                                                                                                                                  >  2:358514/1‑251 (MQ=255)                                                        cgctggctctgATTGCCGTCTCTTCGGCAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCTCACAGCGTGGCGGCGCGAAtt                                                                                                                                                                                              <  2:63037/251‑1 (MQ=255)                                                          ctggcgctgATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATtg                                                                                                                                                                                             >  1:112106/1‑250 (MQ=255)                                                               gctgATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACg                                                                                                                                                                                       <  1:190262/251‑1 (MQ=255)                                                                     tGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAAc                                                                                                                                                                                 >  1:311306/1‑251 (MQ=255)                                                                      gCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAAc                                                                                                                                                                                 <  1:305072/250‑1 (MQ=255)                                                                        cGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCAATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATAGATGACGCTAAACCgg                                                                                                                                                                              >  2:120134/1‑251 (MQ=255)                                                                           ctctTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGc                                                                                                                                                                           >  2:7519/1‑251 (MQ=255)                                                                            tctTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCt                                                                                                                                                                          >  1:345103/1‑251 (MQ=255)                                                                             ctTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTg                                                                                                                                                                         >  1:151198/1‑251 (MQ=255)                                                                             ctTCGGAAACAGAACAAGCGCGCAACTTTACCTTTCGTCAGCTGCATGACGACGTGAACGCGGCGGCCTCCATGTTGCGTTCATTGGGTGTGCAGCGCGGCGCTCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGAGCATATTACTCTGCTGGACTGCGCGAGCATTGACGCTATTCACTCGGTGGTGATTGGTGGATTTGCCTCGCACAGCGTGTCGGCGCGAATTGATGCAGCTAAACCGGTGCTg                                                                                                                                                                         >  2:48381/1‑251 (MQ=255)                                                                                   aaaCAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGc                                                                                                                                                                           >  1:227442/1‑243 (MQ=255)                                                                                   aaaCAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGc                                                                                                                                                                           <  2:227442/243‑1 (MQ=255)                                                                                     aCAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGtctc                                                                                                                                                                 >  1:115212/1‑251 (MQ=255)                                                                                           gaagaGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTg                                                                                                                                                                         >  2:272762/1‑237 (MQ=255)                                                                                           gaagaGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTg                                                                                                                                                                         <  1:272762/237‑1 (MQ=255)                                                                                             agagCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATg                                                                                                                                                         <  2:98472/251‑1 (MQ=255)                                                                                             agagCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATg                                                                                                                                                         <  1:387030/251‑1 (MQ=255)                                                                                               agCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGcc                                                                                                                                                       >  1:379396/1‑251 (MQ=255)                                                                                               agCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGcc                                                                                                                                                       >  2:141054/1‑251 (MQ=255)                                                                                                gcgcACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCg                                                                                                                                                      <  2:379396/251‑1 (MQ=255)                                                                                                        ttacctttCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAgcgcgc                                                                                                                                              >  1:237396/1‑251 (MQ=255)                                                                                                        ttacctttCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAgcgcg                                                                                                                                               <  2:324588/250‑1 (MQ=255)                                                                                                           ccgttCGTCAGCCTCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTTGGTGTGCTGCGCGGCGATCGGGGGCTGGTGTATATGCCGATGAGTAGCGAAGCGCATATTACTCTGCTGGCCTGCCCTCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGCTTTGATTCCGCTAAACTGGGGCTGATTGTCTCGCCTGATGCCGGGTCGCGCGGt                                                                                                                                           <  2:17896/248‑1 (MQ=255)                                                                                                             tttCGTCAGCTGCATGACGAAGTGACCGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTgg                                                                                                                                         >  2:305129/1‑251 (MQ=255)                                                                                                                  tCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCaaa                                                                                                                                     >  1:290842/1‑250 (MQ=255)                                                                                                                   cAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTGTTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCaaaa                                                                                                                                    <  1:316379/250‑1 (MQ=255)                                                                                                                      cTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTg                                                                                                                                          <  1:78579/241‑1 (MQ=255)                                                                                                                      cTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTg                                                                                                                                          >  2:78579/1‑241 (MQ=255)                                                                                                                          aTGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTccc                                                                                                                            >  1:274608/1‑251 (MQ=255)                                                                                                                                             ggCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTcgacg                                                                                                         <  1:365881/251‑1 (MQ=255)                                                                                                                                                  cAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCg                                                                                                    <  2:93767/251‑1 (MQ=255)                                                                                                                                                             tCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAgg                                                                                          <  1:102930/250‑1 (MQ=255)                                                                                                                                                                ttGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGgcgcc                                                                                      >  2:70007/1‑250 (MQ=255)                                                                                                                                                                ttGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGgcgc                                                                                       >  1:28092/1‑250 (MQ=255)                                                                                                                                                                         cAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGcc                                                                             >  2:192343/1‑251 (MQ=255)                                                                                                                                                                                      gggTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTg                                                                <  2:189399/251‑1 (MQ=255)                                                                                                                                                                                      gggTGCTGGTGGATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTg                                                                <  2:81732/251‑1 (MQ=255)                                                                                                                                                                                       ggTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGc                                                               <  1:387935/251‑1 (MQ=255)                                                                                                                                                                                          gCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCtgg                                                            >  1:254463/1‑251 (MQ=255)                                                                                                                                                                                          gCTGGTCTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACCGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAACTCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCCCCAGCCACGCCATGTTTTGCtgg                                                            >  2:353937/1‑251 (MQ=255)                                                                                                                                                                                                gtATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATc                                                      <  1:69296/251‑1 (MQ=255)                                                                                                                                                                                                   tatGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCCTGTTTTGCTGGTGGATcgcg                                                   >  2:119659/1‑251 (MQ=255)                                                                                                                                                                                                          atgatTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCCTGTTTTGCTGGAGGATCGCGGGCTGGc                                            >  2:344484/1‑251 (MQ=255)                                                                                                                                                                                                                          atatTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTc                            <  1:358514/251‑1 (MQ=255)                                                                                                                                                                                                                          atatTACTCTGCTGGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTc                            >  1:147927/1‑251 (MQ=255)                                                                                                                                                                                                                            tttACTCTGCTGGCCTTCCCTCGCATGTTCTCTATTCACTCGGTGGTGTTTGGTGGATGTGCCTCGCCCATCGGGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGGTGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGTCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCa                           <  2:340195/249‑1 (MQ=255)                                                                                                                                                                                                                              tACTCTGCTGGCCGTCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTGGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGAGGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAgc                         <  2:4225/250‑1 (MQ=255)                                                                                                                                                                                                                                  cggcgggCCTGCGCGCGCATTGACGCTATTCCCTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCTGGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAgcggg                      <  2:313324/244‑1 (MQ=255)                                                                                                                                                                                                                                   tgctgGCCTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAgcgggcgg                   >  1:262422/1‑251 (MQ=255)                                                                                                                                                                                                                                         ccTGCGCGCGCATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTc             >  2:368780/1‑251 (MQ=255)                                                                                                                                                                                                                                          cggcgcgggcATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGGGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCg            <  2:345103/249‑1 (MQ=255)                                                                                                                                                                                                                                           tGCGCGCGCATTGGCGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCg            >  1:44674/1‑250 (MQ=255)                                                                                                                                                                                                                                             ggcgcgcATTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGAtt         <  1:349705/250‑1 (MQ=255)                                                                                                                                                                                                                                                    aTTGACGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGTCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGcgtc  <  1:211083/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       CGAACCAACTTGTGCCACAACGCCATCGACCGCTGGCTGGAGAAACAGCCAGAGGCGCTGGCGCTGATTGCCGTCTCTTCGGAAACAGAAGAAGAGCGCACCTTTACCTTTCGTCAGCTGCATGACGAAGTGAACGCGGTGGCCTCAATGTTGCGTTCATTGGGTGTGCAGCGCGGCGATCGGGTGCTGGTGTATATGCCGATGATTGCCGAAGCGCATATTACTCTGCTGGCCTGCGCGCGCATTGGCGCTATTCACTCGGTGGTGTTTGGTGGATTTGCCTCGCACAGCGTGGCGGCGCGAATTGATGACGCTAAACCGGTGCTGATTGTCTCGGCTGATGCCGGAGCGCGCGGTGGCAAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTC  >  NC_000913/352853‑353346

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