BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,098,522	G→C	5.4%	intergenic (+119/+36)	yggW → / ← yggM	HemN family putative oxidoreductase/DUF1202 family putative secreted protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,098,522	0	G	C	5.4%	525.3 / 15.7	149	intergenic (+119/+36)	yggW/yggM	HemN family putative oxidoreductase/DUF1202 family putative secreted protein
Reads supporting (aligned to +/- strand): ref base G (71/70); new base C (3/5); total (74/75)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 155 total aligned reads displayed.

GTATACCGGGCTTTGCGAAGATGTGATTCGCCCACAGTTAGACGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCATCGTTAACGGATATTGCCCTTTTTCGAGCGTTACCCAGTCATCCTGCAATTTACGGCTAAGTTTTTG  >  NC_000913/3098253‑3098749

gTATACCGGGCTTTGCGAAGATGTGATTCGCCCACAGTTAGACGAGGCGATTGCCCCGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGACCGCCTTATCCAGCCGACATGTTGCAGCGGTTGTGGGTCTGAACAGACGCGCAAGCGTCGCATCCGAGGTTGATTGCCCGATGCGGGGGCCAccccgc                                                                                                                                                                                                                                 >  4:113067/1‑274 (MQ=255)

   taCCGGGCTTTGCGAAGATGTGATTGGACCACAGTTAGACGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATc                                                                                                                                                                                                                              <  2:158307/274‑1 (MQ=255)

                   gATGTGATTCGCCCACAGTTAGACGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCa                                                                                                                                                                                                             <  2:113067/275‑1 (MQ=255)

                              cccACAGTTAGACGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCaaa                                                                                                                                                                                                   <  2:160066/274‑1 (MQ=255)

                                cacaGTTAGACGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAact                                                                                                                                                                                                 <  2:284923/274‑1 (MQ=255)

                                         aCGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCTGATGCGGGGTCAACGCCGCATCCTGCTACCAATCGTGCACTATATCAAACTTACCTAATc                                                                                                                                                                                        >  2:263423/1‑274 (MQ=255)

                                             ggCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACTTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAACAcccc                                                                                                                                                                                   >  2:48699/1‑275 (MQ=255)

                                                               tctcACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTAGATCAAACTTACTTAATCAACCCATTGCGCCGCAGTAGTTc                                                                                                                                                                 >  2:185762/1‑275 (MQ=255)

                                                                    ccGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCCGACGTTGATTGCCGGATGCGTGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTAt                                                                                                                                                            >  2:12990/1‑275 (MQ=255)

                                                                                      ggCAGATAACGGAACATGGGAATCTGTTTTTAAATTCGTTGCTGGGGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCAt                                                                                                                                          <  4:263423/275‑1 (MQ=255)

                                                                                                                           gctgctGGATCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                                                                                                                   <  6:1071788/149‑1 (MQ=255)

                                                                                                                                tgGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                                                                                                              <  5:618024/149‑1 (MQ=255)

                                                                                                                                tgGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                                                                                                              <  7:210043/149‑1 (MQ=255)

                                                                                                                                 ggAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATcc                                                                                                                                                                                                                             >  7:470389/1‑149 (MQ=255)

                                                                                                                                    gCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGc                                                                                                                                                                                                                          <  5:1431836/149‑1 (MQ=255)

                                                                                                                                    gCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGc                                                                                                                                                                                                                          <  8:186601/149‑1 (MQ=255)

                                                                                                                                     cTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCt                                                                                                                                                                                                                         >  6:63403/1‑149 (MQ=255)

                                                                                                                                      tttttCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGCGTTGACCCCGCATc                                                                                                                                                                                                                              <  5:383073/143‑1 (MQ=255)

                                                                                                                                            tGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATc                                                                                                                                                                                                                  <  6:240591/149‑1 (MQ=255)

                                                                                                                                            tGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATc                                                                                                                                                                                                                  <  6:483830/149‑1 (MQ=255)

                                                                                                                                               cTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTg                                                                                                                                                                                                               <  8:829005/149‑1 (MQ=255)

                                                                                                                                                tGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGc                                                                                                                                                                                                              >  7:149841/1‑149 (MQ=255)

                                                                                                                                                  aGTAAACTTGTATTGCCTGATGCGGCCTGCACGCCTTATCCAGCCGACATGTGGCAGCGGCTGTAGGTCTGATAAGACGCGCAAGCGTCGCAGCAGACGTTGCTTGCCGGATCCGGGTTCAACGCCGCATCCTGCTACAAACCGTTCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCCGCGTTAATTAGCTTGGTTGATATCGTCAATTTGTATCAAAACCACGACTTTTTTACTGTTCTTTTCtt                                                                              >  4:328638/1‑274 (MQ=255)

                                                                                                                                                  aGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTatat                                                                                                                                                                                                       >  3:21810/1‑154 (MQ=255)

                                                                                                                                                      aaCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCaaa                                                                                                                                                                                                   <  1:21810/154‑1 (MQ=255)

                                                                                                                                                         ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCa                                                                                                                                                                                                     >  6:516447/1‑149 (MQ=255)

                                                                                                                                                           gTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCaaa                                                                                                                                                                                                   <  6:1159434/149‑1 (MQ=255)

                                                                                                                                                            tATTGCCGGATGCGGCGTGAACGCCTTCTCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGCTAAGACGCGCCAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCAATTGCGCCGCAGTTGTTCTGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTAATATCGTCAATTTTCTTCAAATTCGCTACTTTTTTCCTTTTTTTTTCTGAACCtttaac                                                                    >  4:201488/1‑272 (MQ=255)

                                                                                                                                                               tGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAactta                                                                                                                                                                                               >  6:598807/1‑149 (MQ=255)

                                                                                                                                                                    gATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTaa                                                                                                                                                                                          <  5:1329794/149‑1 (MQ=255)

                                                                                                                                                                     aTGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAt                                                                                                                                                                                         <  5:286677/149‑1 (MQ=255)

                                                                                                                                                                     aTGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAt                                                                                                                                                                                         <  7:654221/149‑1 (MQ=255)

                                                                                                                                                                      tGCGGCGTGAACGCCTTACCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATAAGACGTTGATTTCCGGAAGCGGGGTCAACGCCGCAGCCTGCTACAAACCGTGCACTATAGCATANTTACTTAATCACCCCATTTCGCCCCATTTGTTCGGTATCTTTTTTCAAGCGTCCATTTCCCTGGTTTACATCGTAACTTCTCAGCAAAATCACCACTTTTGTAcacntnn                                                                                    >  3:64297/1‑243 (MQ=255)

                                                                                                                                                                             tGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAt                                                                                                                                                                                 >  5:731607/1‑149 (MQ=255)

                                                                                                                                                                              gAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                                                                                                                >  8:1427415/1‑149 (MQ=255)

                                                                                                                                                                                  gCCTTTTCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATAAGACGTTGATTGCCGGATGCGGGGTCAACGCGGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATAACCCCATTgcgc                                                                                                                                                                            >  5:425440/1‑149 (MQ=255)

                                                                                                                                                                                    cTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGcgccg                                                                                                                                                                          >  7:314948/1‑149 (MQ=255)

                                                                                                                                                                                      tATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCa                                                                                                                                                                        <  7:672394/149‑1 (MQ=255)

                                                                                                                                                                                        tCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgt                                                                                                                                                                      >  6:841733/1‑149 (MQ=255)

                                                                                                                                                                                            gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                                                                                                                >  5:580252/1‑135 (MQ=255)

                                                                                                                                                                                            gCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCAtt                                                                                                                                                                                <  6:580252/135‑1 (MQ=255)

                                                                                                                                                                                             ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                                                                                                                 <  7:322593/149‑1 (MQ=255)

                                                                                                                                                                                             ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                                                                                                                 <  5:63403/149‑1 (MQ=255)

                                                                                                                                                                                             ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                                                                                                                 <  8:1053043/149‑1 (MQ=255)

                                                                                                                                                                                              cGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                                                                                                                <  5:1170119/149‑1 (MQ=255)

                                                                                                                                                                                                    gtgGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCt                                                                                                                                                          <  8:679334/149‑1 (MQ=255)

                                                                                                                                                                                                      ggCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCttt                                                                                                                                                        >  6:76008/1‑149 (MQ=255)

                                                                                                                                                                                                       gCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                                                                                                                   >  5:766969/1‑73 (MQ=255)

                                                                                                                                                                                                       gCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                                                                                                                   <  6:766969/73‑1 (MQ=255)

                                                                                                                                                                                                         aGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGACGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCtttttt                                                                                                                                                     <  6:1181840/149‑1 (MQ=255)

                                                                                                                                                                                                          gCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCa                                                                                                                                                                                       >  7:1371099/1‑114 (MQ=255)

                                                                                                                                                                                                          gCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCa                                                                                                                                                                                       <  8:1371099/114‑1 (MQ=255)

                                                                                                                                                                                                               tGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTCTATTTCTGAAACTTTCATCGTTAACGGATATTGCCCTTTTTCGAGCGTTACCCAGTCATCCTGCAAttt                 <  2:84912/275‑1 (MQ=255)

                                                                                                                                                                                                                gTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGt                                                                                                                                              >  6:359643/1‑149 (MQ=255)

                                                                                                                                                                                                                          tAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCtt                                                                                                                                    <  8:340231/149‑1 (MQ=255)

                                                                                                                                                                                                                            aGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGc                                                                                                                                      <  1:52440/145‑1 (MQ=255)

                                                                                                                                                                                                                              aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                                                                                                                <  5:1430397/149‑1 (MQ=255)

                                                                                                                                                                                                                              aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                                                                                                                <  5:767519/149‑1 (MQ=255)

                                                                                                                                                                                                                              aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                                                                                                                <  7:608876/149‑1 (MQ=255)

                                                                                                                                                                                                                              aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                                                                                                                 <  6:376457/148‑1 (MQ=255)

                                                                                                                                                                                                                              aCGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCATCGTTAACGGATATTGCCCTTTTTCGAGCGTTACCCAGTCATCCTGCAATTTACGGCTAAGTTTTtg  >  2:229809/1‑275 (MQ=255)

                                                                                                                                                                                                                               cgcgCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTg                                                                                                                               <  8:699686/149‑1 (MQ=255)

                                                                                                                                                                                                                                gcgcAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGa                                                                                                                              <  8:564841/149‑1 (MQ=255)

                                                                                                                                                                                                                                    aaGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                                                                                          <  7:1093069/149‑1 (MQ=255)

                                                                                                                                                                                                                                      gCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGt                                                                                                                        >  8:59868/1‑149 (MQ=255)

                                                                                                                                                                                                                                      gCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGt                                                                                                                        <  5:1091279/149‑1 (MQ=255)

                                                                                                                                                                                                                                         tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                                                                                                                 >  7:935630/1‑137 (MQ=255)

                                                                                                                                                                                                                                         tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                                                                                                                 <  8:935630/137‑1 (MQ=255)

                                                                                                                                                                                                                                         tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCaa                                                                                                                     <  6:884240/149‑1 (MQ=255)

                                                                                                                                                                                                                                             aTCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAAtttt                                                                                                                 >  6:915898/1‑149 (MQ=255)

                                                                                                                                                                                                                                              tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtc                                                                                                                >  7:854585/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttactt                                                                                                                                                                                            <  5:32571/70‑1 (MQ=255)

                                                                                                                                                                                                                                                 gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttactt                                                                                                                                                                                            >  6:32571/1‑70 (MQ=255)

                                                                                                                                                                                                                                                 gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatc                                                                                                             >  6:1120600/1‑149 (MQ=255)

                                                                                                                                                                                                                                                 gACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatc                                                                                                             <  6:84808/149‑1 (MQ=255)

                                                                                                                                                                                                                                                  aCGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGatat                                                                                                                           <  8:58163/134‑1 (MQ=255)

                                                                                                                                                                                                                                                  aCGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcatca                                                                                                            >  8:52792/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaaa                                                                                                          >  7:35060/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCCTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCaaa                                                                                                          >  8:729073/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       gATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATc                                                                                                       >  7:835918/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        aTTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATcg                                                                                                      >  8:1209291/1‑149 (MQ=255)

                                                                                                                                                                                                                                                         ttGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATcgc                                                                                                     <  8:470389/149‑1 (MQ=255)

                                                                                                                                                                                                                                                            ccGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGAc                                                                                                  >  7:623707/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                                                                                                                 <  8:1317063/117‑1 (MQ=255)

                                                                                                                                                                                                                                                             cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGt                                                                                                                                 >  7:1317063/1‑117 (MQ=255)

                                                                                                                                                                                                                                                             cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                                                                                                 >  5:204205/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             cGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACt                                                                                                 >  5:311568/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               gATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttt                                                                                               >  5:708425/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               gATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttt                                                                                               <  7:404708/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               gATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttt                                                                                               >  7:1188297/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               gATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATc                                                                                                                          >  8:1115117/1‑122 (MQ=255)

                                                                                                                                                                                                                                                                  gCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATcgc                                                                                                     <  4:169504/140‑1 (MQ=255)

                                                                                                                                                                                                                                                                    ggggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTAc                                                                                          <  7:52792/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     gggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACt                                                                                         >  8:1448588/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                        tCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt                                                                                      >  8:354401/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         cAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                                                                                                                 >  6:937657/1‑73 (MQ=255)

                                                                                                                                                                                                                                                                         cAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                                                                                                                 <  5:937657/73‑1 (MQ=255)

                                                                                                                                                                                                                                                                          aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa                                                                                    >  7:861705/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                          aaCGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTa                                                                                    >  7:457669/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                             gccgcATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTAttt                                                                                 <  7:876605/149‑1 (MQ=255)

GTATACCGGGCTTTGCGAAGATGTGATTCGCCCACAGTTAGACGAGGCGATTGCCCAGGGTTATCTCACCGAATGTGCGGATTACTGGCAGATAACGGAACATGGGAAGCTGTTTTTAAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTTTATTTCTGAAACTTTCATCGTTAACGGATATTGCCCTTTTTCGAGCGTTACCCAGTCATCCTGCAATTTACGGCTAAGTTTTTG  >  NC_000913/3098253‑3098749

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 26 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`