BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,385,516	A→T	5.5%	intergenic (‑129/+301)	tldD ← / ← yhdP	putative peptidase/DUF3971‑AsmA2 domains protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,385,516	0	A	T	5.5%	710.3 / 2.1	279	intergenic (‑129/+301)	tldD/yhdP	putative peptidase/DUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand): ref base A (126/123); new base T (3/23); total (129/150)
Fisher's exact test for biased strand distribution p-value = 1.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

CTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTT  >  CP009273/3385368‑3385665

  tACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGc                                                                                                                                                    <  1:1252020/150‑1 (MQ=255)

  tACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGc                                                                                                                                                    <  2:3175849/150‑1 (MQ=255)

   aCCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGGTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCa                                                                                                                                                   >  2:1483006/1‑150 (MQ=255)

     cAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAAc                                                                                                                                                 <  2:198150/150‑1 (MQ=255)

      aGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGATAGTAGGCCTGATAACACCCACAGCGTCCCAACAGGCAAcc                                                                                                                                                >  2:4509732/1‑149 (MQ=37)

          tAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAATACGCACAGCGTCGCAGCAGGCAACGGCTg                                                                                                                                            >  2:2511453/1‑150 (MQ=255)

          tAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                                                                                                            <  2:4465151/150‑1 (MQ=255)

           aaGACTCATCGTTTTTGCTACTCGTTAGTTACTGCCGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTga                                                                                                                                           >  2:4374548/1‑149 (MQ=255)

           aaGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGt                                                                                                                                           >  2:4114443/1‑150 (MQ=255)

            aGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGATAGTAGGCCTTATAAAACGAAAATCCTCTCATAATGCAACgcataac                                                                                                                                          >  1:1670940/1‑144 (MQ=25)

              aCTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCgg                                                                                                                                        >  1:2868951/1‑150 (MQ=255)

               cTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTTGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGGCGGa                                                                                                                                       >  1:962501/1‑150 (MQ=255)

               cTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGa                                                                                                                                       >  1:3178963/1‑150 (MQ=255)

               cTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGa                                                                                                                                       >  1:3642554/1‑150 (MQ=255)

                  aTCGTTTTTGCTACTCGTTAGTTACTGCATTAGAGTATTTTTTACGGCTGCCGGATGCGGCGTTAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                                    <  2:1704407/150‑1 (MQ=255)

                  aTCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                                    >  2:1762696/1‑150 (MQ=255)

                  aTCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATgc                                                                                                                                    >  1:46972/1‑150 (MQ=255)

                       ttttGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGt                                                                                                                               <  2:1559229/150‑1 (MQ=255)

                        tttGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTa                                                                                                                              <  2:667498/150‑1 (MQ=255)

                        tttGCTACTCGTTAGTTACGGCAGTAGCGGTTTTTTTCCGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTa                                                                                                                              <  1:1055762/150‑1 (MQ=255)

                          tGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTaaa                                                                                                                            >  1:3930462/1‑150 (MQ=255)

                           gCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAAc                                                                                                                           <  2:3156085/150‑1 (MQ=255)

                             tACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAATACGCACAGCGACGCATCAGGCAACGCCTGACGGATGCGGCGTAAACGc                                                                                                                         >  1:3504456/1‑150 (MQ=255)

                             tACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGc                                                                                                                         >  1:3504408/1‑150 (MQ=255)

                              aCTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGTCTGTCGGATGCGGCGTAAACGcc                                                                                                                        >  2:429933/1‑150 (MQ=255)

                              aCTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCAGGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGATGACGGATGCGGCGTAAACGcc                                                                                                                        >  2:2752001/1‑150 (MQ=255)

                               cTCGTTAGTTACTGCAGTAGTGGATTTTTTCCGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTCGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCt                                                                                                                       <  1:246610/150‑1 (MQ=255)

                                tCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCtt                                                                                                                      <  2:3457431/150‑1 (MQ=255)

                                tCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAAACGAAAAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGccat                                                                                                                      >  2:3370779/1‑148 (MQ=255)

                                  gttagttaCTGCTGTATAGGATTTTTTACGGCTGCCGGATGCGGCGTTAACGCCATAACCAGCCTACGGTTATGTACCGGTTTGTAGGCCTTATAATAGGCACCGCTTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTAt                                                                                                                    <  1:2756080/150‑1 (MQ=18)

                                    tagttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATcc                                                                                                                  <  2:3974697/150‑1 (MQ=255)

                                       ttaCTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCCGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGAc                                                                                                               <  1:1048263/150‑1 (MQ=255)

                                          cTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTa                                                                                                            >  1:2391884/1‑150 (MQ=255)

                                           tGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                                                                           <  2:3359409/150‑1 (MQ=255)

                                           ggCAGTAGTGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCATATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAATCCGCACCGCGTCGCAACAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTAc                                                                                                           <  2:28055/149‑1 (MQ=18)

                                              tgTAGAGGATTTTTTCGGGCTGCCGGATGCGGCGTTAACGCCCTATCCATCCTAGTGTTATTTTCGGTTTTTTAGGCCTGATATGACGCACCGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt                                                                                                        <  1:2132971/149‑1 (MQ=17)

                                              agatgtggATTTTTCACGGCTGCCGGAGGCGGCGTGAACGCCCTATCCTCCCTACGTTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACCGCGGCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt                                                                                                        <  1:1469121/144‑1 (MQ=17)

                                               gTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt                                                                                                       <  2:3798328/150‑1 (MQ=255)

                                                 agagGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTTATAAGACTCACAGCGTCGCATCAGGCAACGGCTGTCGGTTGCGGCGTAAACGCGTTATCCGACCTACGGTTag                                                                                                     >  1:628029/1‑149 (MQ=255)

                                                    ggATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGtt                                                                                                  >  1:2683964/1‑150 (MQ=255)

                                                     gATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTc                                                                                                 >  1:374448/1‑150 (MQ=255)

                                                      aTTTTTTACGGCTGCCGGATGCGGCGTGAACCCCCTATCCAGCCTACGTTTATGTTCCGGTTTGTTGGCCTTATAAGACGCACTGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTcc                                                                                                <  1:1333146/150‑1 (MQ=17)

                                                       ttttttACGGGTGCCGGTTGGGGGGTTAACGCCCTATCCCGCCTTAGGTTTTTTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGCCGGATGCGGCGTGAACGCCCTATCCGACCTACGGTTATGTTctg                                                                                               <  2:3628948/150‑3 (MQ=11)

                                                        tttttCCGGCGGCCGGATGCGGCGTGAACGCCCTACCCAGCCTACGGTTATGTTCCGGTTTTTCGGCATGATAAGACGCACCGCGTCGCATCCGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGt                                                                                              <  1:213119/150‑1 (MQ=17)

                                                         ttttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAATCGCCTTATCCGACCTACGGTTATGTTCCGtt                                                                                             >  2:2117222/1‑150 (MQ=255)

                                                         ttttACGGATGCCGGATGCGGCGTGACCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt                                                                                             <  2:665290/150‑1 (MQ=255)

                                                          tttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCGGTTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGttt                                                                                            <  1:1317477/150‑1 (MQ=255)

                                                          tttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATACGACGCACAGCGTCGCATCAGGCAACGCCTGCCGGATGCGGCGTAAACGCCGAATCCGACCTACGGTTCTGTGCCggat                                                                                            >  2:703495/1‑147 (MQ=14)

                                                           ttACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCAGGTTTGTAGGCCTGATAAGACGAACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACATACGGTTATGTTCCGGTTg                                                                                           >  1:1587871/1‑150 (MQ=17)

                                                            tACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGt                                                                                          <  1:4591245/150‑1 (MQ=255)

                                                              cGGCTGCCGGATGCGGCGTGAACGCCCTATCCCGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGTGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAg                                                                                        <  1:3440641/150‑1 (MQ=255)

                                                                 cTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGcc                                                                                     >  1:2569471/1‑150 (MQ=255)

                                                                   gCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAACGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTg                                                                                   <  1:3527022/150‑1 (MQ=255)

                                                                     cGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGTTTTGTAGGCCTGATAAGACGCACCGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGAt                                                                                 <  1:450963/150‑1 (MQ=255)

                                                                      ggATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATa                                                                                <  1:3758512/150‑1 (MQ=255)

                                                                          ggggcgTTAACGCCCTTGCCTGGCTACGTTTATGTTCCGGTTTGTCGGCCTGATAATACGCCAAGCGTCGCATTAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                            <  1:3836238/148‑1 (MQ=14)

                                                                          gcggcgTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGAc                                                                            <  2:353588/150‑1 (MQ=255)

                                                                               gTGACCGCCCTACCCCGCCTCCGGTTTTGTACCGGTTTGTTGGCCTTATAAGCCGCACAGCGTCGCATCTGCCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                       <  1:12646/150‑1 (MQ=14)

                                                                               gTGAACGCCCTATCCAGCCTACGGTTTTGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                       <  1:1085514/150‑1 (MQ=255)

                                                                               gTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGcaca                                                                       <  2:2365539/150‑1 (MQ=255)

                                                                                tGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAg                                                                      >  1:4075263/1‑150 (MQ=255)

                                                                                 gAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGc                                                                     >  2:270613/1‑150 (MQ=255)

                                                                                  aaCGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGACGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGTTTGTAGGCCTGATAAGACGCACAGCg                                                                    >  2:4091419/1‑150 (MQ=255)

                                                                                      cccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGc                                                                <  2:1178917/150‑1 (MQ=255)

                                                                                       ccTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCGGGATGTTGGCCTGATAAGACGCACAGCGTCGCa                                                               >  1:4245019/1‑150 (MQ=255)

                                                                                         tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                             <  2:305709/150‑1 (MQ=255)

                                                                                         tATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCAGGTTGTAGGCCTGATAAGACGCACAGCGTCGCATc                                                             >  1:228129/1‑150 (MQ=255)

                                                                                          aTCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCAAAGCGTCGCATCa                                                            >  2:262991/1‑150 (MQ=255)

                                                                                                 cTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACg                                                     <  2:1537299/150‑1 (MQ=37)

                                                                                                 cTACGGTTATGTTCCGGTTGGTAGGCCTGATAAGACGCACAGCGTCGCCTCAGGCCACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACg                                                     <  2:1255173/150‑1 (MQ=255)

                                                                                                  tACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACgg                                                    >  2:1550464/1‑150 (MQ=37)

                                                                                                    cGGTTTTGTTCCGGTTTGTTGGCCTTATAAGACGCACAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCt                                                  <  1:3793597/150‑1 (MQ=255)

                                                                                                    cGGTTATGTTCCGTTTTGTAGGCCTGATAAGACGCCCAGCGTCGCATCTGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCt                                                  <  1:1918524/150‑1 (MQ=255)

                                                                                                     ggTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTg                                                 >  2:2017450/1‑150 (MQ=37)

                                                                                                           gTTCCGTTTTGTTGGCCTGTTAAGACGCACTGGTTCGCATCTGGCAACGGCTGTCGGATGCGGCTTAAACGCCTTATCCGACCTACGTTTATGTTCCGTTTGTAGGCCTTATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                           <  2:3302873/150‑1 (MQ=11)

                                                                                                           gTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGAt                                           >  2:93462/1‑150 (MQ=37)

                                                                                                                ggTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggc                                      >  1:671444/1‑150 (MQ=32)

                                                                                                                 ttttGTCGGCCTTATAAGACGCACCGCGTCCCATCTGGAAACGGCTGTCGGTTGCGGCTTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATgcggcg                                     <  2:3887086/149‑1 (MQ=11)

                                                                                                                 ttttGTCGGCCTGATACGCCGCACCGCGTCGCATCCGGCAACGGCTGTCGTTTGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTCATAAGACGCACCGCGTCGCATCAGGCAACGGCTGCCGGATgcggcg                                     <  2:81467/149‑1 (MQ=11)

                                                                                                                    tGTTGGCATGATACGACGCACCGCGTCGCAACAGGAAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGAATCAGGCAACGGCTGCCGGATGCGGCGTaa                                  <  2:2436430/150‑1 (MQ=11)

                                                                                                                      tAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAAc                                >  2:2942937/1‑150 (MQ=34)

                                                                                                                         gCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGAGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGcc                             >  2:1630537/1‑150 (MQ=21)

                                                                                                                           cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                           >  2:2962213/1‑150 (MQ=35)

                                                                                                                           cTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCtt                           >  1:185055/1‑150 (MQ=35)

                                                                                                                                aaGCCGCCCAGCGTCGCATCTGGCAACGGCGGTCGGATGCGGCTTAAACGCCTGATCCGACCTACGGTTAGGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCAACAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCg                      <  2:1794084/150‑1 (MQ=11)

                                                                                                                                 aGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGa                     >  2:2808214/1‑150 (MQ=37)

                                                                                                                                  gACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCGTATCCGAc                    >  2:908847/1‑150 (MQ=21)

                                                                                                                                      cacaGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTAc                >  1:4414029/1‑150 (MQ=37)

                                                                                                                                        caGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGGAAACGCCTTATCCGACCTACgg              >  1:914234/1‑150 (MQ=37)

                                                                                                                                         aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt             <  2:710324/150‑1 (MQ=37)

                                                                                                                                         aGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCCCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCAGACCTACggg             >  2:46530/1‑149 (MQ=255)

                                                                                                                                         aGCGCCGCATCTGGCAACGGCTGTCGGATGCGGCTTACACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTCGGCCTGATAAGACGCCCAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGt             <  2:4347108/150‑1 (MQ=11)

                                                                                                                                          gcgcCGCATATGTCACCGGCTGTGGTTTGCGGCTTAAACGCCTTACCCGACCTCCTTTTATTTTCCGTTTGTCGGCCTTATAAGACGCACAGCGTCGAATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt            <  1:279282/150‑1 (MQ=11)

                                                                                                                                          gCGTCGCATCTGTCACCGGCGGTGGGATGCGGCGTAACCGCCTAATCCGACCTCCGGTTATGTTCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGtt            <  2:883436/150‑1 (MQ=11)

                                                                                                                                            gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt          >  1:1358967/1‑150 (MQ=37)

                                                                                                                                            gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt          >  1:3678930/1‑150 (MQ=37)

                                                                                                                                            gTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTAt          >  1:3011830/1‑150 (MQ=37)

                                                                                                                                              cGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGt        <  2:2196656/150‑1 (MQ=37)

                                                                                                                                               gCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCCTATCCGACCTACGGTTATGtt       >  1:4177747/1‑150 (MQ=17)

                                                                                                                                                    aGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGtt  >  2:1339690/1‑150 (MQ=40)

CTTACCAGGTTAAGACTCATCGTTTTTGCTACTCGTTAGTTACTGCAGTAGAGGATTTTTTACGGCTGCCGGATGCGGCGTGAACGCCCTATCCAGCCTACGGTTATGTTCCGGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGTCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTTTGTAGGCCTGATAAGACGCACAGCGTCGCATCAGGCAACGGCTGCCGGATGCGGCGTAAACGCCTTATCCGACCTACGGTTATGTTCCGTT  >  CP009273/3385368‑3385665

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