Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,298,401	G→A	L322L (CTC→CTT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,298,401	0	G	A	93.3%	116.1 / ‑0.3	46	L322L (CTC→CTT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (25/17);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

GCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAAC‑GATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAA  >  NC_000913/4298161‑4298637                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              gCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCgg                                                                                                                                                                                                                                     >  1:151927/1‑251 (MQ=255)      tAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGc                                                                                                                                                                                                                                >  1:157198/1‑251 (MQ=255)          gTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATcg                                                                                                                                                                                                                            <  1:103107/251‑1 (MQ=255)               cgtcacgcacTTAGCAATACGCTGCGCGGTGCGGAAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCGCGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGGACTGGTATCCCGGATCCGTAACCGTCAGAGCGCCCATATCGCAGGAACCCTGTACGTTGTTCTGACCACGAACAGGGTTATAACCCGCATGCGGCTTACCAAGGTTACCGGTccacctccctagg                                                                                                                                                                                                                       >  2:118663/4‑239 (MQ=255)               cgGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAgg                                                                                                                                                                                                                       <  1:137099/251‑1 (MQ=255)                 gcacgcacTTCGCCATGCGCTGCGCGGTGCGCAAGCTCGCTGACGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTTGCGAACTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCCCGAACCGGGTTCACACCCGCCTGCGGATTCCCACGGGTACCGGTCATCACCGCGCGGCt                                                                                                                                                                                                                     >  2:371643/1‑251 (MQ=255)                  cacgcacTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTg                                                                                                                                                                                                                    >  1:215604/1‑251 (MQ=255)                           caccaggcgctgcgcGTGGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGATCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAgaga                                                                                                                                                                                                            <  2:122794/244‑1 (MQ=255)                             ccATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCTCGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACg                                                                                                                                                                                                         <  1:32261/251‑1 (MQ=255)                               aTGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGc                                                                                                                                                                                                        >  1:181283/1‑250 (MQ=255)                                          ggTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACg                                                                                                                                                                                            >  1:296474/1‑251 (MQ=255)                                                gCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTgg                                                                                                                                                                                      <  1:47990/251‑1 (MQ=255)                                                  cgCTCGGTGGTGCTATAGCAGCTATGCGGTGGCAGACTTTCCACGCCCCAGCCTTCTGCGAATTTCTCGTGGTTAGCCGGATCTTTCGCGTACTGGTACCCCGGATACGGAGCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACTCGCAGGCGGCGTACCAAGGTTTCCGGTCAGCATCGCGAGGCTGTTCAGTGATCGCCCGGTTTCCACGCCCTGGTa                                                                                                                                                                                    <  2:57287/250‑1 (MQ=255)                                                      cactgttgcgtccgccggtattgtttGGCAGGCGGTCCACGCCCCGTGTTTTGGCGAATTTCTCGCGGTTAGCCTGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCTCCCTGAACGTTGTTCGTACCTCGAACCGGGTTAACTCCCGCATGCGGCTTATCTAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGATAACGCTCGTTTTCCACGCCCTGGTAGaa                                                                                                                                                                                 <  2:222578/225‑1 (MQ=255)                                                      cGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGaa                                                                                                                                                                                 <  1:382209/250‑1 (MQ=255)                                                          gATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTggg                                                                                                                                                                            >  1:256392/1‑251 (MQ=255)                                                           aTGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGTGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGt                                                                                                                                                                           <  1:353499/251‑1 (MQ=255)                                                                aTAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGCCCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGCACGCACGGATTCCACGCCCTGGTAGAACTGGGTTaaac                                                                                                                                                                      >  2:225806/1‑248 (MQ=255)                                                                    ccGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCAt                                                                                                                                                                  >  1:281953/1‑251 (MQ=255)                                                                        gtggcgcGGGCAGGCTTTCCATGTCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACCTGTATCCCGGATACGGACACGGCAGCGCGCCCATAACGCAGGCACCCTTAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGccc                                                                                                                                                              <  2:178863/248‑1 (MQ=255)                                                                             gcTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGcccc                                                                                                                                                             <  1:356720/247‑1 (MQ=255)                                                                             gcTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATg                                                                                                                                                                 <  1:314056/243‑1 (MQ=255)                                                                             gcTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGGTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTCAACCCATg                                                                                                                                                                 >  2:314056/1‑243 (MQ=255)                                                                             gcTGGCAGGCTTTACACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATATTTCACGTACTGGTAACCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTAAACGTTGTTCTGACCACGACCCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGACAGCATCGCGAGGCTGGTCAGAGAACGCCCGGTTTCCACGCCCTGGTAGAACTGGGTTacacacctgaccc                                                                                                                                                             >  2:356720/1‑238 (MQ=255)                                                                                    ggCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGCACTGGCTTACACCCCTCCCCCACAgggcggcg                                                                                                                                                  >  2:116269/1‑245 (MQ=255)                                                                                        ttCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATggcggc                                                                                                                                                <  1:225806/249‑1 (MQ=255)                                                                                        gtCCACGCCCCAGGCTTTGGCGAATTTCTCGCTGTTAGCCGGATCTTTCACGTACTAGTATCCCGTTTACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTTAACGTTGTTCTGACCACGAACCGTGTTAACTCCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCTCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATggcggct                                                                                                                                               <  1:243420/249‑2 (MQ=255)                                                                                           cACGCCCCATGCTTTGGCGAAGTGCTAGTGTGTAGCCGGATCTTTCGCGTACTGGTATCCCGTATACGTAGCCTGCTGCGCGCCCATAGCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCAGGCGGCGTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACTCCCATGCCCCACAGGATGGCGGCGCttt                                                                                                                                           <  2:285998/251‑1 (MQ=255)                                                                                           cACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGATCTGACCACGAACCGGGTTAACACCCGCATGCGGCATACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCttt                                                                                                                                           >  1:41001/1‑251 (MQ=255)                                                                                              gCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCg                                                                                                                                        >  1:263960/1‑251 (MQ=255)                                                                                                                        ggTTAACCGAATCTTTCACGTCCTGGTATCCCGGATACGTCTCAGGCAGCGCGCCCATATCGCAAGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGACGCTGGTCAGAGAACGCAAGGTTTCCACGCCCTGGTAGAACTGGGTTACAACCATGCCCAACACGATGGCGGCGCTTTTCGCCTGGGCATACATCCGAgcctcacg                                                                                                               >  2:70744/1‑245 (MQ=255)                                                                                                                                  aTCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCAc                                                                                                    >  1:372798/1‑251 (MQ=255)                                                                                                                                      ttCACGTACTGGTCTCCCGGATCCGTATCCGGCAGCGCGCCCATATCGCCGGTCCCCTGAACGTTGGTCTGACCACGAACCGGGGTAACTCCCGCATTCGGCTTACCAATGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTTGGTTACACCCAGGCCCCACAGGATGTCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCa                                                                                                     <  2:215604/246‑1 (MQ=255)                                                                                                                                                   aTCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTg                                                                                   >  1:112803/1‑251 (MQ=255)                                                                                                                                                         gATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCt                                                                             >  1:70402/1‑251 (MQ=255)                                                                                                                                                          aTACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGACCCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGCCAGAGAACGCACGGTTTCAACGCCCTGGGAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCACCCTGGCGCTGACGCCGGGCATATCtt                                                                            >  2:145862/1‑251 (MQ=255)                                                                                                                                                                            gcCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCgtgt                                                          <  1:283215/251‑1 (MQ=255)                                                                                                                                                                                  atCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCtt                                                    >  1:117184/1‑251 (MQ=255)                                                                                                                                                                                               tGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTCGAACTGGGTTACACCCATGCCCCACAGGATGGCGACGCTTTTCGCCTGGGCATACCTCCTTGCCGCCTGACGAATCTCACTGGCGCTGTCGCCGGTGATATCTTCAACCGACTCCGGCGTGATGCCTTCACC‑GATTTTagc                                      >  2:311235/1‑249 (MQ=255)                                                                                                                                                                                               ggaaagttgttCTGACCCCGAACCGGGCTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGTCTCCGGCGTGTAGCCTTCAAC‑GATTTTACg                                      <  2:256392/246‑1 (MQ=255)                                                                                                                                                                                                              ccACTAACCGGGTTAACACCCGCATGCGTCTTACCAAGGTTACCGGTCAGCATCGCGAGGCAGGTCAGAGAACGCACGGTTTCCACGCCCTGGGAGAACTGGGTTACAACCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCCCACTGGCGCTGCAGCCGGTGATATCTTCAACCGACTCCGGAGTGGAGCCTT‑AACAGATTTTACGATACTTTTCaaatcc                       >  2:220614/1‑248 (MQ=255)                                                                                                                                                                                                               cACGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAAC‑GATTTTACGATACTCTTCAAAGCCt                      >  1:120361/1‑251 (MQ=255)                                                                                                                                                                                                                aCGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATggcggc                                                                                                                                                >  1:299982/1‑129 (MQ=255)                                                                                                                                                                                                                aCGAACCGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCAACGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATggcggc                                                                                                                                                <  2:299982/129‑1 (MQ=255)                                                                                                                                                                                                                      cGGGTTAACACCCGCATGCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAAC‑GATTTTACGATACTCTTCAAAGCCTTCTGTAc               >  1:98354/1‑251 (MQ=255)                                                                                                                                                                                                                                   gcaggCGGCTTACCAAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAAC‑GATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGaa  <  1:193149/247‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              GCCCATAATGTACGCGGCACGCACTTCGCCATGCGCTGCGCGGTGCGGCAGCTCGCTGATGCGATAGCCGGTATGCGCTGGCAGGCTTTCCACGCCCCAGGCTTTGGCGAATTTCTCGCGGTTAGCCGGATCTTTCACGTACTGGTATCCCGGATACGTATCCGGCAGCGCGCCCATATCGCAGGCACCCTGAACGTTGTTCTGACCACGAACCGGGTTAACACCCGCATGCGGCTTACCGAGGTTACCGGTCAGCATCGCGAGGCTGGTCAGAGAACGCACGGTTTCCACGCCCTGGTAGAACTGGGTTACACCCATGCCCCACAGGATGGCGGCGCTTTTCGCCTGGGCATACATCCGTGCCGCCTGACGAATCTCACTGGCGCTGACGCCGGTGATATCTTCAACCGACTCCGGCGTGTAGCCTTCAAC‑GATTTTACGATACTCTTCAAAGCCTTCTGTACGTGAAGCGACGAA  >  NC_000913/4298161‑4298637

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