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	95.2%	132.3 / NA	44	A42T (GCA→ACA)	tsuB	putative thiosulfate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (23/17);  minor base G (0/2);  total (25/19)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

TGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATACAAGGCC  >  NC_000913/2084084‑2084565                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            tGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCCCTTTGGCAGTACGCGGAAGAGATTCTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCATCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTgg                                                                                                                                                                                                                                         >  1:409554/1‑251 (MQ=255)    tgatgaATTGCAGAAGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAACACTTTGGCAGTACGCGGAAGAGATTGGCTACTTTGAAGCCTCCGCCGCTGGTACGGCGGAGG‑CAACTGTTGAGTCAGGCTTTTTGAACGGAGATGCTCCACGCTGCATCGCCCACATGCTGAAAATCGTGGGTGGCATGTCCTTCTTCTCCCGCCCCCTGTGGGCTAGCTTCCGTTGACTGGGTg                                                                                                                                                                                                                                      >  2:226538/1‑251 (MQ=255)         aaTTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTTCAATc                                                                                                                                                                                                                                 <  1:191829/251‑1 (MQ=255)              cAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTGAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTAGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCAGGTCCTTCTTCTGCCGCCCACTGTGGGATGGCTTCAGTTGTCTGGGTTCAATTAAACt                                                                                                                                                                                                                            <  2:558076/251‑1 (MQ=255)                gATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGCCAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTc                                                                                                                                                                                                                           >  1:219558/1‑250 (MQ=255)                         acaacaATATGCCGCTGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAAc                                                                                                                                                                                                                 >  1:479311/1‑251 (MQ=255)                                 atGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCACTCAAACTCAATCACTACCTCATcgcc                                                                                                                                                                                                         >  1:375155/1‑251 (MQ=255)                                   gCCGATGACAAAGAGAAAGTGGCGCTGGTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCGCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATcgccg                                                                                                                                                                                                        <  2:375155/250‑1 (MQ=255)                                      gATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCGAATTTGAAGCCGCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGa                                                                                                                                                                                                    <  2:317001/251‑1 (MQ=255)                                                 gaAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCGCCACGCTGCATCGCCAATCTGCTGATAAGCGGTGATGGCATGTCCTTCTTCTGCCGCCCACGGTTGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGcc                                                                                                                                                                                         <  2:571468/251‑1 (MQ=255)                                                     gTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGGCTGGGTGCAAGCAAACTCAATCACTAACTCATCGCCGCTTACCATTTCTGCCagag                                                                                                                                                                                     <  2:546674/251‑1 (MQ=255)                                                        gcgcTGTTAAAAGCACTTTGGCAGTACGCGTATGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTTAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTCCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCCTCGCCGCTGACCATGTCTGCCAGAgctg                                                                                                                                                                                  <  1:479433/251‑1 (MQ=255)                                                        gcgcTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAgctg                                                                                                                                                                                  <  2:325680/251‑1 (MQ=255)                                                        gcgcTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAgctg                                                                                                                                                                                  >  1:422282/1‑251 (MQ=255)                                                          gcTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAgctgct                                                                                                                                                                                >  2:371940/1‑251 (MQ=255)                                                                    gCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTACTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGGGCAATCCAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCtt                                                                                                                                                                      >  2:6271/1‑250 (MQ=255)                                                                                          agaTTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGCGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGCCCTg                                                                                                                                                >  2:178410/1‑251 (MQ=255)                                                                                              tGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGGCGCGGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACTGACAGACCTGCGt                                                                                                                                             <  2:345329/250‑1 (MQ=255)                                                                                                     tttGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGGTGAGGCAGGCTTTTTGAACGGAGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGGGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGCGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTATGCCAGAGCTGCTTTTGCTTCTATCAGAGGAAACGGACAGACCGGCGCCACCACAt                                                                                                                                     >  1:60074/1‑251 (MQ=255)                                                                                                              tCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCCGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCtttt                                                                                                                            >  1:311775/1‑251 (MQ=255)                                                                                                              tCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCtttt                                                                                                                            >  1:526815/1‑251 (MQ=255)                                                                                                               ccgtcgctgtttcggcggTGGACAA‑TGTTGAGTAAGGCTTTTTTAACGGTGATGTTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTTCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGCCCATTTCTGCCAGAGCTGCTTTGTCTTCGATGAGTGGAAACGGACAGACCTGCGTCTCCACATCCAGCtttt                                                                                                                            <  1:454753/238‑1 (MQ=255)                                                                                                                     cTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCAc                                                                                                                     >  2:37859/1‑251 (MQ=255)                                                                                                                          ccggcggAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGTATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATca                                                                                                                <  2:471118/250‑2 (MQ=255)                                                                                                                                gaggaCAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGACAATCTGCTGATATTCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGCAACGGACAGACCTGCGTCACCCCATCCAGCTTTTTACTCACCATCTTTacac                                                                                                          >  2:67299/1‑251 (MQ=255)                                                                                                                                  ggaCAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACct                                                                                                        >  1:159744/1‑251 (MQ=255)                                                                                                                                     cAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACcttac                                                                                                     >  1:462073/1‑251 (MQ=255)                                                                                                                                                  aGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCgc                                                                                         >  2:23352/1‑251 (MQ=255)                                                                                                                                                           gAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTg                                                                                 <  1:292193/250‑1 (MQ=255)                                                                                                                                                             aCGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGcc                                                                               >  1:436835/1‑250 (MQ=255)                                                                                                                                                                     gCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGcc                                                                      >  1:88547/1‑251 (MQ=255)                                                                                                                                                                         cACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTAc                                                                  >  1:513584/1‑251 (MQ=255)                                                                                                                                                                            gCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCt                                                               >  1:324782/1‑251 (MQ=255)                                                                                                                                                                               gCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCg                                                            <  2:575852/251‑1 (MQ=255)                                                                                                                                                                               gCATCGCCAATCTGCTGATAATCGGGGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACGTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTg                                                              <  2:64892/249‑1 (MQ=255)                                                                                                                                                                                            gctgATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATa                                                  <  1:371940/248‑1 (MQ=255)                                                                                                                                                                                                         gATGGCATGTCCGTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCa                                  <  2:508180/251‑1 (MQ=255)                                                                                                                                                                                                            ggCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACCGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAAtt                                                                                               >  1:237446/1‑187 (MQ=255)                                                                                                                                                                                                            ggCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAAtt                                                                                               <  2:237446/187‑1 (MQ=255)                                                                                                                                                                                                                  tcattcttcctccgccCACTGTGAGATGGCTTCCGTTGTCTGAGTGCAATCAAACTCCATCACTACCTCATCGCCGCTCACCATTTCTGCCAGAGCTGCTTTTGCTTATATGAGTGGAACCGGACAGACCTGCGTCACCACCTCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCCCCTCACGCTGCGGTCGAACATACACAAGCAAGAACGACGTCCAGACt                         >  2:378094/12‑251 (MQ=255)                                                                                                                                                                                                                         tctGCCGCCCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCgag                    <  2:88547/249‑1 (MQ=255)                                                                                                                                                                                                                            gccgccCACTGTGGGATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATc                >  2:248929/1‑250 (MQ=255)                                                                                                                                                                                                                               gccCTCTGTGGGATTGCTTCAGTTGGCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTGGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATaa            <  2:450562/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggATTGCTTCAGTTGTCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATACAAGGcc  >  1:519375/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            TGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAGAGATTGTCTAATTTGAAGCCTCCGCCGCTGGTACGGCGGAGGACAA‑TGTTGAGTCAGGCTTTTTGAACGGTGATGCTCCACGCTGCATCGCCAATCTGCTGATAATCGGTGATGGCATGTCCTTCTTCTGCCGCCCACTGTGGGATTGCTTCAGTTGCCTGGGTGCAATCAAACTCAATCACTAACTCATCGCCGCTGACCATTTCTGCCAGAGCTGCTTTTGCTTCTATGAGTGGAAACGGACAGACCTGCGTCACCACATCCAGCTTTTTAATCACCATCTTTACACCTTACTTAATTTGCCGCAGCAGTTGCCAGTCGCGCCTTACGCTGCGGTCGAACATACACAAGCCAGGACGCAGTCCAGACTCCGAGAATCATAAATACAAGGCC  >  NC_000913/2084084‑2084565

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