Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	2,133,554	T→G	10.6%	L78R (CTG→CGG)	comEC →	ComE operon protein 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	2,133,554	0	T	G	10.6%	258.3 / 2.8	92	L78R (CTG→CGG)	comEC	ComE operon protein 3
Reads supporting (aligned to +/- strand):  ref base T (72/7);  new base G (2/11);  total (74/18)
Fisher's exact test for biased strand distribution p-value = 3.91e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGGCGCGCCGCT  >  ACB143/2133421‑2133683                                                                                                                                      |                                                                                                                                  gcTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTgg                                                                                                                                 >  4:162386/1‑136 (MQ=255) gcTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTg                                                                                                                                  >  5:20517/1‑135 (MQ=255) gcTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTg                                                                                                                                  >  5:334063/1‑135 (MQ=255) gcTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTg                                                                                                                                  >  1:38576/1‑135 (MQ=255) gcTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCGTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTgg                                                                                                                                 >  6:99204/1‑136 (MQ=255)   tgtgGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTGCTGCTGGGCTGGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAg                                                                                                                               <  4:387743/136‑1 (MQ=255)     tCGGCAGCCTGTTTACCCGGGGCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGGCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAg                                                                                                                               <  3:95922/134‑1 (MQ=255)        gCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcg                                                                                                                          >  6:171230/1‑136 (MQ=255)        gCAGCCTGTTTACCCGGGTCTGGCCGCGGGGTTGGTTGCGGCGGGGCTGGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTGGAGGCTCGTCTGGCGCCCGGCCGGGAGGGGCGCACCGTGGGGCGGGAGGggcg                                                                                                                          >  5:233960/1‑136 (MQ=255)         cAGCCTGTTTACCCGGGTCGGGCCGCTGGGTTGGTTGAGGCGGGGCTTGGGTGGGGCGGGCGGGTCTGCGCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcgg                                                                                                                         <  5:65546/136‑1 (MQ=255)          aGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcggg                                                                                                                        >  4:119843/1‑136 (MQ=255)               gTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcgggtggcg                                                                                                                   >  4:392312/1‑136 (MQ=255)                  tACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcgggtggcgggt                                                                                                                >  3:314595/1‑136 (MQ=255)                  tACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcgggtggcgggt                                                                                                                >  3:19318/1‑136 (MQ=255)                  tACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGggcgggtggcgggt                                                                                                                >  5:105350/1‑136 (MQ=255)                       gggcgtGGCCGCGGGGTTGGTTGCGACTGTTCTTGTGTTGGCCGGGCGGGGCTGCTCACCTGGCCGTGGATGTTCGTCTGGCTACCGGCCTGGAGGGGCGCCCCTTGTGGTTGTAGGGGCGGGTGGTGGGTTTGcc                                                                                                           <  2:266037/131‑1 (MQ=255)                           cTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCCGCAGGCCCTTGATGAGCGGCTGGCGACCGGCCTGGAGGGGCGCACCTTGTGGCGGGAGGGGCGGGTGGCGGGGTTGCCGGcc                                                                                                       >  4:312839/1‑136 (MQ=255)                           cTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGcc                                                                                                       >  1:308290/1‑136 (MQ=255)                           cTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGcc                                                                                                       >  3:7688/1‑136 (MQ=255)                           cTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACATTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGcc                                                                                                       >  7:128562/1‑136 (MQ=255)                           cTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCATGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGACGGcc                                                                                                       >  5:80980/1‑136 (MQ=255)                             ggCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCg                                                                                                     >  8:16789/1‑136 (MQ=255)                                 ggtGGGTTGGTTGCTGCGGGGCTGGTGTGGGGCGGGCGGGACTGCGCCGCCGGCCCTGGATCATCGGCGGCCTCCGGGCCGGGTGGGGCGCACCTTGTGGCGGGCGGGGCGGGGGGGGGGTTTGCCGGCCCGCACt                                                                                                 <  2:129295/134‑1 (MQ=255)                                 gCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACt                                                                                                 >  5:326904/1‑136 (MQ=255)                                      gttggttTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGcaca                                                                                            >  4:243398/1‑136 (MQ=255)                                      gttggttTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGcaca                                                                                            >  4:24599/1‑136 (MQ=255)                                      gttggttTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGcaca                                                                                            >  5:9351/1‑136 (MQ=255)                                      gttggttTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCAGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGATTGCCGGCCCGCACTGcaca                                                                                            >  8:153800/1‑136 (MQ=255)                                      gttggttTCTGCTGGGCTTGTGTTGGGCGTGCGCGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGcaca                                                                                            >  7:63046/1‑136 (MQ=255)                                             tCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGtgtg                                                                                     >  8:340076/1‑136 (MQ=255)                                                gctgGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGt                                                                                  >  3:205528/1‑136 (MQ=255)                                                gctgGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGt                                                                                  >  7:288190/1‑136 (MQ=255)                                                  tgGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGTCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGttt                                                                                >  5:90865/1‑136 (MQ=255)                                                       ttGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGc                                                                           >  8:225825/1‑136 (MQ=255)                                                            ttGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACCCGGTGTGCGTTTCGAGCTgg                                                                        >  5:294186/1‑134 (MQ=255)                                                                  gTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGTGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCGGGAAGCGCCa                                                                >  1:329669/1‑136 (MQ=255)                                                                  gTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGACGCGccc                                                                >  1:132013/1‑135 (MQ=255)                                                                  gTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCa                                                                >  4:35465/1‑136 (MQ=255)                                                                  gTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCa                                                                >  7:139193/1‑136 (MQ=255)                                                                  gTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCa                                                                >  5:220908/1‑136 (MQ=255)                                                                     cGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGc                                                             >  1:65542/1‑136 (MQ=255)                                                                      gggTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCt                                                            >  7:186946/1‑136 (MQ=255)                                                                       gggcggCTCACCAGGCCCTGGATGAGCGTCGGGCTCCCGGCCTGGCGGGGCGCACCTTGTGGCGGGCGGGGCGGGGGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTc                                                           <  3:225887/131‑1 (MQ=255)                                                                       ggTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTTGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGTTGTGCGTTGCGAGCTGGAAGCGCCACGCTc                                                           >  2:324357/1‑136 (MQ=255)                                                                       ggTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACTGTGTGCGTTTCGAGCTGGAAGCGCCACGCTc                                                           >  6:47588/1‑136 (MQ=255)                                                                       ggTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTc                                                           >  7:339290/1‑136 (MQ=255)                                                                       ggTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTc                                                           >  7:349979/1‑136 (MQ=255)                                                                       ggTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTc                                                           >  1:19162/1‑136 (MQ=255)                                                                            gCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCgcggc                                                      >  6:168976/1‑136 (MQ=255)                                                                              tcggcagGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCCCGCTCGcggcgg                                                    >  6:140984/4‑136 (MQ=255)                                                                              tCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGcggcgg                                                    >  4:154384/1‑136 (MQ=255)                                                                                 gctggCCGTGGATGAGCGGCTGGCTCCGGGCCGGGCGGGGCGCACCTGGGGGCGGTTGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGcc                                                 <  8:52686/133‑1 (MQ=255)                                                                                 ccagGCCCTTGATGACCGGCGGGCTACCGCCCTGGCGGGGCGCACCTTGTGGCGGGCGGGGGGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGcc                                                 <  8:288586/135‑1 (MQ=255)                                                                                     gCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAAc                                             >  4:30136/1‑136 (MQ=255)                                                                                     gCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAAc                                             >  6:216783/1‑136 (MQ=255)                                                                                      cccTGGATGAGCGGCTGGCTACCGGCCGGGAGGGGCGCACCTTGTGGCGGGCGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACt                                            <  6:173674/136‑1 (MQ=255)                                                                                        cTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGc                                          >  4:299947/1‑136 (MQ=255)                                                                                          tgatgaTCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCt                                        >  8:274074/1‑136 (MQ=255)                                                                                          tgatgaTCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCt                                        >  1:127753/1‑136 (MQ=255)                                                                                          tgatgaTCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGACCGAACTGCCt                                        >  6:39708/1‑136 (MQ=255)                                                                                          tgatgaTCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGCAGGGCCGAACTGCCt                                        >  8:361598/1‑136 (MQ=255)                                                                                          tgatgaTCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTGTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCt                                        >  2:173942/1‑136 (MQ=255)                                                                                           gaggagcGTGTGGCTACCGGCCTGGAGGGGCGAACCTGGTGGCTGGAGGGGGGGGGGGGGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGcc                                         <  2:304715/128‑1 (MQ=255)                                                                                            atgatCGTCTGGCACCCGGCCTGGCGGGGCGCCCCTTGGGGCGGTTGGGGCGGGGGGCGGGTTTGCCGGCCCGCACTGCACACGGGGTGCGTTTCGAGCGGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTc                                       <  3:342785/135‑1 (MQ=255)                                                                                             tgatCGTCTGGCTACCGGCCTGGAGGGTATCAAAATGTGGCTGGAGGGCCGGGTGGCGGCTATGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCCGCGGGCCGAACTGCCTCAg                                     >  8:18638/1‑136 (MQ=255)                                                                                             tgatCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGCACTGCCTCAg                                     >  2:405703/1‑136 (MQ=255)                                                                                               atCGTCTGGCTCCCGGCCGGGAGGGGCGCACCTGGGGGCGGGTGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAgc                                    <  5:369466/135‑1 (MQ=255)                                                                                                    caggCTACCGGCCTGGAGGGGCGCACCATGTGGATGGAGGGGCGGGAGGCGGCTTACCCGGCCCGCACTGAACACGCTGAGCGTTACGAGCTCGAAGCGCCACGCTCTCGGCGGGCCGAACTGCCTCAGCCCCTGc                              >  2:229168/3‑136 (MQ=255)                                                                                                    cTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCCGCGCCTg                               >  1:209427/1‑135 (MQ=255)                                                                                                    cTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTg                               >  1:230409/1‑135 (MQ=255)                                                                                                    cTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTg                               >  6:280634/1‑135 (MQ=255)                                                                                                    cTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTg                               >  6:381705/1‑135 (MQ=255)                                                                                                    cTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTg                               >  5:366475/1‑135 (MQ=255)                                                                                                    cTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTg                               >  3:318499/1‑135 (MQ=255)                                                                                                      ggCTACCGGCCTGGAGGGGCGCCCCTTTTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACATGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCGCAGCGCCTGCag                            >  4:3986/1‑136 (MQ=255)                                                                                                      ggCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCag                            >  7:226229/1‑136 (MQ=255)                                                                                                           ccGCCCTGGCGGGGCGCCCCGTGTGGGTGGCGGGGGGGGGGCCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCagctg                         <  8:301705/134‑1 (MQ=255)                                                                                                               ccTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTgg                   >  8:365833/1‑136 (MQ=255)                                                                                                                  ggAGGGGCGCACGTTGTGGCTGGCGGGGCGGGGGGCGGGTTTGCCGG‑CCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCg               <  3:26082/136‑1 (MQ=255)                                                                                                                      gggCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACg            >  8:91496/1‑136 (MQ=255)                                                                                                                      gggCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACg            >  1:385944/1‑136 (MQ=255)                                                                                                                      gggCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACg            >  5:32534/1‑136 (MQ=255)                                                                                                                      gggCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACg            >  2:161166/1‑136 (MQ=255)                                                                                                                         cgcCCCTGGGGGCGGGAGGGGGGGGGGGCGGTGTTGCCGGCCCGCACTGCACACGGGGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGgcg         <  3:324933/136‑1 (MQ=255)                                                                                                                         cgcACCTTGTGGGGGGCGGGGGGGGGGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGgcg         <  4:48272/136‑1 (MQ=255)                                                                                                                           cACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGgcgcg       >  1:177596/1‑136 (MQ=255)                                                                                                                            aCCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTTCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGgcgcgc      >  4:121196/1‑136 (MQ=255)                                                                                                                            aCCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGgcgcgc      >  3:371065/1‑136 (MQ=255)                                                                                                                            aCCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGgcgcgc      >  4:279591/1‑136 (MQ=255)                                                                                                                             ccTTGGGGCGGGCGCGGCGGGGGGCGGGGTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGGCGcgcc     <  6:159254/136‑1 (MQ=255)                                                                                                                                tgtgGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGGCGCGCCGCt  >  6:164170/1‑136 (MQ=255)                                                                                                                                 ggggCGGGAGGGGCGGGGGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGGCGCGCCGCt  <  4:158036/133‑1 (MQ=255)                                                                                                                                      |                                                                                                                                  GCTGTCGGCAGCCTGTTTACCCGGGTCTGGCCGCTGGGTTGGTTTCTGCTGGGCTTGTGTTGGGCGTGCGGGTCTGCTCAGCAGGCCCTTGATGATCGTCTGGCTACCGGCCTGGAGGGGCGCACCTTGTGGCTGGAGGGGCGGGTGGCGGGTTTGCCGGCCCGCACTGCACACGGTGTGCGTTTCGAGCTGGAAGCGCCACGCTCGCGGCGGGCCGAACTGCCTCAGCGCCTGCAGCTGAGCTGGTTCGACGGCGCGCCGCT  >  ACB143/2133421‑2133683

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