Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,602,509	C→T	W117* (TGG→TAG)	yjjP ←	putative succinate exporter YjjP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,602,509	0	C	T	97.9%	168.5 / ‑3.3	48	W117* (TGG→TAG)	yjjP	putative succinate exporter YjjP
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (20/27);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

GGCCATCGCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAG  >  NC_000913/4602265‑4602750                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  ggcCATCGCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGTCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTg                                                                                                                                                                                                                                             <  2:208545/251‑1 (MQ=255)                     tATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                                                                                                                                                                         >  2:70743/1‑251 (MQ=255)                        tGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTgg                                                                                                                                                                                                                      >  2:179948/1‑251 (MQ=255)                           tGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTCAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGAAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCt                                                                                                                                                                                                                    >  2:361211/1‑250 (MQ=255)                           ggAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGTTGGCGGCGAAAGCGGTAAGGCAAAAGTGGATCTGGGGATGAAGATGACGTGGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCAGTAAGGCTACCAGCTATCTTGGGTAACTTAATGGCTGAATTTGGCTa                                                                                                                                                                                                                   <  2:354361/250‑1 (MQ=255)                             aaaGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTaa                                                                                                                                                                                                                  >  2:122407/1‑250 (MQ=255)                                   gggAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGttt                                                                                                                                                                                                           >  2:196299/1‑251 (MQ=255)                                               aGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCttt                                                                                                                                                                                               <  1:53110/251‑1 (MQ=255)                                                 cAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGt                                                                                                                                                                                             <  2:17681/251‑1 (MQ=255)                                                  aaTCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAGCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGATGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTa                                                                                                                                                                                            <  2:339366/251‑1 (MQ=255)                                                  aaTCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGTTATCTTGGGTAACGTAATGGCTg                                                                                                                                                                                                                             <  2:132608/218‑1 (MQ=255)                                                  aaTCCGGAAATGGTCGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGCTATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTg                                                                                                                                                                                                                             >  1:132608/1‑218 (MQ=255)                                                     ccGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAt                                                                                                                                                                                          <  2:330633/250‑1 (MQ=255)                                                     ccGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATc                                                                                                                                                                                         <  1:190489/251‑1 (MQ=255)                                                            tggtggtggCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTataccagc                                                                                                                                                                                     >  2:131633/1‑241 (MQ=255)                                                                  gggcggcGAAAGCGGTAAGGCAAAAGTTGATCGGTGGATGAAGATGACGTTGGGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTGTCTCTACGCCTTTGTAATCCAGCAGATGATGc                                                                                                                                                                            <  2:321590/250‑1 (MQ=255)                                                                    gcggcgAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTc                                                                                                                                                                          >  1:253257/1‑251 (MQ=255)                                                                     cggcgAAAGCGGTAAGGCAAAAGTTGATCGGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTGGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGAGGCTcc                                                                                                                                                                         <  2:289145/251‑1 (MQ=255)                                                                               ggTAAGGCACCAGTTGATCTGAGGCTGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGAAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATc                                                                                                                                                                >  2:347401/1‑250 (MQ=255)                                                                                  aaGGCGAAAGTTGAGCTGTGGATGAAGATGACGTTGTGCCAGCAGCTAGCGGATATACATAGCGGGCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGCAAGGCCAACCATTAAGGCGACCAGCTATCTTGGATAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATTCTCCGCAAGAATcaca                                                                                                                                                             >  2:14776/1‑250 (MQ=255)                                                                                   aGGCAAAAGCTGATCTGTGGATGAATATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTTCAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTGGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGTCTTTGTAATCCAGCAGATGATGCACCGCTAGAATCACAAt                                                                                                                                                           <  2:25795/251‑1 (MQ=255)                                                                                    ggCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg                                                                                                                                                          >  1:53539/1‑251 (MQ=255)                                                                                     gCAAAAGTGGATCTGGGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtgt                                                                                                                                                         <  2:137981/251‑1 (MQ=255)                                                                                          aGTTGATCTGTGGATGAAGATGACGTGGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCGTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCGTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTgg                                                                                                                                                    <  2:100846/251‑1 (MQ=255)                                                                                                                            cagcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGAt                                                                                                                  <  1:299759/251‑1 (MQ=255)                                                                                                                              gcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCg                                                                                                                <  1:122407/251‑1 (MQ=255)                                                                                                                                        ataCATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTa                                                                                                       <  1:129001/250‑1 (MQ=255)                                                                                                                                                          ggCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCAtt                                                                                    <  2:88292/251‑1 (MQ=255)                                                                                                                                                           gCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTg                                                                                   >  1:92837/1‑251 (MQ=255)                                                                                                                                                            cAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGc                                                                                  >  1:3721/1‑251 (MQ=255)                                                                                                                                                                  agGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCt                                                                            <  1:153537/251‑1 (MQ=255)                                                                                                                                                                      gatgaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTa                                                                         >  2:193396/1‑250 (MQ=255)                                                                                                                                                                      gatgaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTa                                                                         >  1:61075/1‑250 (MQ=255)                                                                                                                                                                         gaCGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAg                                                                     <  1:185568/251‑1 (MQ=255)                                                                                                                                                                            ggCACCATCCCCGACACCGTTCTTGAGTTTACAGAAACCGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtc                                                                                                                                                                                                          >  2:233247/1‑115 (MQ=255)                                                                                                                                                                            ggCACCATCCCAGCCACCGTTATTGAGGTTACAGAAACAGGCGCAAGAAAGGCCAACCAGTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtc                                                                                                                                                                                                          <  1:233247/115‑1 (MQ=255)                                                                                                                                                                                       agaaaacGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTTGGTACCGTAATGGCTGAATCTGGCGAAATCGTCGCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAGGTGCTGGACCTCAGGCACCACATGCATATTAATGCCGTGATCGTGATTTTTACGTGTCGGTGTCAGGCATCGCCCATCTGTATTAGTTGTTAGCACTATg                                                       <  2:159041/245‑1 (MQ=255)                                                                                                                                                                                                             gAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCTTTTCTCTACTCCTTTGTAATCCAGCAGATGATTCTCCGCAAGAATCGCAATGTGCTGGGCTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTAGGGGGCGATGTCGGGCATTGCCCATCTGTAATAGTGGGCAGCACTAGGGCGTTCGTAGTGATTGAACttt                                 <  2:331307/251‑1 (MQ=255)                                                                                                                                                                                                              aaaCAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTc                                <  1:361211/251‑1 (MQ=255)                                                                                                                                                                                                                  aGGCGCAAGAAAGGCCAACCATTAAGGCTACAAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACg                            >  1:146467/1‑251 (MQ=255)                                                                                                                                                                                                                     cgcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTg                         >  1:337828/1‑251 (MQ=255)                                                                                                                                                                                                                      gcAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGt                        <  2:337828/251‑1 (MQ=255)                                                                                                                                                                                                                         aGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCa                     <  1:196299/251‑1 (MQ=255)                                                                                                                                                                                                                           aaaGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGTCTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGGTCGAAGAGATAGAACTTTCGACGCTGTCCAtt                   <  1:267684/251‑1 (MQ=255)                                                                                                                                                                                                                              ggCCAACCATTAAGGCGACCAGCTCTCTTGGGTAACGTAATGGCTGAATTTGAATAAATCGGTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGTCTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCTTTccc                >  2:84465/1‑251 (MQ=255)                                                                                                                                                                                                                                 cAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTcacca                                                                                                                                       >  1:320907/1‑129 (MQ=255)                                                                                                                                                                                                                                 cAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGGGCTGGACTTCAGTcacca                                                                                                                                       <  2:320907/129‑1 (MQ=255)                                                                                                                                                                                                                                            gCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAg  <  1:282959/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  GGCCATCGCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAG  >  NC_000913/4602265‑4602750

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