Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,160,157	G→A	L212L (CTG→TTG)	sthA ←	soluble pyridine nucleotide transhydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,160,157	0	G	A	98.4%	199.9 / NA	64	L212L (CTG→TTG)	sthA	soluble pyridine nucleotide transhydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (32/30);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCC  >  NC_000913/4159912‑4160401                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGc                                                                                                                                                                                                                                                 <  1:428510/251‑1 (MQ=255)  aGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCg                                                                                                                                                                                                                                                <  1:316634/251‑1 (MQ=255)      gTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATcg                                                                                                                                                                                                                                            <  1:475315/251‑1 (MQ=255)      gTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCGTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATcg                                                                                                                                                                                                                                            <  2:399987/251‑1 (MQ=255)           cgGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGgt                                                                                                                                                                                                                                       <  2:464020/251‑1 (MQ=255)            ggCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGgtg                                                                                                                                                                                                                                      <  1:253027/251‑1 (MQ=255)                  cAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGAt                                                                                                                                                                                                                                 <  1:464284/250‑1 (MQ=255)                  cAGTTTATAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATc                                                                                                                                                                                                                                >  1:467319/1‑251 (MQ=255)                           gCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACt                                                                                                                                                                                                                       >  2:458337/1‑251 (MQ=255)                            cccAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCTGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACtt                                                                                                                                                                                                                      >  1:261334/1‑251 (MQ=255)                                   ttCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATc                                                                                                                                                                                                               >  1:227369/1‑251 (MQ=255)                                           cGCCAACGAATCAGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCCCTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTTACATCTCTTGATCGAGACATTCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACcgc                                                                                                                                                                                                       >  2:479242/1‑251 (MQ=255)                                                 cGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGa                                                                                                                                                                                                    >  1:224380/1‑248 (MQ=255)                                                    aTCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCAGAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATTTCTTGATCGAGAAATGCCAACAGGCGAGCGCGTGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCg                                                                                                                                                                                              <  1:222515/251‑1 (MQ=255)                                                         tATTACCGGAGCGACAGTTGGCATAGAGCAGGCAGTCTGCTTTCAGTTTTTTAACCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGCCGAATCACTACGCCACTGTTCCAGAAGTGCTAGGTGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTCCATCCATACCGCGGAAGagcgccgc                                                                                                                                                                                          >  2:560624/1‑243 (MQ=255)                                                          aTTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCat                                                                                                                                                                                        >  1:183785/1‑251 (MQ=255)                                                          aTTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCGTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCa                                                                                                                                                                                         <  1:554265/250‑1 (MQ=255)                                                                   gCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAg                                                                                                                                                                                <  2:286168/250‑1 (MQ=255)                                                                    cGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                                                                                                                                                                              <  1:149472/251‑1 (MQ=255)                                                                    cGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGcc                                                                                                                                                                              <  2:261334/251‑1 (MQ=255)                                                                      cccGTTGGCAGAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGAGCGCGGGTGTTGATCAGATCCACTTTTACATCCCTACCGCGGAAGCTCGACGCATATTCACAGCCg                                                                                                                                                                             >  2:482038/2‑250 (MQ=255)                                                                                        gcagTCAGCTTTCAGTTTTTTAACCGTCTTCAGATGCATGATCACACAATCGTCACAGCCTTCGATCTTCTCGTAATCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAAAAGGCGATCGCGGGTGTTGATCAGATCAACTTTTACATCCATACCGCGGAAGATCGCCGCATATTCACAGCCGATCACTCCAGCAcaat                                                                                                                                                             >  2:534694/1‑245 (MQ=255)                                                                                          agTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATa                                                                                                                                                        <  2:224380/251‑1 (MQ=255)                                                                                              cgCTTTCAGTTTTTTACCCGACTTCAGATGAATGATCACACCCTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGCCATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGACGATAGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGt                                                                                                                                                     >  2:504417/2‑250 (MQ=255)                                                                                                         ttttACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCgg                                                                                                                                          <  1:286250/250‑1 (MQ=255)                                                                                                            tACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCg                                                                                                                                      >  1:489798/1‑251 (MQ=255)                                                                                                                 gACTTCAGATGCCTGATCACACCTTCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTTATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCGCTTTTATATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCATCATAGATAAGTACATGGCGCGGTTCgtgttg                                                                                                                                 <  1:334504/251‑4 (MQ=255)                                                                                                                 gACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTGCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtggtg                                                                                                                                 <  1:334486/251‑1 (MQ=255)                                                                                                                 gACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCGTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCgtggtg                                                                                                                                 <  2:24313/251‑1 (MQ=255)                                                                                                                      cataTGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCGTGATCGAGAAAGGCCAACAGGCGAGCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCTGAAGATCGACGCATATTCACAGCCGATCACTGCAGCACCATAGATAGGTACATTGCGGGGTTCGTAGTGCATg                                                                                                                             <  2:380224/247‑1 (MQ=255)                                                                                                                        gATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTg                                                                                                                          >  2:199429/1‑251 (MQ=255)                                                                                                                             aTGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAAt                                                                                                                     >  1:279967/1‑251 (MQ=255)                                                                                                                             aTGATAACACAATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCCCTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAAt                                                                                                                     >  2:494728/1‑251 (MQ=255)                                                                                                                                    cacCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCg                                                                                                              >  1:321742/1‑251 (MQ=255)                                                                                                                                            tCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGt                                                                                                       >  1:181311/1‑250 (MQ=255)                                                                                                                                                 gCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCAATGTTCAAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCCTAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGc                                                                                                 >  1:544177/1‑251 (MQ=255)                                                                                                                                                      cGATCGTCTCGTTCTCTTCGTTGTGACGAATCACTACGCCACTGTGCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGAGCGCGGGTGTTGAGCAGATCCACTGTTACATCCATACCGCGGAAGATCGACGCATATTCACAGACGATCACTCCAGCACCATAGATAAGTACATGGTGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGa                                                                                            <  1:479242/251‑1 (MQ=255)                                                                                                                                                                  aCTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATc                                                                                 >  2:202187/1‑250 (MQ=255)                                                                                                                                                                       tcgtgtttaCGAATCACTACGCCACTGTTCCAGAAGAGATAGGAGAGAGAATCTGACATCTCTTTATCGAGAAATGCCAACAGGCGAGCTCTGGTGTTGTGCAGTTCCACCTTTACATCCATACCGCGGAAGATCGACGCATATGCACAGCCGATCGCTCCAGCACCATAGATAAGTACATGGCGCGTTTTGTGGTGCATTCTGTGAATTGATTCGCTGTCGTAAATGCGTGGAGGTGGGAAATTAACAt                                                                            <  2:525991/242‑1 (MQ=255)                                                                                                                                                                       tCGTTGGGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGAGAAGTACATGGCGCGGTTCGTGGTGCAGGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAAGCAACAt                                                                            <  2:82910/250‑1 (MQ=255)                                                                                                                                                                         gttgttcagAATCACTACACCCCTGTTCCAGAAGAGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTg                                                                         >  2:244783/9‑251 (MQ=255)                                                                                                                                                                         gTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATGCACAGCCGATCACTCCAGCACCAGAGATAAGGACATGGCGCTGTTCGTGGTGCATGCGGAGAATTGAGTCGCTGTCGTAAAGGCGTGGAGTGGTTAAAGCAACATCTg                                                                         <  2:262199/251‑1 (MQ=255)                                                                                                                                                                                cGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCt                                                                                                            <  1:304267/209‑1 (MQ=255)                                                                                                                                                                                cGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCt                                                                                                            >  2:304267/1‑209 (MQ=255)                                                                                                                                                                                       cTACGCCACTGTTCAAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGAACATGGCGCGGTTCGTGGAGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGAGAAATCAACATCTGTTGGATGATATgg                                                            >  2:395576/1‑250 (MQ=255)                                                                                                                                                                                                    ccAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCACCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTAAATCCATACCGCGGACGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGTTGATCTGGACGAGAGCCGCAgg                                              >  2:381163/1‑251 (MQ=255)                                                                                                                                                                                                     cAGAAGTGATAGGAGAGAGAATCTGACAGCTCGTGATCGAGAAAGGCCAACAGGCGATCGCTGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGTATAGTCACAGCCGATCACTCCAGCACCATAGATAATTACATGGCTCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGc                                             <  1:535693/251‑1 (MQ=255)                                                                                                                                                                                                            gATAGGAGAGAGAATCTGACCTCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCCTCGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTCAATGCGTGGATGGGTGAAATCAACATCTGTTGGCTGATATGGACGAGAGCCGCAGGCAATAAc                                       >  1:397149/1‑250 (MQ=255)                                                                                                                                                                                                                  agagagaATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAAtttt                                <  2:183785/251‑1 (MQ=255)                                                                                                                                                                                                                   gagagaATCTGACATCTCTTGCTCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGCCCCACTTTTACATCCATACCGCGGAAGATCGACGAATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGGGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAACttttt                               >  2:97519/1‑251 (MQ=255)                                                                                                                                                                                                                   gagagaATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAAttttt                               <  1:458337/251‑1 (MQ=255)                                                                                                                                                                                                                     gagaATCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCTGTGCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGGTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCa                             <  2:227369/251‑1 (MQ=255)                                                                                                                                                                                                                       gaagcTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACAGCTGTTGGATGATATGTACGAGAGCCGCAGGCAATAACAAATTTTTCAGc                           <  1:347896/247‑1 (MQ=255)                                                                                                                                                                                                                         aTCTGACATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGACGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCgg                         >  1:198644/1‑251 (MQ=255)                                                                                                                                                                                                                               cATCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAgtg                   >  2:19311/1‑251 (MQ=255)                                                                                                                                                                                                                                aTCTCTTGATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAgtgt                  >  1:161760/1‑251 (MQ=255)                                                                                                                                                                                                                                   tctTGATCGAGAAATGCCAACAGGCGATCGCGGGAGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTc               >  1:589381/1‑251 (MQ=255)                                                                                                                                                                                                                                       gATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGCATTGAGTCGCTGTCGTAAATg                                                                                                  >  2:214262/1‑164 (MQ=255)                                                                                                                                                                                                                                       gATCGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATg                                                                                                  <  1:214262/164‑1 (MQ=255)                                                                                                                                                                                                                                          cGAGAAATGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGCTGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACCCTg        >  2:314755/1‑251 (MQ=255)                                                                                                                                                                                                                                               aaTGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGGGCATGCTGAGAATTGAGTCGATGTCGTAAATGCGTGGATGGGTGAAATCAAAATCTGTTGGATGATATGGACGATCGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTc   >  1:566103/1‑251 (MQ=255)                                                                                                                                                                                                                                                aggCCAACAGGGGAACGCGGGTGTTGATCAGATCCTCTTTTACATCCATACCGCGGAAGAGCGACGAATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTcc  <  2:321742/249‑1 (MQ=255)                                                                                                                                                                                                                                                aTGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACAt                                                                            >  2:315184/1‑177 (MQ=255)                                                                                                                                                                                                                                                aTGCCAACAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACAt                                                                            <  1:315184/177‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAGCCTTCGATCTTCTCGTACTCTTCGTTGTGACGAATCACTACGCCACTGTTCCAGAAGTGATAGGAGAGAGAATCTGACATCTCTTGATCGAGAAATGCCAGCAGGCGATCGCGGGTGTTGATCAGATCCACTTTTACATCCATACCGCGGAAGATCGACGCATATTCACAGCCGATCACTCCAGCACCATAGATAAGTACATGGCGCGGTTCGTGGTGCATGCTGAGAATTGAGTCGCTGTCGTAAATGCGTGGATGGGTGAAATCAACATCTGTTGGATGATATGGACGAGAGCCGCAGGCAATAACAAATTTTTCAGCGGTTAGTGTTTCAACGCTGCCGTCC  >  NC_000913/4159912‑4160401

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