Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,297,782	G→A	P529S (CCT→TCT)	fdhF ←	formate dehydrogenase H

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,297,782	0	G	A	96.2%	136.4 / NA	54	P529S (CCT→TCT)	fdhF	formate dehydrogenase H
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (28/22);  minor base C (1/1);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCAC  >  NC_000913/4297538‑4298026                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     aTTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACtg                                                                                                                                                                                                                                                <  1:209262/251‑1 (MQ=255)      tgCGTAGCCAGGTTAATCTGCAAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATAGGGTACTCGTCGGTGCGATTGTCGCTTGGCGCTACCCCATCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAACCAGATAAGACGACCCCTGATCGGCATCTTAAGCATGCCGGCAAGACCACtgaata                                                                                                                                                                                                                                            >  2:449828/1‑249 (MQ=255)        cGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGAGAAGAACTGCGCCAGACAGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCACTGATCGGCATCTGAAGTATCGCGGCAAGACCACtgaatga                                                                                                                                                                                                                                           >  1:609461/1‑249 (MQ=255)          cggctcgtgtcaccaccaagtcccgccAGTTCCGCACAGTTACCGGTCATCGAACGGCATGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTGTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGGCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGcc                                                                                                                                                                                                                                       <  2:24466/229‑1 (MQ=255)                 gTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGGCCCCTGATCGGCATCTGAAGTATCGCGGCACGACCACTGAATGAAGcc                                                                                                                                                                                                                                       >  1:24466/1‑244 (MQ=255)                    catcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGCACGGCAAGAGTAGTGCCCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCACTTCGccc                                                                                                                                                                                                                             >  2:156606/1‑251 (MQ=255)                    catcaGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGCACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCACGACCACTGAATGAAGCCCAGTTCGccc                                                                                                                                                                                                                             >  1:494942/1‑251 (MQ=255)                          ccagcGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGTCCACTGAATGAAGCCCAGTTCGCCCATTTtc                                                                                                                                                                                                                       <  1:574315/251‑1 (MQ=255)                            agcGCCGCCAGTGCCGCACAGTGACCGGTCAGCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCACCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGACGCAGGGGAAGAACTGCGCCAGACCGTTCGTGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGACCAGTTCGCCCATTTtctc                                                                                                                                                                                                                     <  2:378422/251‑1 (MQ=255)                             gcgcCGCCAGTGCCGCACAGGTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCAGCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCg                                                                                                                                                                                                                    <  2:599506/251‑1 (MQ=255)                               gccgccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGATTGAAGCCCAGTTCGCCCATTTTCTCGt                                                                                                                                                                                                                   >  1:24056/1‑250 (MQ=255)                               gccgccAGTGCCGCACAGTTACCGATCATCGACCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTCCCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGCGAAGAACTGCGCCAGACCGTTCGGGGTCTCAAACTTCTCTTTAAACAGATAAGAAGCCCCCTGATCGGCCTCTGATGTATCGCGGCAAGAACACTGAATGAAGCCCAGTTCGCCCCTTTTCTCGtc                                                                                                                                                                                                                  >  2:576852/1‑250 (MQ=255)                                   ccAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTc                                                                                                                                                                                                              <  1:642205/251‑1 (MQ=255)                                                    ccGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCgg                                                                                                                                                                                             <  1:206931/251‑1 (MQ=255)                                                     cGGTCAGCGAACGGCAATAGTAGGGACCAACATCCCTAGCCGTTGACAGTACCAGCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTGCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAATTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCggg                                                                                                                                                                                            <  2:559895/251‑1 (MQ=255)                                                          aTCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAg                                                                                                                                                                                       <  1:48523/251‑1 (MQ=255)                                                               aCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGAc                                                                                                                                                                                  <  1:513791/251‑1 (MQ=255)                                                               aCGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGAc                                                                                                                                                                                  <  2:113618/251‑1 (MQ=255)                                                                cGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAACCTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCACGCCCCCTGAATGAAGCCAAGTTCGCCCCTTTTCTCGTACGTCGCACCGTAGAAATCCGGGcccactttacc                                                                                                                                                                                 >  2:326786/1‑241 (MQ=255)                                                                    aaGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAACCTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGCATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAAc                                                                                                                                                                             >  1:569892/1‑251 (MQ=255)                                                                          agtGACCAACTTCCCGCCCCGTCGACAGTACCATCGGGTACTCGTCGGGCAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCTTTCGGGGTATCAAACTTCTCTTTAAACTGATAAGAAGTCCCCTGATCGGCAACTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCTATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATc                                                                                                                                                                        <  2:434760/250‑1 (MQ=255)                                                                                      cACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACGGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTgg                                                                                                                                                           <  1:236944/251‑1 (MQ=255)                                                                                         gCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCTCAGATGACGCAACTCATCCCAGATCTCCTGGgtg                                                                                                                                                        >  1:556285/1‑251 (MQ=255)                                                                                                ggACAGTTCCATCGGTGACTTGTCGGTGAGTTTGTCGAGTGGCGCTACCCAGTTGCAGGTGAAGGACTGCGCCAGACCGTTCGGGGTATCCAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGTGGCAAGACCACTGAATGAAGTCCAGTTCGCCCAGTTTCTCGTACGTCTCACCGTAGAAATCCGGGCTCAGATTACGCAACTCATCCCATATCTCCTGGGtgttgttt                                                                                                                                                  <  2:423682/249‑2 (MQ=255)                                                                                                 gacagtacgttcgggttcgtggcgGTGAGTTTGGCGATTGGCGCGTCCCAGTCGCTGGTGAAGAACTGCGCCAGACCGTTCGTGTTATCAAACTTCTCTTTAAACAGATAAGAAGGCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCTCCGTAGAAATCCTGGCACAGATGGCGCAACTCATCCCTGATCTCCTGGGTGTTGTTGTa                                                                                                                                                <  2:623941/229‑1 (MQ=255)                                                                                                     gTACCATCGGGTAATCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGAAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGc                                                                                                                                            >  1:504666/1‑251 (MQ=255)                                                                                                       aCCATCGGGTACTCGTCGGTGAGTTTGTCGATTGACGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGACCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGCAAGTCGACCCGTAGAAATCAGGGCACAGAAGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGAGCAt                                                                                                                                          >  2:596581/1‑251 (MQ=255)                                                                                                         cATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCg                                                                                                                                        >  1:557880/1‑251 (MQ=255)                                                                                                            cGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCACACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCAGCAAGACCACTGCATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTAATCCCAGATCTCCTGGTAGTTGTTGTAGTGCATCGGa                                                                                                                                      >  2:32377/1‑250 (MQ=255)                                                                                                               gTACTCGACGGTGAGTTTGTCGATTGTCGTGACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTGCTCTTTAAACAGATAAGAAGTCGCCTGATCGGCCTCTGAAGTATCGCGGCAAGACCACTGAATGAAGTCCAGTTCGCCCATTTTCTCGTAAGTCTCCCTGTAGAAATCCGGGCACAGATGACGCTACTCATCCCAGATCTCCTGGGTGTTTTTGTAGTGCATCGTATAAc                                                                                                                                  <  2:297257/251‑1 (MQ=255)                                                                                                                        ggTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACggg                                                                                                                         >  1:213161/1‑251 (MQ=255)                                                                                                                                 gTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGTTGGCGAtt                                                                                                                 <  2:337989/250‑1 (MQ=255)                                                                                                                                              aCCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTg                                                                                                   >  1:680946/1‑251 (MQ=255)                                                                                                                                                         ggTGTAGAACTGCGCCAGACCGTTCGGGGCATCAACCTTCTCTTTAAACAGATAGGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGATTGAAGCCCCGTTCGCCCATTTTCTCGTAAGTCGCACCGAAGAACTCCGTGCACAGATGACGCACCTCATCCCAGATCTCCTGGGTGTTGATGTAGAGCATCAGATAGCCCATACGGGTGGAGCTTTCACTGATGATTGTCCAGTCCGttt                                                                                        >  2:588726/1‑251 (MQ=255)                                                                                                                                                          gTGACGAACTGCGCCAGACCGTCCGGGGTATCAAACTTCTCTTTAAACATATAAGAAGTCCCCTGAGCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCAACGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGTTGTTGTTGTAGTGCATCGGATAACCCATTCGGGTGGCGATTTCACTGATGATTTGCCAGTCCTtttt                                                                                       <  2:541670/251‑1 (MQ=255)                                                                                                                                                               aaaCTGCGCCAGACCGTTTTGGGTATCAAACTTATCTTTATACAGAGAAGAAGTCCCCTGATCGGCATCTGAAGTAGCGCGGCAAGACCACTTAATTAATCACAGTTCGCCTATTTTCTCGTATGTCGTATCGTAGATAGCCGGGCACAGATGACGTAACTCATTCCAGATCTCCTGGTTGTTGTTGTAGTGCATCCGATAATCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGAt                                                                                  <  2:517085/250‑1 (MQ=255)                                                                                                                                                                aaCTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTc                                                                                      >  1:52423/1‑246 (MQ=255)                                                                                                                                                                 aCTGCGCTAGACCGTTCTGGGTATCAAACTTCTCTTGAAACAGAGAAGAAGTCCACTGAGCGGCATCTGAAGTATCGCGGCAAGAACACTGAATGAAGCCTAGTTCGCCCATTTGCTCGTAAGGCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGAGTTGCCAGTCCGTTTTCAGagcc                                                                                <  2:494942/251‑4 (MQ=255)                                                                                                                                                                  cTGCGCCAGACCGTTCGGGGTATCAAACTGCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCAGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGGAGTGCATCGGCTAACCCATACGGGGGGCGATTTCACTGATGTATTGCCAGGCCGtttccaagtcc                                                                                >  2:152220/1‑242 (MQ=255)                                                                                                                                                                             cGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTTAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCa                                                                    <  1:204735/251‑1 (MQ=255)                                                                                                                                                                                  ggggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGc                                                               >  2:571700/1‑251 (MQ=255)                                                                                                                                                                                    ggTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCt                                                             >  2:519149/1‑251 (MQ=255)                                                                                                                                                                                       atcaaaactctcTTTAAACAGAAAAGAAGTCCCCTGAGCGGCATCTGAAGTATCGCGGCAAGTCCTCTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACTCAACTCACCCCAGATCTCCTTGGTTTTTGTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCTGTACTTTTTCAGATCCCATTTCGGTTCAACCGCCTTga                                                          <  2:95712/243‑1 (MQ=255)                                                                                                                                                                                                                 aaggccccggaccGGCATCGGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGGTGTTGTAGTGCTTCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGcc                                <  1:576852/239‑1 (MQ=255)                                                                                                                                                                                                                  agTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGTTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCg                               <  1:24191/251‑1 (MQ=255)                                                                                                                                                                                                                        cTGATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTa                          >  1:392992/1‑250 (MQ=255)                                                                                                                                                                                                                          gATCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGACACGCTGGAAGCCACGGTCAGCCGCAGTaaa                        >  1:280947/1‑250 (MQ=255)                                                                                                                                                                                                                           aTCGGCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAac                       >  1:509520/1‑250 (MQ=255)                                                                                                                                                                                                                               gCATCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGcc                  >  2:194010/1‑251 (MQ=255)                                                                                                                                                                                                                                 aTCTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCtt                >  1:277914/1‑251 (MQ=255)                                                                                                                                                                                                                                   cTGAAGTATCGCGGCAAGCCCACTGAATGAAGCCCATTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGACCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGaa                                             <  2:67127/220‑1 (MQ=255)                                                                                                                                                                                                                                   cTGAAGTATCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGaa                                             >  1:67127/1‑220 (MQ=255)                                                                                                                                                                                                                                          aTCGCGGCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGc       >  1:528731/1‑251 (MQ=255)                                                                                                                                                                                                                                               ggCAAGACCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCac  >  1:396575/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     ATTTGTGCGTAGCCAGGTTCATCAGCCAGCGCCGCCAGTGCCGCACAGTTACCGGTCATCGAACGGCAAGAGTAGTGACCAACTTCACGCACCGTTGACAGTACCATCGGGTACTCGTCGGTGAGTTTGTCGATTGGCGCTACCCAGTCGCAGGTGAAGAACTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGACGCAACTCATCCCAGATCTCCTGGGTGTTGTTGTAGTGCATCGGATAACCCATACGGGTGGCGATTTCACTGATGATTTGCCAGTCCGTTTTCAGATCCCATTTCGGTTCAACCGCCTTGAAGAAACGCTGGAAGCCACGGTCAGCCGCAGTAAACACGCCTTCATGCTCGCCCCAC  >  NC_000913/4297538‑4298026

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