BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	576,183	G→T	8.8%	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
*	CP009273	576,183	0	G	T	8.8%	432.5 / 2.3	150	intergenic (‑282/‑107)	ylcI/nohD	DUF3950 family protein, DLP12 prophage/DLP12 prophage; DNA packaging protein
Reads supporting (aligned to +/- strand): ref base G (87/42); new base T (2/19); total (89/61)
Fisher's exact test for biased strand distribution p-value = 5.43e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

GTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTG  >  CP009273/576038‑576329

 tgACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtat                                                                                                                                                 <  1:815124/148‑1 (MQ=255)

    cATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtattta                                                                                                                                              <  1:2659261/148‑1 (MQ=255)

      tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                            >  2:807968/1‑148 (MQ=255)

      tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                            <  1:814794/148‑1 (MQ=255)

      tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                            >  2:317062/1‑148 (MQ=255)

      tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                            <  2:3266090/148‑1 (MQ=255)

      tGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGtatttatt                                                                                                                                            <  2:3266139/148‑1 (MQ=255)

          gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                        <  1:8922/148‑1 (MQ=255)

          gatgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTaaa                                                                                                                                        >  1:2177157/1‑148 (MQ=255)

           atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                       >  2:3155333/1‑148 (MQ=255)

           atgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTaaaa                                                                                                                                       >  1:2510401/1‑148 (MQ=255)

            tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                                      <  1:829074/148‑1 (MQ=255)

            tgatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTATTAGTCGTAACTTAATGTATTTATTTAAAAt                                                                                                                                      >  1:1498142/1‑148 (MQ=255)

             gatAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATa                                                                                                                                     <  1:2646447/148‑1 (MQ=255)

              atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                    <  2:3264059/148‑1 (MQ=255)

              atAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAc                                                                                                                                    <  2:3252394/148‑1 (MQ=255)

               tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                   >  1:117509/1‑148 (MQ=255)

               tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                   >  1:1732076/1‑148 (MQ=255)

               tAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcc                                                                                                                                   >  1:1891762/1‑148 (MQ=255)

                aaTTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAccc                                                                                                                                  >  2:784455/1‑148 (MQ=255)

                 attattATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAcccc                                                                                                                                 <  2:2270873/148‑1 (MQ=255)

                       atCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                           <  2:3348195/148‑1 (MQ=255)

                       atCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                           <  2:486470/148‑1 (MQ=255)

                       atCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGaaa                                                                                                                           <  2:486318/148‑1 (MQ=255)

                         cATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaag                                                                                                                         <  1:1408818/148‑1 (MQ=255)

                           ttttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaa                                                                                                                       >  1:2779042/1‑148 (MQ=255)

                            tttGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAaaagaaa                                                                                                                      >  2:1405804/1‑148 (MQ=255)

                                tgGTTCCTTTGCGGCGATCCGACAGTTTACGTGTCGGCGAGCTCGCTGGTTTTAGCTTTTTATTTAAATTTTCGGT‑TTTTA‑GGTGTTTCAG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAATAGAAAGTaa                                                                                                                  <  2:1421082/148‑1 (MQ=14)

                                tgGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaa                                                                                                                  <  2:2537774/148‑1 (MQ=40)

                                 gggTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGaaa                                                                                                                 >  2:522715/1‑148 (MQ=39)

                                     ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAc                                                                                                             >  1:3207198/1‑148 (MQ=39)

                                     ccTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGAc                                                                                                             >  1:3207683/1‑148 (MQ=39)

                                      cTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACa                                                                                                            >  1:1730771/1‑148 (MQ=39)

                                       tttCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAg                                                                                                           >  1:2827239/1‑148 (MQ=39)

                                                aTCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGaaa                                                                                                  >  2:3255603/1‑148 (MQ=37)

                                                 tCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGACACGACAGGTGCTGAAAg                                                                                                 >  2:573697/1‑148 (MQ=38)

                                                 tCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAg                                                                                                 >  2:573731/1‑148 (MQ=38)

                                                 cccGACGGGTTCGGGGGCGGCGACCTCGCGGTTTTTCGATATTTATTAAAATTTTCCGG‑TTTAC‑GGTGTTACCG‑TTCTTCTTCGTCGTAAATTATTTTATTTATTTAAAATCCCCCCTCAAAAGAAAGTAAACGACAGGTGCTGAAAg                                                                                                 <  2:3484435/147‑1 (MQ=255)

                                                          ttcgggggggGCGACCTCGCGGGTTTTCGCTATTTATAAAAATTTTCGGT‑TTTAT‑GTTGTTTCCT‑TTATGCTCCTTCTAAAATTATTGTATTTATTAAAAATACCCCCTCAAACGAATGTAAACGACAGGTGCTGAAAGCGAGCtttt                                                                                        <  2:2269300/144‑1 (MQ=9)

                                                                ggcggcGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCct                                                                                  >  1:2459201/1‑148 (MQ=255)

                                                                    gcgACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGt                                                                              >  2:3476654/1‑148 (MQ=255)

                                                                       aCCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATAACCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                           >  1:778615/1‑148 (MQ=255)

                                                                       aCCGCGGGGGTTTGCGCTATTTATGAAAATTTTCCGG‑TTTCA‑GGGGTTTCCG‑TTCTTCGTCGTCGTAACTTATTGTATTTATTTAAAATCCCCCCTGAAAACAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGt                                                                           <  2:1920536/148‑1 (MQ=14)

                                                                        ccTCGCGGGTTTTCGCTATTTATGAAAATTTTACGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGtt                                                                          >  1:1072455/1‑148 (MQ=255)

                                                                         cTCGCGGGTTTTCGCTATTTATTAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCACTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGttt                                                                         >  2:1297531/1‑148 (MQ=255)

                                                                            gcgGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATAAACCCTGAAAAGAAAGGAAACGACAGGTGCTGAAATCGAGCTTTTTGTCCTCTTTCGTTTcat                                                                      >  2:97358/1‑146 (MQ=14)

                                                                            gcgGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTcgt                                                                      >  1:692756/1‑146 (MQ=255)

                                                                             cgGGTTTGCTTTATTTAATAAATTTTTCCGT‑TTTAT‑GGTTTTTCTG‑TTGTACTTCGTCGTAACTTAATGTTTTTATTTAAAATCCCCCCCAAAAATAATGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCtt                                                                     <  2:1207596/148‑1 (MQ=2)

                                                                                gTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTtct                                                                  >  1:2600968/1‑148 (MQ=255)

                                                                                gTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCGTtct                                                                  >  1:2600976/1‑148 (MQ=255)

                                                                                gTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTAGTCGTAACTTAATGTATTTATTTAAAATACACCCTGAAAAGAAAGGAAAAGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTACTTtct                                                                  >  1:2688367/1‑148 (MQ=255)

                                                                                  tttCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGTTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                >  2:2356511/1‑148 (MQ=255)

                                                                                  tttCGATATTTATTAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCTTCGTAACATAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTctct                                                                <  2:581995/148‑1 (MQ=255)

                                                                                   ttctttttttATGAAATTTTTCGTT‑TTTAT‑GGTGTTCCGT‑TTTTTCTTCGACTTAATTTATTTTTTTAATTAAAAAAACCCCCCCAAAATAACCCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTg                                                               <  1:2427398/141‑1 (MQ=2)

                                                                                       cTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTGTCCTTTCTCTGtttt                                                           >  1:714541/1‑148 (MQ=255)

                                                                                         aTTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGtttttt                                                         >  1:2582837/1‑147 (MQ=255)

                                                                                         aTTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTg                                                         >  1:823834/1‑148 (MQ=255)

                                                                                         aTTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTg                                                         >  2:2514862/1‑148 (MQ=255)

                                                                                          tttATTAAATTTTTCGGG‑TTTAT‑GGTGTTTCCG‑TTCTTATTCGTCGTAAATTAATTTATTTATTTAAAATACCCCCTCAAAACAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                        <  2:1088270/148‑1 (MQ=14)

                                                                                          tttATGAAAATTTTCGGT‑TTTAT‑GGTGTTTCCT‑TTCTTCTTCGGCGTAACTAAATTTATTTATTTAAAAACCCCCCTCAAAATAACGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGt                                                        <  2:1018556/148‑1 (MQ=14)

                                                                                           ttattaAATTTTTCGGT‑TTTAT‑GGTGTTTCCT‑TTTTTTTTCGTCTTAATTTATTTTATTTATTTAAAATCCCCCGAGAAAATACTGGAAACGACAGGTGCTGAAAGCGAGCTTTTTTGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTc                                                       <  2:2234791/148‑1 (MQ=9)

                                                                                              tGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTGCTTTCTCTTTTTTTGTGCGt                                                    >  2:723861/1‑148 (MQ=255)

                                                                                              tGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                    >  2:2774158/1‑148 (MQ=255)

                                                                                              tGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCGTTCTCTGTTTTTGTCCGt                                                    >  2:3280678/1‑148 (MQ=255)

                                                                                               ttaaaTTTTTCGG‑TTTAATGGTTTTTCGT‑TTTTTCTTCGTCTTAATTTTATTTATTTATTAAAAATACCCCCTAAAAATAAAGCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGt                                                    <  1:1665208/146‑1 (MQ=14)

                                                                                                aaaaTTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACTACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTgg                                                  >  2:2311668/1‑148 (MQ=255)

                                                                                                  aaTTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTGCggggta                                                >  1:408306/1‑146 (MQ=255)

                                                                                                      tttcGG‑TTTTT‑GGTGTTTTCGTTTGTTCTCCGTATTAAATTATTTTATTTATTAAAAAAACCCCCTCAAAACAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgt                                             <  1:1291586/145‑1 (MQ=9)

                                                                                                      ttCCGGTTTTAT‑GGTGTTTCCG‑TTCTACTTCGTCTTAAATTAATGTTTTTATTAAAAACCCCCCCTCAAAACAAACCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgt                                             <  1:3469055/148‑1 (MQ=14)

                                                                                                      ttCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                            >  2:1557635/1‑148 (MQ=255)

                                                                                                      ttCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAAtgtg                                            >  2:1557891/1‑148 (MQ=255)

                                                                                                      ttCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTGTACTTTATCTGTTTTTGTGCGGGGGAtttg                                            >  1:1799508/1‑148 (MQ=14)

                                                                                                         ggg‑TTTAA‑GGTGTTTGCT‑TTGTTCTGGGTCGTAACTTATTGTTTATTTTTAAAATACCCCCTTAAAAGAACGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                         <  2:2984811/147‑1 (MQ=14)

                                                                                                         ctg‑TTTTA‑GTTGTTTCGT‑TTTTGCTTCGTCTTAACTTATTGTTTTTATTTAAAATACCCCCTTAAAACAATGAACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                         <  2:1237082/146‑1 (MQ=14)

                                                                                                         cGG‑TTTAT‑GTTGTTTCGT‑TTCTACTTGTACTTATATTATTGTTTTTTTTTAAAATACCCCATGAAAATAAAGGACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCaa                                         <  2:443067/148‑1 (MQ=14)

                                                                                                             tttAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATTTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTGCTTTGTATTTTTTTTTCCGTTGAATGTGCaaatg                                      >  1:558139/1‑145 (MQ=17)

                                                                                                             tttAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTTTGCGTGGCATGTGCAAttg                                      >  1:3364571/1‑146 (MQ=255)

                                                                                                                      gtTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTATCTGTTTTTTTGCGTGGTATGTTCCATTGAAGTgaac                              >  1:2229587/1‑148 (MQ=14)

                                                                                                                      gtTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTGTCCTTTCTATGTTTTTTTGGGTGGGATTTTCAATGGAAGTcaac                              >  1:2229713/1‑148 (MQ=14)

                                                                                                                        ttCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTGCTTTATCAGTTTTTGTGCGGGGGATGTTCAATGGGAGTcaaaaa                            >  1:2370093/1‑148 (MQ=14)

                                                                                                                          ccG‑TTCTTCTTCGTCTTAACTTATTGTTTTTATTTAAAATACCCCCTCAAAACAACGTAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACaaaa                          <  2:2779042/148‑1 (MQ=14)

                                                                                                                                 gtctgcGTCTAAACTTATTTTTTTTATTAAAAAAACCCCCTAAAAACACAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAagcagc                    <  1:2135483/147‑1 (MQ=14)

                                                                                                                                  tcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                   >  2:668063/1‑148 (MQ=255)

                                                                                                                                  tcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                   >  2:2655775/1‑148 (MQ=255)

                                                                                                                                  tcttcGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgct                   >  2:2655785/1‑148 (MQ=255)

                                                                                                                                    ttcGTCTTAACTAAATTTATTAATTTAAAATCCCCCCTCAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctgg                 <  2:356092/148‑1 (MQ=14)

                                                                                                                                       ttcttAACTTATTTTATTTATTAAAAACCCCCCCAAAAAAGAACGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg              <  1:422702/144‑1 (MQ=14)

                                                                                                                                       gtcgaAACTTAATTTATTTATTTAAAAACCCCCCTCAAAACAAAGCACACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAgctggctg              <  2:2124768/148‑1 (MQ=14)

                                                                                                                                         cgtAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGAc            >  2:1792526/1‑148 (MQ=255)

                                                                                                                                           tAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt          >  1:2908248/1‑148 (MQ=255)

                                                                                                                                           tAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCGTTCTCTGTTTTTGTGCGTGGAAAGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAt          >  2:1003066/1‑148 (MQ=255)

                                                                                                                                            aaattAATTTATTAATTAAAAACACCCCCTCAAAACAAACCAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtt         <  1:2335763/145‑1 (MQ=14)

                                                                                                                                            aaCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTGCGTGGTATGTTCAATTGAAGTCAACAAAAAGCAGCTGGCTGACAtt         >  1:1637892/1‑148 (MQ=255)

                                                                                                                                            aaCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtt         >  2:888301/1‑148 (MQ=255)

                                                                                                                                             aCTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATTGAATTCAACAAAAAGCAGCTTGCTGACAttt        >  2:1818745/1‑148 (MQ=255)

                                                                                                                                              cTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACCGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt       >  2:3403533/1‑148 (MQ=255)

                                                                                                                                              cTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt       >  2:1357073/1‑148 (MQ=255)

                                                                                                                                              cTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACAtttt       >  2:3104176/1‑148 (MQ=255)

                                                                                                                                                  aTGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGt   >  2:1746319/1‑148 (MQ=255)

                                                                                                                                                   tGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTGCGTTATCTGTTTTTGTGCGTGGTATGTTCAAATGCAGTCAACAAAAAGCAGCTGGCGGAAATTTTCGGTg  >  1:1843191/1‑148 (MQ=14)

GTGACATGTAGATGATAATTATTATCATTTTGTGGGTCCTTTCCGGCGATCCGACAGGTTACGGGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGG‑TTTAA‑GGTGTTTCCG‑TTCTTCTTCGTCGTAACTTAATGTATTTATTTAAAATACCCCCTGAAAAGAAAGGAAACGACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTGGAATGTGCAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTG  >  CP009273/576038‑576329

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