Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,823,537	A→G	L78P (CTG→CCG)	recA ←	DNA recombination/repair protein RecA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,823,537	0	A	G	95.3%	125.8 / ‑1.5	44	L78P (CTG→CCG)	recA	DNA recombination/repair protein RecA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (12/29);  total (13/31)
Fisher's exact test for biased strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

GATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTC  >  NC_000913/2823290‑2823777                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  gATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTgcggc‑g                                                                                                                                                                                                                                               <  2:281175/251‑1 (MQ=255)     tCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCTGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGcggc‑gtc                                                                                                                                                                                                                                             >  1:91128/1‑249 (MQ=255)                  ccgcgtgcggggTCAGTGCCGCCACGGAGTCAACGACGATAACGGCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCGGCGCTTCGGCGATCACCTGCGGC‑GTCAGCGTGGTTCTACCg                                                                                                                                                                                                                              <  2:344905/240‑1 (MQ=255)                                 gTGCCGCCTCGGAGGCAACGACGATACCGTCTCCTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGTCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGGTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCg                                                                                                                                                                                                               <  2:41363/251‑1 (MQ=255)                                   gccgccACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTa                                                                                                                                                                                                             <  1:79633/251‑1 (MQ=255)                                              gTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGCTTCCGGTCCGTAGAttttgccg                                                                                                                                                                                                   >  2:297579/1‑245 (MQ=255)                                                cAACGACGATAACGTCTACGGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGAt                                                                                                                                                                                                 <  2:219023/250‑1 (MQ=255)                                                     acgaTAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGcc                                                                                                                                                                                           >  1:107412/1‑251 (MQ=255)                                                          aaCGTCTACTGCGCCATAACCCGCCAGGGCGTCACAGATTTCCCGTGCCGTCTCGCCTGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGGCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGAGTCCGTTCCGTCGATTTCGACGATACGGCCCATc                                                                                                                                                                                       <  2:28149/250‑1 (MQ=255)                                                            cGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGc                                                                                                                                                                                    >  1:163739/1‑251 (MQ=255)                                                                         aGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGccc                                                                                                                                                                       <  1:266828/251‑1 (MQ=255)                                                                          gAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGCCCACCTGcccc                                                                                                                                                                      >  1:206442/1‑251 (MQ=255)                                                                             cgcgCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGccccc                                                                                                                                                                     >  2:70953/1‑248 (MQ=255)                                                                             cgcgCCAGGGCGTCACAGATTTCCAGTGCCTGCGCGCTGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAg                                                                                                                                                                   <  1:71048/251‑1 (MQ=255)                                                                              gcgcCAGGGCGTCACAGATTGCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAgc                                                                                                                                                                  <  2:62817/251‑1 (MQ=255)                                                                                                tttCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCACGCGCGATATCCAGTGaacgc                                                                                                                                                >  2:364995/1‑248 (MQ=255)                                                                                                   ccAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCTCAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGaa                                                                                                                                             <  1:56101/251‑1 (MQ=255)                                                                                                         ccTGCTCGCCGGTGTCCGGCTGGGATCACAGCAGGGTGTCGCTATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGAACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTa                                                                                                                                       <  2:258317/251‑1 (MQ=255)                                                                                                               cGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAgaga                                                                                                                                   <  1:58681/249‑1 (MQ=255)                                                                                                                      gtCCGGCTGGGAGAACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGAGCGTGTGCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATCCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTcc                                                                                                                           <  2:217137/250‑1 (MQ=255)                                                                                                                              gggAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATg                                                                                                                  >  2:30702/1‑251 (MQ=255)                                                                                                                                  gCACAGCAGGTTGGCGATATCGACGCCCAGGTTACGTGCGTAGATTGGGTCCAGGGCGTGTGCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAAc                                                                                                              <  2:342707/251‑1 (MQ=255)                                                                                                                                         agGTTGTCGAGATCGACCCCCAGTTTACGGTCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCTCAGGTTTTACCTTCCCGCTGCGCTGCGGCGATCCCCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTc                                                                                                       <  2:315961/251‑1 (MQ=255)                                                                                                                                              gTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCa                                                                                                  <  1:284734/251‑1 (MQ=255)                                                                                                                                                 gATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGCTAAACGCCCAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGgc                                                                                               <  2:71049/251‑1 (MQ=255)                                                                                                                                                 gATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGgc                                                                                               >  1:44266/1‑251 (MQ=255)                                                                                                                                                 gATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGgc                                                                                               <  1:237011/251‑1 (MQ=255)                                                                                                                                                                  tACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCttt                                                                               <  1:61241/250‑1 (MQ=255)                                                                                                                                                                  tACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCttt                                                                               >  1:54176/1‑250 (MQ=255)                                                                                                                                                                   aCGTGCGTAGATGTGGTCCAGCTCGTGTTCCGCATCTATAACCGCACAGGTTTTCCCTTCCCGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCAGCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTAc                                                                             <  2:323098/251‑1 (MQ=255)                                                                                                                                                                     gTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTAcc                                                                            >  1:228902/1‑250 (MQ=255)                                                                                                                                                                           tgATTGTGTCCAGCGCGTCTTCAGCATCGAGATACGCAGTGGTTTTCCCTTCCCGCGGCGCTGCGGCGTTCACCTGC‑GCAGTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGAGACGGCCCATCTGCAGACCTCCTGCACCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGGTTCCACACCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTg                                                                      <  2:321253/249‑1 (MQ=255)                                                                                                                                                                            gATTGGGTCCTGCGCGTGTTCAGCATCGATAAACGCCCAGGGTTTACCTTCACGCTGCGCGGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGGCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGtt                                                                    <  2:302868/251‑1 (MQ=255)                                                                                                                                                                               tGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTtc                                                                  <  1:27438/250‑1 (MQ=255)                                                                                                                                                                                       gcgcggtgtCAGCATCGAAAAACGCACAGGTTTTCCCTTCCCGCTGCCGTGCGGCCACTATCGGCGGC‑GTTAGCGTGGTTTTCCCGGAAGATTCCGGTCCGTAGATTTCGACTATACGGCCCAACGGCAGACCTCCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCTGTAGGGATGGTTTCCACATCCCTGGAACGGTCTTCACGCAGGCGCATGATGGAGCCTTTACCAAATTCGTTCGCAATCTgg                                                         <  2:124521/243‑1 (MQ=255)                                                                                                                                                                                               tCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGtg                                                  <  1:121626/250‑1 (MQ=255)                                                                                                                                                                                                      cGATAAACGCCCAGGTTTTCCCTTCCCGCGGCGCTGCGGCGATCACCTGCAGC‑GTCAGCGGTGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATCGTTGCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTgccgc                                           <  2:85008/250‑1 (MQ=255)                                                                                                                                                                                                      cGATAAACGCACAGGTTTTACCTTAACGCTGCGATGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAACCGGGATATACAGTGAAGGCGGACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCATCACCCCGGAGAATGATGGAGCCTTTCCCAAAATTTTTCTCAATTTGGCACAGTGCTgccggc                                          >  2:62790/1‑249 (MQ=255)                                                                                                                                                                                                          aaaCGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGTGCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCATACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAAc                                       <  2:240154/250‑1 (MQ=255)                                                                                                                                                                                                          aaaCGCACAGGTTTTACCGTCACGCGGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATGCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACg                                      <  1:20172/251‑1 (MQ=255)                                                                                                                                                                                                                   ggTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCtgtt                             <  1:372569/251‑1 (MQ=255)                                                                                                                                                                                                                          ccTTCCCGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTtc                      <  2:44266/251‑1 (MQ=255)                                                                                                                                                                                                                            ttCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCTGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCt                                                           <  2:327663/212‑1 (MQ=255)                                                                                                                                                                                                                            ttCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCTGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCt                                                           >  1:327663/1‑212 (MQ=255)                                                                                                                                                                                                                             tCACGCTGCGCTGCGGCGATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTtcg                     <  1:57223/249‑1 (MQ=255)                                                                                                                                                                                                                                              gATCACCTGCGGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTc  >  1:329490/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  GATTTCGCCTTCGATTTCCGCTTTCGGCGTCAGTGCCGCCACGGAGTCAACGACGATAACGTCTACTGCGCCAGAACGCGCCAGGGCGTCACAGATTTCCAGTGCCTGCTCGCCGGTGTCCGGCTGGGAGCACAGCAGGTTGTCGATATCGACGCCCAGTTTACGTGCGTAGATTGGGTCCAGCGCGTGTTCAGCATCGATAAACGCACAGGTTTTACCTTCACGCTGCGCTGCGGCGATCACCTGCAGC‑GTCAGCGTGGTTTTACCGGAAGATTCCGGTCCGTAGATTTCGACGATACGGCCCATCGGCAGACCACCTGCCCCAAGCGCGATATCCAGTGAAAGCGAACCGGTAGAGATGGTTTCCACATCCATGGAACGGTCTTCACCCAGGCGCATGATGGAGCCTTTACCAAATTGTTTCTCAATCTGGCCCAGTGCTGCCGCCAACGCTTTCTGTTTGTTTTCGTCGATAGCCATTTTTACTC  >  NC_000913/2823290‑2823777

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