BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,662,505	G→T	5.9%	intergenic (‑374/+78)	iscR ← / ← trmJ	isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,662,505	0	G	T	5.9%	486.2 / 11.4	150	intergenic (‑374/+78)	iscR/trmJ	isc operon transcriptional repressor; suf operon transcriptional activator; oxidative stress‑ and iron starvation‑inducible; autorepressor/tRNA mC32,mU32 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base G (84/55); new base T (7/4); total (91/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

AGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCT  >  NC_000913/2662358‑2662653

cggCCTGATACCACGCATTACGCGGCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGaa                                                                                                                                                    <  2:462542/149‑1 (MQ=255)

  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGTTGCGGCGTGACCGCGTGCTCCGACCTACGGTACAGTTCACTGTTGGACGGATTAGTCGCATGACGGGTCGCATCAGGCAACGGCTGTCGGAGGGGGCGTGAACg                                                                                                                                                  >  2:2393576/1‑150 (MQ=11)

  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                                                                                  >  2:1063214/1‑150 (MQ=255)

  gCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACg                                                                                                                                                  >  1:1691264/1‑150 (MQ=255)

   ccTGATTAGACGCTTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGCCCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGc                                                                                                                                                 <  2:578705/150‑1 (MQ=255)

   ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGc                                                                                                                                                 <  2:197849/150‑1 (MQ=255)

           gACGCATGACGCGTCGCACCAGGCAACGGCTGTCGGAGGCGGCGGGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                         <  2:78258/150‑1 (MQ=255)

            aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                                                        >  1:1162072/1‑150 (MQ=255)

            aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                                                        >  1:217347/1‑150 (MQ=255)

            aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGa                                                                                                                                        >  1:2090010/1‑150 (MQ=255)

                aTTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTa                                                                                                                                    >  2:2041519/1‑150 (MQ=255)

                       ggcGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCt                                                                                                                             <  2:2326869/148‑1 (MQ=255)

                       gTCGCATCTGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCt                                                                                                                             <  2:2041673/150‑1 (MQ=255)

                       gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCt                                                                                                                             >  1:2372080/1‑150 (MQ=255)

                       gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCt                                                                                                                             >  1:2372074/1‑150 (MQ=255)

                          gCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTTTTCACTGTAGGCCTGATAAGACGCATTACGCGACGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGtt                                                                                                                          <  2:634065/150‑1 (MQ=255)

                           caccAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                                                         >  1:1020779/4‑150 (MQ=255)

                           cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                                                         >  1:1411959/1‑150 (MQ=255)

                           cATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGGGAACGCCTTATCCGACCTACGGTTCTGTTc                                                                                                                         >  2:1117896/1‑150 (MQ=255)

                           cATCAGGCAACGGCTGTCGGATGCCGCGTGAACGCCTAATCCGACCAACGGTTCTGTTCACTGTAGGCCTGATAAGAAGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGAACTACGGTTCTGTTc                                                                                                                         >  1:1824256/1‑150 (MQ=17)

                            aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCa                                                                                                                        >  2:41533/1‑150 (MQ=255)

                            aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCACATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCa                                                                                                                        >  1:258699/1‑150 (MQ=255)

                            aTCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGGTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCACCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCa                                                                                                                        >  1:1581467/1‑150 (MQ=255)

                              cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCATGAACGCCTTATCCGACCTACGGTTCTGTTCACt                                                                                                                      <  2:189258/150‑1 (MQ=255)

                                 gaaaCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                                                                   <  2:1171625/148‑1 (MQ=255)

                                 ccaccGGGTGTGGGTTGCGGCGTTAACGCCTTATCCGACCTTGGTTTCTGTTCACTGTAGCCCTTATAAGCCGATTTACGCGCCGAATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTa                                                                                                                   <  1:2231334/146‑1 (MQ=11)

                                  cAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAg                                                                                                                  <  2:243099/150‑1 (MQ=255)

                                           gcGGAGGCGGCGGGACCGCCTTCTCCGCCCTCCGGTTCGGTTCACTGTAGGCCTGATAAGACGCATTACGCGGCGCATCAGGCAACGGCGGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATaa                                                                                                         <  2:1674371/149‑1 (MQ=11)

                                              gataaggcgTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATAAGGCAACGGCTGTAGGGTGCGGCGTGAACCCCTAATCCGACCAACGGTTGAGTTCACAATAGGCCGGATAATAc                                                                                                      >  2:315724/6‑150 (MQ=14)

                                              gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGAc                                                                                                      >  1:1206496/1‑150 (MQ=255)

                                                    gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGCTTCTGTTCACAGTAGGCCTCATAAGACTCATTa                                                                                                >  2:1028777/1‑150 (MQ=255)

                                                    gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCAGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                                                                                >  2:470463/1‑150 (MQ=255)

                                                    gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTAATCCGACCAACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTa                                                                                                >  2:1908396/1‑150 (MQ=255)

                                                          aCGCCTTATCCGACCTTCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                          <  2:805808/150‑1 (MQ=255)

                                                          aCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                          >  2:1479040/1‑150 (MQ=255)

                                                          aCGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTc                                                                                          >  2:1139580/1‑150 (MQ=255)

                                                              cTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCAt                                                                                      >  1:1680602/1‑150 (MQ=40)

                                                               ttAACCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTAGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATc                                                                                     <  2:1928832/150‑1 (MQ=255)

                                                                tATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATGCGACGCGTAATGCGTCTTATCa                                                                                    <  2:928358/150‑1 (MQ=14)

                                                                    cGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCa                                                                                >  2:2299532/1‑150 (MQ=38)

                                                                     gcccTGCGGTTCTGGTCACTGTAGGCCTGAGTCGACGCTTTAGGGGTGGGAGCTGGCAACGGCTGTCGGAGGCGGCGTGAACGCCTTATCCGACCTTCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                               <  1:567245/148‑1 (MQ=11)

                                                                     gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                               >  2:820601/1‑150 (MQ=38)

                                                                     gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                               >  1:414756/1‑150 (MQ=38)

                                                                     gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                               >  2:2284616/1‑150 (MQ=38)

                                                                     gACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCaa                                                                               >  2:399307/1‑150 (MQ=38)

                                                                         tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                           >  1:1521531/1‑150 (MQ=34)

                                                                         tACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGc                                                                           >  1:1521422/1‑150 (MQ=34)

                                                                            gggtCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTAATCCGACCTACGGTTCTGTACACTGTAGGCCTGATAATACGCATTACGCGTCGCATCAGGCAACGGCTGt                                                                        >  2:175536/4‑150 (MQ=255)

                                                                              ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCACCGGCGGTTa                                                                      >  2:1765108/1‑150 (MQ=18)

                                                                              ttCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCGACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTa                                                                      >  2:2253478/1‑150 (MQ=18)

                                                                                 tGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACACATTACGCGTCGCATCAGGCAACGGCTGTTATTc                                                                   >  2:1929341/1‑150 (MQ=25)

                                                                                   ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCTTGAACGCCTTATCCGACCAACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATACgg                                                                 >  2:1494131/1‑150 (MQ=255)

                                                                                   ttCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCgg                                                                 >  2:1494522/1‑150 (MQ=37)

                                                                                    tCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGc                                                                >  1:1602824/1‑150 (MQ=37)

                                                                                     cACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGcc                                                               >  1:499495/1‑150 (MQ=37)

                                                                                              ccTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCt                                                      <  2:1602824/150‑1 (MQ=255)

                                                                                               cTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCtt                                                     >  2:379408/1‑150 (MQ=255)

                                                                                                    aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTc                                                >  1:1238498/1‑150 (MQ=255)

                                                                                                    aaGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTc                                                >  1:258384/1‑150 (MQ=255)

                                                                                                     aGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTct                                               >  2:813170/1‑150 (MQ=255)

                                                                                                     aGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTct                                               >  1:2341568/1‑150 (MQ=255)

                                                                                                     aGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTct                                               >  2:1022281/1‑150 (MQ=255)

                                                                                                       aCGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgc                                             >  1:84484/1‑150 (MQ=255)

                                                                                                         gCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                           >  2:92616/1‑150 (MQ=255)

                                                                                                         gCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTctgctg                                           >  2:1163982/1‑150 (MQ=255)

                                                                                                          cATTCCGCGGCGCAGCAGGCCACGGGTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                                          <  2:1235244/150‑1 (MQ=17)

                                                                                                          cATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                                          <  2:2372074/150‑1 (MQ=255)

                                                                                                          cATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCtgctgc                                          <  2:2372073/150‑1 (MQ=255)

                                                                                                            ttACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTc                                        >  1:182064/1‑150 (MQ=255)

                                                                                                                         cAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGa                           >  2:1708918/1‑150 (MQ=255)

                                                                                                                          aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                          >  2:2199844/1‑150 (MQ=255)

                                                                                                                          aGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGaa                          >  2:338474/1‑150 (MQ=255)

                                                                                                                            gCATCGGCTGTCGGATGCGGCGTGAACGCCTTTTCCGACCTACGGTTCTGTTCACTGTTGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATc                        <  2:1284111/150‑1 (MQ=255)

                                                                                                                            gCAACGGCTGTCGGATGCGGCGTTAACGCATTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATc                        <  2:222542/150‑1 (MQ=255)

                                                                                                                               ccGGGTGTCGGTTGGGGGGTGAACGCCTTAGCCGACCTCGGGTTCTGTTCAGTGTTGGCCTGATAAGACGCATTACGGGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                     <  2:1045887/149‑1 (MQ=255)

                                                                                                                               aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                     <  2:2102839/150‑1 (MQ=255)

                                                                                                                               aCGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCcg                     <  2:1222746/150‑1 (MQ=255)

                                                                                                                                     gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAgg               <  2:720707/150‑1 (MQ=255)

                                                                                                                                     gTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAgg               <  2:199585/150‑1 (MQ=255)

                                                                                                                                      tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGTCGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTAGCTTTATTCTGCTGCTCAATAGAAGACAGAATCCCGCGCAGGa              >  2:1664273/1‑150 (MQ=21)

                                                                                                                                      tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa              >  1:403649/1‑150 (MQ=25)

                                                                                                                                      tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa              >  1:621140/1‑150 (MQ=25)

                                                                                                                                      tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa              >  1:1018088/1‑150 (MQ=25)

                                                                                                                                      tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa              >  1:403663/1‑150 (MQ=25)

                                                                                                                                      tCGGATAAGGCGTTCACGCCGCATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGa              >  2:2386756/1‑150 (MQ=25)

                                                                                                                                        ggATGCGGCGCGAACGCATTATCCGACCTCCGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCCGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGata            <  1:1989226/150‑1 (MQ=255)

                                                                                                                                         gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat           >  1:1451942/1‑150 (MQ=255)

                                                                                                                                         gATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGatat           >  1:1451661/1‑150 (MQ=255)

                                                                                                                                           tGCGGCTTTAACGCCTTCTCAGACCTAGGTTTTTTTTCACTTTCGGCCTTATAATTCGAATTACGCTTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc         <  1:1862771/150‑1 (MQ=25)

                                                                                                                                           tGCGGCGTTAACGCATTATCCGCCCTTGGTTTGTTTTCAGTGTTGGCTTTATAAGACGAATTACGCGTCGAATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCTGCGCGGGattctg         <  1:768308/150‑5 (MQ=255)

                                                                                                                                           tGCGGCGTGAACGCCTTATCCGACCTTGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc         <  2:698832/150‑1 (MQ=255)

                                                                                                                                           tGCGGCGTGAACGCCTTATCCGACCTCGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc         <  2:2074180/150‑1 (MQ=255)

                                                                                                                                           tGCGGCGTGAACGCCTTATCCGACCTAGGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTc         <  2:1497595/150‑1 (MQ=255)

                                                                                                                                           tGCGGCGTGAACGCCTTACCCGCCCTGGGGTTCTGTACCCTTTAGGCATGATAAGACGAATGCCGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCTGCGCGGGattctg         <  1:2198810/150‑5 (MQ=18)

                                                                                                                                           ggcggcgTTACCGCCTTCTCCGACCTCCGGTTGTGTACAGTGTCGGCGTTATATGACGCATTCCGCGTCGCAGCCGGCCACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCTGCGCGGGattctg         <  1:739879/149‑5 (MQ=255)

                                                                                                                                             cggcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCATCGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCaa       >  2:1300853/1‑150 (MQ=255)

                                                                                                                                               gcgTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACt     >  1:991852/1‑150 (MQ=255)

                                                                                                                                                  tGACCGCCTTATCCGACCTACGGTTCTGTTCAGTGTTGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt  <  1:349641/150‑1 (MQ=255)

                                                                                                                                                  tGAACGCCTTATCCGACCTTGGGTTCTGTACACTGTAGGCCTTATACGACGCATTACGCGTCGCATCTGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt  <  1:349653/150‑1 (MQ=37)

                                                                                                                                                  tGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCt  >  1:2315848/1‑150 (MQ=255)

AGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTGAACGCCTTATCCGACCTACGGTTCTGTTCACTGTAGGCCTGATAAGACGCATTACGCGTCGCATCAGGCAACGGCTGTTATTCGGCCTTGTTACCTTTATTCTGCTGCTCAATAGAAGCCAGAATCCCGCGCAGGATATTCAACTCCT  >  NC_000913/2662358‑2662653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`