BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,123,001	A→C	10.6%	Q247P (CAG→CCG)	yceM →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,123,001	0	A	C	10.6%	852.0 / 89.3	294	Q247P (CAG→CCG)	yceM	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base A (134/129); new base C (17/14); total (151/143)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 303 total aligned reads displayed.

TGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCG  >  CP009273/1122852‑1123150

 gCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAg                                                                                                                                                      <  2:3470860/150‑1 (MQ=255)

 gCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAg                                                                                                                                                      <  2:4926241/150‑1 (MQ=255)

 gCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAg                                                                                                                                                      <  2:618177/150‑1 (MQ=255)

  ccGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAgg                                                                                                                                                     <  2:1308288/150‑1 (MQ=255)

  ccGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAgg                                                                                                                                                     >  2:2925442/1‑150 (MQ=255)

    gAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGc                                                                                                                                                   <  1:3793732/150‑1 (MQ=255)

      gCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGt                                                                                                                                                 >  1:4607336/1‑150 (MQ=255)

        aCCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCCGGGCGTAg                                                                                                                                               <  1:1700120/150‑1 (MQ=255)

        aCCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAg                                                                                                                                               >  2:2867972/1‑150 (MQ=255)

         ccATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGt                                                                                                                                              >  1:1025299/1‑150 (MQ=255)

                 cGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCTGGCAGGGCGTAGTGCATAAAc                                                                                                                                      <  2:1980053/150‑1 (MQ=255)

                 cGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAAc                                                                                                                                      >  2:3790059/1‑150 (MQ=255)

                  ggCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAAcc                                                                                                                                     <  2:2753180/150‑1 (MQ=255)

                  ggCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAAcc                                                                                                                                     >  1:2979562/1‑150 (MQ=255)

                    cTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGa                                                                                                                                   >  1:4252989/1‑150 (MQ=255)

                      ggTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGAtt                                                                                                                                 <  1:4995956/150‑1 (MQ=255)

                       gTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTc                                                                                                                                >  1:643539/1‑150 (MQ=255)

                       gTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTc                                                                                                                                <  2:2606086/150‑1 (MQ=255)

                          cTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTg                                                                                                                             >  2:1594/1‑150 (MQ=255)

                          cTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTg                                                                                                                             >  2:4280796/1‑150 (MQ=255)

                             tGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCCGGGCGTAGTGCATAAACCGATTCCTGGtt                                                                                                                          <  1:1588455/150‑1 (MQ=255)

                                 gATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCa                                                                                                                      >  1:4014382/1‑150 (MQ=255)

                                  aTCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCag                                                                                                                     >  1:2477226/1‑150 (MQ=255)

                                   tCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCaga                                                                                                                    <  1:731323/150‑1 (MQ=255)

                                       caccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTAc                                                                                                                <  2:711527/150‑1 (MQ=255)

                                       caccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTAc                                                                                                                >  1:1333902/1‑150 (MQ=255)

                                         ccTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGc                                                                                                              <  2:1567352/150‑1 (MQ=255)

                                          cTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCt                                                                                                             <  1:5219089/150‑1 (MQ=255)

                                            gTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTg                                                                                                           >  2:1798544/1‑150 (MQ=255)

                                                gCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCa                                                                                                       <  1:2016909/150‑1 (MQ=255)

                                                  aTCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAAc                                                                                                     >  2:548206/1‑150 (MQ=255)

                                                     gCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTg                                                                                                  >  1:403414/1‑150 (MQ=255)

                                                     gCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTg                                                                                                  >  1:3484051/1‑150 (MQ=255)

                                                     gCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTg                                                                                                  >  1:1233466/1‑150 (MQ=255)

                                                        gTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGt                                                                                               <  1:5369708/150‑1 (MQ=255)

                                                         tGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGtt                                                                                              <  1:206433/150‑1 (MQ=255)

                                                         tGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGtt                                                                                              >  2:112324/1‑150 (MQ=255)

                                                            cGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGt                                                                                           <  1:4172304/150‑1 (MQ=255)

                                                                  tCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCtg                                                                                     <  1:3052385/150‑1 (MQ=255)

                                                                    aGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCCGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCtgtg                                                                                   >  2:5233573/1‑150 (MQ=255)

                                                                    aGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCtgtg                                                                                   <  1:4982718/150‑1 (MQ=255)

                                                                        tGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGcgcg                                                                               <  2:3535976/150‑1 (MQ=255)

                                                                            aCCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCCGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCAc                                                                           >  2:3150661/1‑150 (MQ=255)

                                                                            aCCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCAc                                                                           >  1:595454/1‑150 (MQ=255)

                                                                             ccGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACt                                                                          <  2:849200/150‑1 (MQ=255)

                                                                              cGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCCGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACtt                                                                         <  2:1498159/150‑1 (MQ=255)

                                                                              cGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACtt                                                                         <  2:4713465/150‑1 (MQ=255)

                                                                               gTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTc                                                                        >  2:4473027/1‑150 (MQ=255)

                                                                                  cAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCAtt                                                                     <  1:2265830/150‑1 (MQ=255)

                                                                                   aGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTg                                                                    <  1:4515969/150‑1 (MQ=255)

                                                                                     gCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGaa                                                                  <  2:122398/150‑1 (MQ=255)

                                                                                       cGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAAtg                                                                <  2:2439919/150‑1 (MQ=255)

                                                                                        gTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAAtgt                                                               <  2:2975797/150‑1 (MQ=255)

                                                                                         tgactgacGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCCGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAAtgtg                                                              <  2:2954699/150‑1 (MQ=255)

                                                                                           actgacGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAAtgtgtg                                                            <  2:1663251/150‑1 (MQ=255)

                                                                                            ctgacGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGc                                                           <  1:5511062/150‑1 (MQ=255)

                                                                                            ctgacGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGc                                                           >  1:140568/1‑150 (MQ=255)

                                                                                            ctgacGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGc                                                           <  1:3044677/150‑1 (MQ=255)

                                                                                                cGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCaaaa                                                       <  1:949480/150‑1 (MQ=255)

                                                                                                 ggTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAAc                                                      >  2:2057313/1‑150 (MQ=255)

                                                                                                  gTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAAcc                                                     >  2:4635335/1‑150 (MQ=255)

                                                                                                   tGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCa                                                    <  1:4810363/150‑1 (MQ=255)

                                                                                                   tGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCa                                                    >  1:3164173/1‑150 (MQ=255)

                                                                                                          aTCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGtt                                             >  2:908525/1‑150 (MQ=255)

                                                                                                          aTCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGtt                                             >  2:1053185/1‑150 (MQ=255)

                                                                                                          aTCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGtt                                             >  1:5438783/1‑150 (MQ=255)

                                                                                                             gACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCg                                          <  1:5342982/150‑1 (MQ=255)

                                                                                                               cATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCACAACCAGACAGTTCCGCa                                        >  2:2719943/1‑150 (MQ=255)

                                                                                                               cATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCa                                        >  1:4765112/1‑150 (MQ=255)

                                                                                                                aTTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaa                                       >  2:4256971/1‑150 (MQ=255)

                                                                                                                aTTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaa                                       <  1:1810816/150‑1 (MQ=255)

                                                                                                                 ttACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaaa                                      >  2:2272601/1‑150 (MQ=255)

                                                                                                                 ttACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCaaa                                      >  2:3335314/1‑150 (MQ=255)

                                                                                                                  tACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAAc                                     >  1:2305534/1‑150 (MQ=255)

                                                                                                                    cGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAAccg                                   <  1:1154405/150‑1 (MQ=255)

                                                                                                                    cGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAAccg                                   <  1:3281138/150‑1 (MQ=255)

                                                                                                                       atatGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAAccgccg                                <  1:2405583/150‑1 (MQ=255)

                                                                                                                        tatGCGCGAATGGCGTGAGGAGCGCGGGCCGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCgg                               >  1:5097402/1‑150 (MQ=255)

                                                                                                                        tatGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCgg                               <  2:2013589/150‑1 (MQ=255)

                                                                                                                        tatGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCgg                               >  2:3636553/1‑150 (MQ=255)

                                                                                                                        tatGCGCGAATGGCGAGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGAATCCTGGATGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGAGCGCGGCACTTCATTGAATGTGTGCAAAACCCGACAGTTCCGCAAACCGCCgg                               >  1:777774/1‑150 (MQ=255)

                                                                                                                           gcgcGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGa                            >  1:2537643/1‑150 (MQ=255)

                                                                                                                           gcgcGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGa                            >  1:706673/1‑150 (MQ=255)

                                                                                                                              cgAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACa                         >  1:3195524/1‑150 (MQ=255)

                                                                                                                                aaTGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAgg                       >  2:588486/1‑150 (MQ=255)

                                                                                                                                aaTGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAgg                       <  2:2199731/150‑1 (MQ=255)

                                                                                                                                 aTGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGc                      <  2:4445492/150‑1 (MQ=255)

                                                                                                                                  tGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGcc                     >  2:4526429/1‑150 (MQ=255)

                                                                                                                                     cGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTg                  >  1:2168203/1‑150 (MQ=255)

                                                                                                                                       tGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCt                >  2:3085252/1‑150 (MQ=255)

                                                                                                                                        gaggagCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTg               >  1:3369621/1‑150 (MQ=255)

                                                                                                                                         aggagCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTgg              <  1:4897830/150‑1 (MQ=255)

                                                                                                                                            agCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGgcgc           <  2:4119233/150‑1 (MQ=255)

                                                                                                                                            agCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGgcgc           <  2:116137/150‑1 (MQ=255)

                                                                                                                                                cgggccgggcGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACg       <  2:1963802/150‑1 (MQ=255)

                                                                                                                                                cgGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACg       <  1:5591470/150‑1 (MQ=255)

                                                                                                                                                  ggCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTa     <  1:4120589/150‑1 (MQ=255)

                                                                                                                                                     aGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCg  >  2:3769355/1‑150 (MQ=255)

                                                                                                                                                     aGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCg  >  1:2781883/1‑150 (MQ=255)

                                                                                                                                                     aGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCg  >  2:3009082/1‑150 (MQ=255)

TGCCGAGCACCATTTTTCGGCTGGTCCTTTGCAGATCACCACCTGTATGCATCGCCGTGCCGGAAGTCAGCGTGAAACCGTGCAGGCCGTGACTGACGGTGCGCTCATCGACATTACGGATATGCGCGAATGGCGTGAGGAGCGCGGGCAGGGCGTAGTGCATAAACCGATTCCTGGTTGGCAGAGTACGCTTGAGCAACGTGGGTTTGTCGGCTGTGCGCGGCACTTCATTGAATGTGTGCAAAACCAGACAGTTCCGCAAACCGCCGGCGAACAGGCCGTGCTGGCGCAACGTATCG  >  CP009273/1122852‑1123150

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