Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,084,330	C→T	A42T (GCA→ACA)	tsuB ←	putative thiosulfate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,084,330	0	C	T	97.8%	150.4 / NA	46	A42T (GCA→ACA)	tsuB	putative thiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/21);  minor base G (0/1);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATACAAGGCCAA  >  NC_000913/2084092‑2084567                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              aaTTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCACTGTTAAAAGCACTTTGCCCGTACGCGGAAGAGATTGTCTAATTTGAAGACTCCGCCGCTGGTACGGCGGAGGACAAAGCTGCGTCAGGCTTTTTGAACGGTGATGCTCAACCATGCATCGCCAATCTGCTGATAATCGGTGATGGCATTTCCTTCTTCTGCCGCCCAATTTGGGATTGCTTCCGTTGTCTGGGAGCAATc                                                                                                                                                                                                                                   >  1:335444/1‑251 (MQ=255)                 acaacaATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCCCTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTaa                                                                                                                                                                                                                    <  1:184271/250‑1 (MQ=255)                 acaacaATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAAc                                                                                                                                                                                                                   >  1:123450/1‑251 (MQ=255)                      aaTATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCATTACGCGGAAGAGATTTTCTAATTTGAAGGCTCCGCCGCTGGTACTGCGGAGGACAATGTTGATTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTTATAATCGGTGATGGCATGTCTTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATc                                                                                                                                                                                                              <  1:267229/251‑1 (MQ=255)                       atatGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATcg                                                                                                                                                                                                             <  1:297642/251‑1 (MQ=255)                        tatGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAAATCAATCACTAACTCATcgc                                                                                                                                                                                                            >  1:560820/1‑251 (MQ=255)                            ccGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCt                                                                                                                                                                                                        >  1:329507/1‑251 (MQ=255)                                  aCAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCCCTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGCCCAt                                                                                                                                                                                                  >  1:52834/1‑251 (MQ=255)                                           aaGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGGGTCAGGCTTTTTGCACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTTTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGcc                                                                                                                                                                                           >  2:278393/1‑249 (MQ=255)                                             gTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCagag                                                                                                                                                                                       <  1:542603/251‑1 (MQ=255)                                                       tAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGCGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTTCCAGAGCTGCTTTTGc                                                                                                                                                                             >  1:477238/1‑251 (MQ=255)                                                               cTTTGGCCGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATg                                                                                                                                                                      >  1:158711/1‑250 (MQ=255)                                                                    gCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTGCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCGCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAAGCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGa                                                                                                                                                                <  1:505430/251‑1 (MQ=255)                                                                     cAGTACGAGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCCTTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGa                                                                                                                                                                >  2:582225/1‑250 (MQ=255)                                                                                agagATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGAcc                                                                                                                                                    <  2:304057/251‑1 (MQ=255)                                                                                                          ccgcTGGTACGGCGGAGGACAATGTTGAGTCAGGTTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCCGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCATTTGTCTGCGTGCAATTATACTCAATCTCTATCTCATCGCCGCTGACCATTTCTGCCCGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCTCCACATCCAGCTTTTTAAt                                                                                                                          <  1:517481/251‑1 (MQ=255)                                                                                                            gcTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATcc                                                                                                                        >  1:517299/1‑250 (MQ=255)                                                                                                               ggTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCAc                                                                                                                       <  1:278393/249‑1 (MQ=255)                                                                                                                gTACGGCGGAGGACAAGGTTGAGGCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCAt                                                                                                                    <  2:379725/251‑1 (MQ=255)                                                                                                                   cggcggAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCtt                                                                                                                 >  2:118010/1‑251 (MQ=255)                                                                                                                       ggaggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTac                                                                                                              >  1:497409/1‑250 (MQ=255)                                                                                                                        gaggaCAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTgcac                                                                                                            >  2:294954/1‑251 (MQ=255)                                                                                                                             cAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCACACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACCTCCCGCTTTTTCATCACCATCTTTACACcttac                                                                                                       >  2:255832/1‑251 (MQ=255)                                                                                                                                        cAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGcc                                                                                             <  2:329507/250‑1 (MQ=255)                                                                                                                                                              tCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCtt                                                                      >  1:98844/1‑251 (MQ=255)                                                                                                                                                               ccACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCtt                                                                      >  1:22525/1‑250 (MQ=255)                                                                                                                                                                   gCTGCATCGCCAATCTTCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACAGCCAGCTTTGTAATCACCATCTTTACACCTTACTTAATTTTCCGCAGCAGTTGCCAGTcgcg                                                                          <  2:275337/242‑1 (MQ=255)                                                                                                                                                                   gCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTcgcg                                                                          >  1:275337/1‑242 (MQ=255)                                                                                                                                                                      gcagcGCCAATCTGCTGAGAATCGGTGATGGCATGGCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAAGCAAACTCAATCACTAAATCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACAGCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCCGCAGTTGCCAGTCGCGCCTAACGCTg                                                                <  2:112012/245‑1 (MQ=255)                                                                                                                                                                          cGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTc                                                           <  2:98844/250‑1 (MQ=255)                                                                                                                                                                                 ctgctgATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATac                                                   >  2:318005/1‑251 (MQ=255)                                                                                                                                                                                    ctgATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATacac                                                 >  2:105728/1‑250 (MQ=255)                                                                                                                                                                                      gATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATAcaca                                                <  2:158711/249‑1 (MQ=255)                                                                                                                                                                                      gATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACaa                                               <  1:577706/250‑1 (MQ=255)                                                                                                                                                                                        tAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGcc                                            >  1:7204/1‑251 (MQ=255)                                                                                                                                                                                        tAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGcc                                            >  1:12785/1‑251 (MQ=255)                                                                                                                                                                                              gTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGaa                                                                                                                                                               >  1:586797/1‑130 (MQ=255)                                                                                                                                                                                              gTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACg                                      <  2:457259/251‑1 (MQ=255)                                                                                                                                                                                              gTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACGCAGCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGaa                                                                                                                                                               <  2:586797/130‑1 (MQ=255)                                                                                                                                                                                                gATGGCATGGCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTGTAAGCACCATCTTTACACCGTACTTAATTTGCCGCAGCAGTTTCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCa                                    <  1:251251/251‑1 (MQ=255)                                                                                                                                                                                                             tcttctGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTcc                         <  1:118010/249‑1 (MQ=255)                                                                                                                                                                                                                 ctGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAAt                   >  2:8304/1‑251 (MQ=255)                                                                                                                                                                                                                         cACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATa           <  2:149791/251‑1 (MQ=255)                                                                                                                                                                                                                          aCTGTGGGATCGCTTCAGTTGTCTGGGTGCAATCCAACTCACTCAATAACTCATCGCCGCTGACCATTTCTGCCAGCGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATAc          >  2:555343/1‑251 (MQ=255)                                                                                                                                                                                                                             gtgGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATATATACAAg       <  2:382511/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTTGCTTCAGTTGGCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCGTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATACAAGGCCaa  <  1:255832/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              AATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAATGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATACAAGGCCAA  >  NC_000913/2084092‑2084567

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