BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,666,805	A→C	6.7%	S197R (AGT→CGT)	ydgD →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,666,805	0	A	C	6.7%	800.0 / 46.1	282	S197R (AGT→CGT)	ydgD	putative peptidase
Reads supporting (aligned to +/- strand): ref base A (145/118); new base C (6/13); total (151/131)
Fisher's exact test for biased strand distribution p-value = 5.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 287 total aligned reads displayed.

ATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGC  >  CP009273/1666657‑1666953

attgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACa                                                                                                                                                         >  2:3553065/1‑149 (MQ=255)

 ttgattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGAGAAAGCCGCGCTGACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAg                                                                                                                                                        <  2:3669554/149‑1 (MQ=255)

   gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtc                                                                                                                                                      >  1:2090594/1‑149 (MQ=255)

   gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtc                                                                                                                                                      >  1:156529/1‑149 (MQ=255)

     ttgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat                                                                                                                                                    >  1:3873618/1‑149 (MQ=255)

     ttgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat                                                                                                                                                    >  1:3093031/1‑149 (MQ=255)

     ttgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat                                                                                                                                                    >  2:229455/1‑149 (MQ=255)

     ttgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTAACGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcat                                                                                                                                                    >  1:68853/1‑149 (MQ=255)

      tgtgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcatc                                                                                                                                                   <  2:3374111/149‑1 (MQ=255)

       gtgCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTGATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtcatca                                                                                                                                                  >  1:1763381/1‑149 (MQ=255)

         gCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCaaa                                                                                                                                                >  1:1258251/1‑149 (MQ=255)

           tACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAAc                                                                                                                                              >  2:2120994/1‑149 (MQ=255)

                 aTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGaa                                                                                                                                        >  1:508800/1‑149 (MQ=255)

                  tCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAATCTGTGAAg                                                                                                                                       >  2:632891/1‑149 (MQ=255)

                   cccccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGt                                                                                                                                      >  2:178305/1‑149 (MQ=255)

                    ccccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTg                                                                                                                                     <  1:2630905/149‑1 (MQ=255)

                    ccccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTg                                                                                                                                     >  1:2572888/1‑149 (MQ=255)

                     cccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATACCGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGa                                                                                                                                    <  2:1411795/149‑1 (MQ=255)

                     cccTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGa                                                                                                                                    >  2:271696/1‑149 (MQ=255)

                            ggCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGAct‑‑‑ggctg                                                                                                                          >  2:1443741/1‑149 (MQ=255)

                             gCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGAct‑‑‑ggctgg                                                                                                                         >  2:610191/1‑149 (MQ=255)

                               aTTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGgc                                                                                                                       <  1:1504133/149‑1 (MQ=255)

                                ttACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGgcg                                                                                                                      >  2:2639601/1‑149 (MQ=255)

                                    gccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCaaa                                                                                                                  >  1:3921729/1‑149 (MQ=255)

                                     ccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAAc                                                                                                                 >  1:137458/1‑149 (MQ=255)

                                       gttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAAcct                                                                                                               >  1:745931/1‑147 (MQ=255)

                                           ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGt                                                                                                           <  2:516524/149‑1 (MQ=255)

                                           ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGt                                                                                                           >  2:1969592/1‑149 (MQ=255)

                                           ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGt                                                                                                           >  1:2771278/1‑149 (MQ=255)

                                                aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGt                                                                                                      <  2:713337/149‑1 (MQ=255)

                                                aTTTGAGGGAGAGAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGt                                                                                                      <  2:1039874/149‑1 (MQ=255)

                                                      gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATc                                                                                                >  2:3829388/1‑149 (MQ=255)

                                                      gggAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATc                                                                                                >  2:2424667/1‑149 (MQ=255)

                                                        gagaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAg                                                                                              >  2:1018473/1‑149 (MQ=255)

                                                          gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTg                                                                                            >  2:3439115/1‑149 (MQ=255)

                                                             aaaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGa                                                                                         <  2:1756371/149‑1 (MQ=255)

                                                              aaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGAt                                                                                        >  1:3661727/1‑149 (MQ=255)

                                                              aaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACAGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGAt                                                                                        >  1:1201715/1‑149 (MQ=255)

                                                                gCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATAc                                                                                      >  1:408712/1‑149 (MQ=255)

                                                                  cgcgCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCt                                                                                    >  2:3234565/1‑149 (MQ=255)

                                                                   gcgcTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATAc                                                                                      >  2:2782890/1‑149 (MQ=255)

                                                                           cgccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTg                                                                           >  2:1659446/1‑149 (MQ=255)

                                                                                   tAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGt                                                                   <  1:3486436/149‑1 (MQ=255)

                                                                                          ggCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCctc                                                            <  2:510844/149‑1 (MQ=255)

                                                                                             aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGa                                                         <  2:3534061/149‑1 (MQ=255)

                                                                                             aGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGa                                                         >  1:3403523/1‑149 (MQ=255)

                                                                                                tCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGt                                                      <  2:2138063/149‑1 (MQ=255)

                                                                                                 cGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGtt                                                     <  1:1615795/149‑1 (MQ=255)

                                                                                                 cGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGtt                                                     <  1:4069250/149‑1 (MQ=255)

                                                                                                   tAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGc                                                   >  1:440504/1‑149 (MQ=255)

                                                                                                    aaaGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCa                                                  <  2:2350422/149‑1 (MQ=255)

                                                                                                     aaGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCAt                                                 <  2:2910248/149‑1 (MQ=255)

                                                                                                     aaGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCAt                                                 >  2:1895290/1‑149 (MQ=255)

                                                                                                          aCTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGa                                            <  1:83178/149‑1 (MQ=255)

                                                                                                           cTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTTCAGTCATCAACACTGTGATGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGAt                                           <  1:3011517/149‑1 (MQ=255)

                                                                                                             caggcaggCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGa                                         <  1:2165026/149‑1 (MQ=255)

                                                                                                              aggcaggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGAc                                        <  2:2555670/149‑1 (MQ=255)

                                                                                                              aggcaggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGAc                                        <  2:3873444/149‑1 (MQ=255)

                                                                                                                gcaggcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACgg                                      <  2:3101938/149‑1 (MQ=255)

                                                                                                                    gcTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTgg                                  >  1:1956060/1‑149 (MQ=255)

                                                                                                                       aCCCTGAAGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaa                               <  2:1060963/149‑1 (MQ=255)

                                                                                                                       aCCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaa                               <  2:741550/149‑1 (MQ=255)

                                                                                                                       aCCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaa                               >  1:2405471/1‑149 (MQ=255)

                                                                                                                       aCCCTGAAGATCATCTCGATACGTTGTACAGACATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaa                               >  2:1090946/1‑149 (MQ=255)

                                                                                                                         ccTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCactt                             >  1:3115165/1‑149 (MQ=255)

                                                                                                                          cTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCacttc                            >  2:4015765/1‑148 (MQ=255)

                                                                                                                          cTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaatta                            >  2:3844527/1‑149 (MQ=255)

                                                                                                                            gAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaat                          >  2:86042/1‑149 (MQ=255)

                                                                                                                             aaGATCATCTCGATACGTTGTACAGTCATCAAATCTGTGAAGTGACT‑‑‑GGCTGTGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGGCTGGTGACAGCGGTTCGCCTCTGATGTTGCATACAGATGACGGATGGCaattaatt                         <  2:2237472/149‑1 (MQ=255)

                                                                                                                             aaGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCaattaatt                         >  1:877971/1‑149 (MQ=255)

                                                                                                                              aGATCATCTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTg                        >  2:2583723/1‑149 (MQ=255)

                                                                                                                                  catcTCGATACGTTGTACAGTCATCAAAACTGTTAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGt                    >  2:474638/1‑149 (MQ=255)

                                                                                                                                   atcTCGATACGTTGTACCGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTg                   <  1:1798403/149‑1 (MQ=255)

                                                                                                                                    tctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGc                  >  2:2302407/1‑149 (MQ=255)

                                                                                                                                    tctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGc                  <  1:1061369/149‑1 (MQ=255)

                                                                                                                                     ctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCa                 >  2:2370255/1‑149 (MQ=255)

                                                                                                                                     ctcGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCa                 >  1:1663254/1‑149 (MQ=255)

                                                                                                                                       cGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCaaa               >  2:1559127/1‑149 (MQ=255)

                                                                                                                                       cGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCaaa               <  2:3116124/149‑1 (MQ=255)

                                                                                                                                        gATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAg              >  2:121186/1‑149 (MQ=255)

                                                                                                                                        gATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAg              >  2:2200668/1‑149 (MQ=255)

                                                                                                                                         aTACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGt             >  1:1122228/1‑149 (MQ=255)

                                                                                                                                         aTACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGt             <  2:171722/149‑1 (MQ=255)

                                                                                                                                         aTACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGt             >  1:816974/1‑149 (MQ=255)

                                                                                                                                         aTACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGt             >  1:2004967/1‑149 (MQ=255)

                                                                                                                                         aTACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGt             <  1:399458/149‑1 (MQ=255)

                                                                                                                                         aTACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGt             <  2:172265/149‑1 (MQ=255)

                                                                                                                                          tACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGtt            <  2:3093031/149‑1 (MQ=255)

                                                                                                                                          tACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGtt            >  1:3183468/1‑149 (MQ=255)

                                                                                                                                          tACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGtt            >  2:2947749/1‑149 (MQ=255)

                                                                                                                                           aCGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTc           >  2:4055177/1‑149 (MQ=255)

                                                                                                                                           aCGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTc           >  2:3720843/1‑149 (MQ=255)

                                                                                                                                              ttGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc        <  1:3689723/149‑1 (MQ=255)

                                                                                                                                              ttGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc        >  1:730132/1‑149 (MQ=255)

                                                                                                                                              ttGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGc        <  1:962065/149‑1 (MQ=255)

                                                                                                                                                gTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTc      <  1:3553134/149‑1 (MQ=255)

                                                                                                                                                 tACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTcc     <  1:3234565/149‑1 (MQ=255)

                                                                                                                                                 tACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTcc     <  1:3769843/149‑1 (MQ=255)

                                                                                                                                                    aGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGc  >  2:3977034/1‑149 (MQ=255)

                                                                                                                                                    aGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGc  >  2:1483402/1‑149 (MQ=255)

ATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACT‑‑‑GGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTTCGCCTCTGATGTTGCATACCGATGACGGCTGGCAATTAATTGGGGTGCAAAGTTCGGCTCCTGC  >  CP009273/1666657‑1666953

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