BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	580,727	G→T	9.6%	intergenic (‑282/‑107)	ylcI ← / → nohD	DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	580,727	0	G	T	9.6%	1423.2 / 23.6	469	intergenic (‑282/‑107)	ylcI/nohD	DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand): ref base G (228/176); new base T (3/60); total (232/237)
Fisher's exact test for biased strand distribution p-value = 2.16e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 474 total aligned reads displayed.

GTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580579‑580876

   gtgACATGTAGATGATAATTATGATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtatt                                                                                                                                                   <  2:1567398/150‑1 (MQ=255)

     gACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGtattta                                                                                                                                                 <  2:3676610/150‑1 (MQ=255)

         tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTAtttattta                                                                                                                                             <  1:5034431/150‑1 (MQ=255)

         tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTAtttattta                                                                                                                                             >  1:789739/1‑150 (MQ=255)

          gTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGGCAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaa                                                                                                                                            <  2:6291752/150‑1 (MQ=255)

           tAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                           <  2:2364276/150‑1 (MQ=255)

             gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                                         >  1:13125/1‑150 (MQ=255)

                 atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     <  1:7150324/150‑1 (MQ=255)

                 atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                     <  1:6477028/150‑1 (MQ=255)

                  tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcccc                                                                                                                                    <  2:4810317/150‑1 (MQ=255)

                  tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccac                                                                                                                                    >  1:7761784/1‑148 (MQ=255)

                   aaatattatAATGTTGTGGGTCCTTTCCGGCGATCAGACAGGTTACGGCTCGGCGACCTCGCGGGTATACGCTATATATGAAAATTTAACGGTTTAAGGTGATTCCGATCTTCTTCGTCGTAACTTAATGTATATATATAAAATAcaccc                                                                                                                                   >  2:8715651/4‑150 (MQ=14)

                    attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCt                                                                                                                                  >  1:800271/1‑150 (MQ=255)

                     ttattaTCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTg                                                                                                                                 >  2:3794291/1‑150 (MQ=255)

                      tattatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGa                                                                                                                                >  2:5840464/1‑150 (MQ=255)

                        ttatCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                              <  2:4431019/150‑1 (MQ=255)

                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGTGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCTTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:6224046/150‑1 (MQ=255)

                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:9916046/150‑1 (MQ=255)

                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:4353217/150‑1 (MQ=255)

                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  1:4427860/150‑1 (MQ=255)

                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:531500/150‑1 (MQ=255)

                             aTTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                         <  2:7062989/150‑1 (MQ=255)

                              ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              >  2:3310107/1‑128 (MQ=255)

                              ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTatttattt                                                                                                                                              <  1:3310107/128‑1 (MQ=255)

                              ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaag                                                                                                                        >  2:5696848/1‑150 (MQ=255)

                                  gtgGGTCCTTTCCGGCGATCCGACAGGTTTCGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                    <  1:4595617/150‑1 (MQ=255)

                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 <  2:2715713/150‑1 (MQ=39)

                                     ggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGa                                                                                                                 >  2:2535996/1‑150 (MQ=39)

                                         cTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAgg                                                                                                             >  2:1091440/1‑150 (MQ=39)

                                          tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGt                                                                                                            >  2:4187467/1‑150 (MQ=39)

                                           ttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTg                                                                                                           >  2:3496141/1‑150 (MQ=39)

                                               ggCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGa                                                                                                       >  1:10025386/1‑150 (MQ=39)

                                                   aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGc                                                                                                   >  2:1829583/1‑150 (MQ=38)

                                                     ccGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTACCGTTCTTCTTCGTCTTAACTTAATGTATTTATTTAAAATACCCCCTGAAACGAACGGAAACGACAGGTGCTGAAAGCGa                                                                                                 <  2:5079224/150‑1 (MQ=255)

                                                       gACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTGGTAACTTAATGTATTTATTTAATATACCCCCTGAAATGAAAGGAAACGACAGGTGCTGAAAGCGAGc                                                                                               >  2:981811/1‑150 (MQ=255)

                                                            gTTAGGGGGCGGCGACCTCGCGGGTTTTCGATATTTATTAAATTTTTCCGGTTTAAGGTGTTTCCTTTCTTCTTCTTCTAAACTTAATTTATTTATTTAAAATACCCCCTTAAAAGAACGCAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:6777372/150‑1 (MQ=14)

                                                            gTTACGGGGGGGCGACCTCGCGGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTTTTTCCGTACTTCTTCGTCGTAAATTAATTTATTTATTTAAAATACCCCCTCAAAACAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          <  2:9558783/150‑1 (MQ=11)

                                                            gTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCttttt                                                                                          >  1:6715710/1‑150 (MQ=40)

                                                             ttACGGGGCGGCGACCTCGCGGTTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAAATTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                                         <  2:8860275/150‑1 (MQ=255)

                                                             ttACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTg                                                                                         >  1:5296868/1‑150 (MQ=40)

                                                                cGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGcc                                                                                      >  2:7937622/1‑150 (MQ=255)

                                                                 ggggggGCTACCTGTGGGTTTTTCTCTATTTATTAATATTTTCTTGTTTATGGTGTTTCGTTTATTCTTGGTCGTAACTTAATGTATTTATTTATTTTACCCCCTGAAAAGAATGGAATCGACAGGTTCTGAAAGCGTGCTTTTTGGCct                                                                                     <  2:5578348/150‑1 (MQ=255)

                                                                  gggCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCctc                                                                                    >  1:6024246/1‑150 (MQ=255)

                                                                   ggcggcGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCcgct                                                                                   >  1:7008069/1‑150 (MQ=255)

                                                                       gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACACCCGGAAAAGAAATGAAACGACAGGTGCTGAAAGCGAGCATTTTGGCCTCTGTCg                                                                               >  1:6793898/1‑150 (MQ=255)

                                                                             tCGCGGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         <  2:3975227/150‑1 (MQ=255)

                                                                             tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCt                                                                         >  1:9223518/1‑150 (MQ=255)

                                                                             tCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTGTTTGGCCTCTGTCGTTTCCt                                                                         >  2:498967/1‑150 (MQ=255)

                                                                               gcgGGTTTTCGCTATTTATTAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       <  2:3164527/150‑1 (MQ=255)

                                                                               gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTTAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTTGCCTCTGTCGtttggttt                                                                       >  1:7582853/1‑145 (MQ=255)

                                                                               gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAATGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       >  1:2593694/1‑150 (MQ=255)

                                                                               gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       >  2:1483602/1‑150 (MQ=255)

                                                                               gcgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCttt                                                                       >  2:8901537/1‑150 (MQ=255)

                                                                                cgGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtc                                                                      >  2:7732668/1‑150 (MQ=255)

                                                                                 gggTTGTCGCTATTTATGAAAATTTCCCGGTTTATGGTGTTTCGGTGCTTCTTCGTCGTAACTTAATGTATTTATTAAAAATACCCCCTGAAAAGAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  2:9331491/150‑1 (MQ=14)

                                                                                 gggTATTCCCTATTTATTAAATTTTTCCGCTTTAAGCTTTTTCCGTTCTTCTTCGTCGTAACTTAATTTATTTATTAAAAATACCCCCTCAAAACAAACCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                     <  2:8476287/150‑1 (MQ=14)

                                                                                   gTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCAGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCATTctct                                                                   >  1:3717424/1‑150 (MQ=255)

                                                                                             tttATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                                                         <  2:8342556/150‑1 (MQ=255)

                                                                                             tttATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTcc                                                         >  1:4714706/1‑150 (MQ=255)

                                                                                               tATTAAATTTTTCGGGTTTATGGTGTTTCCTTTCTTCTTCGTCTTAACTTAATGTTTTTATTTAAAATACCCCGTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  2:6045713/150‑1 (MQ=14)

                                                                                               tATTAAAATTTTCCGGTTTATGGTGTTCCCTTTATTCTTCGTCGTAACTTAATTTATTTATTTAAAAAACCCCCTTAAAACAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                       <  1:492036/150‑1 (MQ=14)

                                                                                                  gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTCCGTGGa                                                    >  2:4607726/1‑150 (MQ=255)

                                                                                                  gAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGa                                                    >  1:92849/1‑150 (MQ=255)

                                                                                                   aaaaTTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGaa                                                   >  1:9552805/1‑150 (MQ=255)

                                                                                                        tttCCGGTTTATGTTTTTTCCTTTCTTATGCTTCTAAACTTAATTTATTTATTTAAAATCCCCCCTCAAAACAACGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGc                                              <  1:2246024/150‑1 (MQ=14)

                                                                                                          tCCTGTTTATTGTGTTTCCGTTCTTATTCGTCGTAACTTATTTTATTTATTTAAAATACCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                            <  2:1190495/150‑1 (MQ=14)

                                                                                                          tCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                            >  2:638718/1‑150 (MQ=255)

                                                                                                          tCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTATGTTTTTGTCCGTGGAATGTGCaa                                            >  1:2118233/1‑150 (MQ=255)

                                                                                                             gttttATGGTTTTTCGGTTCTTGTTCGTTTTAACTTAATTTTTTTATTTAAAACCCCCCGTTAAACGAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         <  2:4265387/148‑1 (MQ=14)

                                                                                                             ggTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGACTGTGCAATgg                                         >  2:3093544/1‑150 (MQ=255)

                                                                                                             ggTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATgg                                         >  2:9088387/1‑150 (MQ=255)

                                                                                                              gTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                        >  1:843873/1‑150 (MQ=255)

                                                                                                              cataaaGGTGTTTCCGTACTTCTTCGTCCAAAATAAATCAATTAATTAAAAATACCCCCTCAAAACAAAGCAAACGACAGGAGCTGAAAGCGAGATTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGa                                        <  2:6138251/146‑1 (MQ=17)

                                                                                                                  aaGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCGATGGAAGTc                                    >  1:8645929/1‑150 (MQ=255)

                                                                                                                     gtgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTGTTTGGCCTCTGTCGTTTCCTTTCTCTGTGTTTGTCCGGGGACTGTGCACGGGAAGTcaac                                 >  2:7583578/1‑150 (MQ=14)

                                                                                                                      tgtTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTcaaca                                >  1:731625/1‑150 (MQ=255)

                                                                                                                        tttCCTTTATTCTTCGCCGTAACTTATTTTTTTTTTTAAAAATACCCCCTTAAAACAACGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaa                              <  1:3040367/150‑1 (MQ=14)

                                                                                                                         ttCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             >  1:8870262/1‑150 (MQ=255)

                                                                                                                         tccctttcttcttcGTCGTAACTAAATGTATTTATTTAAAACACCCCCTCAAAAGAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                             <  1:170670/145‑1 (MQ=14)

                                                                                                                           ccGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAag                           >  2:7839746/1‑150 (MQ=255)

                                                                                                                            cGTTCTTCTTCGTCGTAACTTAATGTTTTTATTTAAAATACCCCCTGAAAACAACGGAATCGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagc                          <  1:4152697/150‑1 (MQ=255)

                                                                                                                             gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         <  2:7386300/150‑1 (MQ=255)

                                                                                                                             gTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagca                         >  1:5672700/1‑150 (MQ=255)

                                                                                                                               tcttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc                       <  2:540337/150‑1 (MQ=255)

                                                                                                                                ctacttcGTCGAAACTTAATTTTTTTATTAAAAACACCCCCTGAAAAGAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                      <  1:3295802/147‑1 (MQ=14)

                                                                                                                                 ttcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctg                     >  2:2792144/1‑150 (MQ=255)

                                                                                                                                    ttcGTCTTAATTTAATTTATTTATTAAAAATCCCCCCTGAAAATAATGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  <  1:5917456/150‑1 (MQ=14)

                                                                                                                                    ttcGTCTCAACTAAATGTATTATATAAAAACACCCCCTCAAAACAACCCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggct                  <  1:6904266/150‑1 (MQ=14)

                                                                                                                                     tcgtcgTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg                 >  1:3848703/1‑150 (MQ=255)

                                                                                                                                      cgtcttAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGa                <  2:143981/150‑1 (MQ=255)

                                                                                                                                        tcgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACa              <  2:718902/150‑1 (MQ=255)

                                                                                                                                         cttAACTTATTGTATTTATTTAAAATACCCCCTCAAAATAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt             <  2:6909678/148‑1 (MQ=255)

                                                                                                                                          ttAAATTATTGTATTTATTAAAAATACCCCATGAAAACACAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtt            <  2:2579500/149‑1 (MQ=14)

                                                                                                                                          gtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtt            >  2:9709752/1‑150 (MQ=255)

                                                                                                                                            aaCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt          <  2:2973241/150‑1 (MQ=255)

                                                                                                                                            aaCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt          <  2:6627039/150‑1 (MQ=255)

                                                                                                                                                tAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt      >  2:2023160/1‑150 (MQ=255)

                                                                                                                                                 aaTGTATTTATTTAAAATACCCCCTTAAAATAAAGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTg     <  2:4740061/150‑1 (MQ=255)

                                                                                                                                                  aTGTATTTATTTAAAATACCCCCTGAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGc    <  2:8084911/150‑1 (MQ=255)

                                                                                                                                                   tGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCGTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCg   >  1:821327/1‑150 (MQ=255)

GTTGTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGTGTTTCCGTTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGA  >  NC_000913/580579‑580876

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