BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,412,779	T→G	27.2%	V131V (GTT→GTG)	yjfM →	DUF1190 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,412,779	0	T	G	27.2%	363.6 / 74.4	205	V131V (GTT→GTG)	yjfM	DUF1190 family protein
Reads supporting (aligned to +/- strand): ref base T (36/110); new base G (57/2); total (93/112)
Fisher's exact test for biased strand distribution p-value = 2.89e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 208 total aligned reads displayed.

AATAATGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGA  >  NC_000913/4412630‑4412925

   aatGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCt                                                                                                                                                   >  1:2736025/1‑150 (MQ=255)

    atGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTg                                                                                                                                                  >  2:3962605/1‑150 (MQ=255)

      gCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGa                                                                                                                                                >  2:623710/1‑150 (MQ=255)

       ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑t                                                                                                                                              >  2:1985880/1‑150 (MQ=255)

       ccGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑t                                                                                                                                              >  2:1426246/1‑150 (MQ=255)

          aCATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑tttt                                                                                                                                           <  2:4790348/150‑1 (MQ=255)

           cATTTGCGCGCGTGGCTTGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑ttttt                                                                                                                                          <  1:4427065/150‑1 (MQ=255)

           cATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑ttttt                                                                                                                                          <  1:762212/150‑1 (MQ=255)

           cATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑ttttt                                                                                                                                          <  1:4642166/150‑1 (MQ=255)

           cATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑ttttt                                                                                                                                          <  1:3173908/150‑1 (MQ=255)

             tttGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGt                                                                                                                                        <  1:4525529/150‑1 (MQ=255)

             tttGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGt                                                                                                                                        <  1:3089021/150‑1 (MQ=255)

               tGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGG‑ttttttttt                                                                                                                                      >  1:1292311/1‑149 (MQ=255)

                   cgcgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGG‑TTTTTTTTAAGcc                                                                                                                                  >  1:2827957/1‑150 (MQ=255)

                     cgTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCgg                                                                                                                                <  1:4034205/150‑1 (MQ=255)

                        ggCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGtt                                                                                                                             <  1:3265719/150‑1 (MQ=255)

                           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTTAGCCGGGTTAtt                                                                                                                          >  1:2006497/1‑150 (MQ=255)

                           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTAtt                                                                                                                          <  1:2742132/150‑1 (MQ=255)

                           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGGTAtt                                                                                                                          >  2:4019474/1‑150 (MQ=255)

                           tGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTGTTTAGCCGGGTTAtt                                                                                                                          >  1:4312340/1‑150 (MQ=255)

                            ggAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTc                                                                                                                         <  2:1684598/150‑1 (MQ=255)

                            ggAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTc                                                                                                                         <  2:1346544/150‑1 (MQ=255)

                            ggAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTc                                                                                                                         <  2:1088979/150‑1 (MQ=255)

                            ggAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGA‑TTTTTTTTTATCCGGGTTAtta                                                                                                                         >  2:2179902/1‑149 (MQ=255)

                              aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACCGAGTTGGATCCCGGTCGTGTGGGGA‑TTTTTTTTTATCCGGGGTAttggg                                                                                                                       >  1:3708543/1‑147 (MQ=255)

                              aacaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGG‑TTTTTTTTAAGACGGGTTTTTA‑Gt                                                                                                                      >  1:2453611/1‑150 (MQ=255)

                               acaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGT‑TTTTTGTTAATCCGGGGTATTC‑GTa                                                                                                                     >  2:551097/1‑150 (MQ=255)

                               acaacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTATGGA‑TTTTTTTTAAGCCGGGGTATTC‑gga                                                                                                                     >  2:924992/1‑148 (MQ=255)

                                 aacGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGGGGA‑TTTTTGTTAAGCCGGGTTATTA‑GTaaa                                                                                                                   >  2:3906495/1‑150 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGA‑TTTTTGTTTAGCCGGGTTATTA‑GTAAAGAt                                                                                                                >  2:3336773/1‑150 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGAt                                                                                                                <  2:4219525/150‑1 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATTGA‑TTTTTTTTAAACCGGGTTTTTA‑GTAAAGAt                                                                                                                >  2:2593774/1‑150 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGC‑TTTTTGTTTATCCGGGGGATTG‑GTGAAGAt                                                                                                                >  1:1652206/1‑150 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGGGT‑TTTTTGTTTATCCGGGGTATTG‑GTAAAGAt                                                                                                                >  2:87832/1‑150 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGG‑TTTTTGTTAAGCCGGGGTAGTA‑GTAAAGaa                                                                                                                >  2:4642021/1‑149 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTGTTAATACGGGGTATTC‑GTAAAGAt                                                                                                                >  1:3516287/1‑150 (MQ=255)

                                    gCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAatat                                                                                                                >  2:4519127/1‑150 (MQ=255)

                                     ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATc                                                                                                               <  2:3406887/150‑1 (MQ=255)

                                     ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATc                                                                                                               <  2:143208/150‑1 (MQ=255)

                                     ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTTTTAAGCCGGGGTATTC‑GTAAAGATc                                                                                                               >  2:297785/1‑150 (MQ=255)

                                     ccAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGTGG‑TTTTTGTTAAGCCGGGTTATTG‑GGAAAGATc                                                                                                               >  2:4562614/1‑150 (MQ=255)

                                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATCTGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACTATGTTGAACAGAGTTGGATCCCGTTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATcgc                                                                                                             <  2:2392843/150‑1 (MQ=255)

                                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATcgc                                                                                                             <  2:4206196/150‑1 (MQ=255)

                                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATcgc                                                                                                             <  2:1832443/150‑1 (MQ=255)

                                       aaaaCGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATcgc                                                                                                             <  2:4206502/150‑1 (MQ=255)

                                              cATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGCCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TGTTTGTTAAGCCGGGTGATTA‑GTGACTAGCGCGATTAga                                                                                                      >  2:3107713/1‑149 (MQ=255)

                                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGG‑TTTTTTTTTATCCGGGGTTTTCTGTAAA‑ATAGCGAttttcc                                                                                                     >  1:1389174/1‑145 (MQ=255)

                                               aTTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTCGGATCCCGGTCGTGTCGGGATTTTTTGTTAA‑CCGGGGTAATA‑GTAAATATAGCGCTGAGcc                                                                                                     >  2:1732587/1‑150 (MQ=255)

                                                    aTGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCgtttg                                                                                                <  2:4778738/150‑1 (MQ=255)

                                                    aTGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCgtttg                                                                                                <  2:3994849/150‑1 (MQ=255)

                                                    aTGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGGGTGGT‑TTTTTTTTTAGGCCGGTTTTTA‑GTAAATATCGCGATGAGCCctttt                                                                                                >  2:4067029/1‑149 (MQ=255)

                                                     tGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGG‑TTTTTTTTTATACCGGTTATTA‑GTAAAAATAGCGGTTAGCCggtttt                                                                                               >  1:4110861/1‑148 (MQ=255)

                                                      gCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCgtttgtt                                                                                              <  1:2780919/150‑1 (MQ=255)

                                                      gCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTATTGA‑TTTTTTTTAAGACGGGTTATTA‑GTAAAGATCGAGAAGAGCCctttttt                                                                                              >  1:3478225/1‑150 (MQ=255)

                                                      gCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGTGGG‑TTTTTTTTTATCCCGGGTATTA‑GTAAAGATAACGATGAGCCggttttt                                                                                              >  1:516243/1‑150 (MQ=255)

                                                        cGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTTGA‑TTTTTTTTAAGCCCGGTTATTC‑GTAAATATAGCGATGAGCCctttttttt                                                                                            >  2:624775/1‑149 (MQ=255)

                                                        cGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTTTTAAGCCGGGGGAGTC‑GTAAACATCGCGCAGAGCCCGGTGTTta                                                                                            >  1:2663081/1‑150 (MQ=255)

                                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAATCCGGGTTATTC‑GTAAATATCGCGATGAGCCGTTTGTTTATaaaa                                                                                       >  2:743617/1‑148 (MQ=255)

                                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACa                                                                                       <  1:486983/150‑1 (MQ=255)

                                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACa                                                                                       <  1:2210966/150‑1 (MQ=255)

                                                             ttCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACa                                                                                       <  2:3284839/150‑1 (MQ=255)

                                                                cGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTTTTAAGCCGGGGTATTA‑GGAAAGATAGCGATGGTCCGTTTTTTTTTAACAgcg                                                                                    >  1:3098461/1‑150 (MQ=255)

                                                                cGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGC‑TTTTTGGGAAGCCGGGTGATTC‑GTCAAGCTCGCGAGGAGCCGTTTGGTTATAACAgcg                                                                                    >  2:3558706/1‑150 (MQ=255)

                                                                  aagaaTATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAgcggc                                                                                  <  2:1537292/150‑1 (MQ=255)

                                                                       tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTc                                                                             <  2:1171241/150‑1 (MQ=255)

                                                                       tatGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTTTTAAGCCGGGTTATTC‑GTCAAGATCGCGATGAGCCGGTTTTTTATAAAAGCGGCGGTTc                                                                             >  2:1197182/1‑150 (MQ=255)

                                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTGTGGA‑TTTTTTTTAATACGGGTTTTTA‑GTAAAGATAGCGATGATGCGGTTGTTTTTAAAAACAGCGGttggaattt                                                                       >  2:3671174/1‑143 (MQ=255)

                                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCttt                                                                       >  2:3145877/1‑150 (MQ=255)

                                                                             tCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCgtt                                                                       >  2:3511260/1‑150 (MQ=255)

                                                                              cagcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCtttt                                                                      <  2:703556/150‑1 (MQ=255)

                                                                              cagcagAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCtttt                                                                      <  1:1694474/150‑1 (MQ=255)

                                                                                  agAACTGTCAGTCTAAGTATGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTGGTGA‑TTTTTTTTTTTCCGGGGTATTG‑GTTAAGATCGCGGTGAGACGTTTTTTTTTTAAAAGGGCGGTTGCTCATTTGGtt                                                                  >  1:2305604/1‑150 (MQ=255)

                                                                                     aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAATATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                               <  2:836691/150‑1 (MQ=255)

                                                                                     aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                               <  2:932716/150‑1 (MQ=255)

                                                                                     aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                               <  1:3562126/150‑1 (MQ=255)

                                                                                     aCTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTcgc                                                               <  2:1497828/150‑1 (MQ=255)

                                                                                            gTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGt                                                        <  2:4502458/150‑1 (MQ=255)

                                                                                               tAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTg                                                     <  1:1985880/150‑1 (MQ=255)

                                                                                                aaGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTgg                                                    <  2:2265894/150‑1 (MQ=255)

                                                                                                aaGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTgg                                                    <  1:566060/150‑1 (MQ=255)

                                                                                                    aCGAACATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCGGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                <  2:1726709/150‑1 (MQ=255)

                                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                <  1:538373/150‑1 (MQ=255)

                                                                                                    aCGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCgca                                                <  1:560000/150‑1 (MQ=255)

                                                                                                      gAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTGTCGGGG‑TTTTTGTTAAGCCGGGTTATTC‑GGAAAGATAGCGATGAGCCGGTTTTTTATAAAAGCGGCGGTTGCTCTTTTTCGTAGCGCCCCCTCGGGcgccgc                                              >  2:529793/1‑150 (MQ=255)

                                                                                                          aTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACtt                                          <  1:938813/150‑1 (MQ=255)

                                                                                                          aTTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACtt                                          <  2:3625934/150‑1 (MQ=255)

                                                                                                               tactacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGt                                     <  1:4776456/150‑1 (MQ=255)

                                                                                                                 ctacGACACTGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGa                                   <  2:4274082/150‑1 (MQ=255)

                                                                                                                  tacGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGAt                                  <  1:1806971/150‑1 (MQ=255)

                                                                                                                          tGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTg                          <  1:3722028/150‑1 (MQ=255)

                                                                                                                           gTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                         <  2:1044526/150‑1 (MQ=255)

                                                                                                                           gTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                         <  1:2997706/150‑1 (MQ=255)

                                                                                                                           gTTGAACAGAGTTGGATCCCGGTCGTGTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTgg                         <  2:2897012/150‑1 (MQ=255)

                                                                                                                             tGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTTTTAAGACGGGTTATTC‑GTAAATATCGCCAAGTGCCGTTTTTTTATAACAGCGGCGGTTACACTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTGCTGGGGATTAATCCAGgcg                       >  2:2329229/1‑150 (MQ=255)

                                                                                                                                 cAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTcc                   <  1:3362/150‑1 (MQ=255)

                                                                                                                                   gagTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCgg                 <  2:150946/150‑1 (MQ=255)

                                                                                                                                   gagTTGGATCCCGGTCGTGTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTACAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCgg                 <  2:150788/150‑1 (MQ=255)

                                                                                                                                            cccGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaa        <  1:87832/150‑1 (MQ=255)

                                                                                                                                                gTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCaaaaaa    <  1:1002527/150‑1 (MQ=255)

                                                                                                                                                  cGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAaga  <  2:2106980/150‑1 (MQ=255)

AATAATGCCGACATTTGCGCGCGTGGCTGGAACAACGCCAAAACGGCATTTTATGCCGATGTTCCGAAGAATATGACTCAGCAGAACTGTCAGTCTAAGTACGAAAATTGCTACTACGACAATGTTGAACAGAGTTGGATCCCGGTCGTTTCTGGA‑TTTTTGTTAAGCCGGGTTATTC‑GTAAAGATCGCGATGAGCCGTTTGTTTATAACAGCGGCGGTTCCTCTTTTGCTTCGCGCCCAGTCTGGCGCAGCACTTCTGGTGATTACTCCTGGCGCTCCGGTTCTGGCAAAAAAGA  >  NC_000913/4412630‑4412925

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: `‑`