BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,457,499	C→A	90.3%	R20S (CGT→AGT)	mgtA →	magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,457,499	0	C	A	90.3%	830.3 / 82.3	307	R20S (CGT→AGT)	mgtA	magnesium transporter
Reads supporting (aligned to +/- strand): ref base C (14/16); new base A (137/140); total (151/156)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 307 total aligned reads displayed.

TCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGA  >  CP009273/4457353‑4457644

tCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGt                                                                                                                                                 >  1:4617256/1‑149 (MQ=255)

tCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGt                                                                                                                                                 >  2:826638/1‑149 (MQ=255)

     gCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATcc                                                                                                                                            >  2:2651315/1‑149 (MQ=255)

      cTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCa                                                                                                                                           <  2:3013099/149‑1 (MQ=255)

       tCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAt                                                                                                                                          >  2:318425/1‑149 (MQ=255)

       tCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAt                                                                                                                                          <  1:1900943/149‑1 (MQ=255)

        cAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCAtt                                                                                                                                         <  1:2003938/149‑1 (MQ=255)

           cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGcc                                                                                                                                      <  1:2031491/149‑1 (MQ=255)

           cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGcc                                                                                                                                      <  2:1964161/149‑1 (MQ=255)

           cTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGcc                                                                                                                                      <  2:1053247/149‑1 (MQ=255)

               attaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgc                                                                                                                                  <  2:3141988/149‑1 (MQ=255)

                ttaCCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcg                                                                                                                                 >  2:4702134/1‑149 (MQ=255)

                   ccTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGgcgcgc                                                                                                                              >  1:1650395/1‑149 (MQ=255)

                         ggTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGa                                                                                                                        >  1:1975341/1‑149 (MQ=255)

                            aaGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCGTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCCGCAGACAg                                                                                                                     >  1:2580714/1‑149 (MQ=255)

                             aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGt                                                                                                                    >  1:3432930/1‑149 (MQ=255)

                             aGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGt                                                                                                                    <  2:2304720/149‑1 (MQ=255)

                               gCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGa                                                                                                                  >  2:3348691/1‑149 (MQ=255)

                                  tCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATa                                                                                                               >  1:3071901/1‑149 (MQ=255)

                                      cACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggt                                                                                                           >  1:4124581/1‑149 (MQ=255)

                                      cACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggt                                                                                                           >  1:3496617/1‑149 (MQ=255)

                                      cACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACggt                                                                                                           <  1:1013022/149‑1 (MQ=255)

                                          ccTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTc                                                                                                       <  2:2686102/149‑1 (MQ=255)

                                               aGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCccgcc                                                                                                  >  1:1784996/1‑149 (MQ=255)

                                                 aCATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCAATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGt                                                                                                <  1:590582/149‑1 (MQ=255)

                                                        tGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                         >  2:2519652/1‑149 (MQ=255)

                                                        tGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                         >  1:4700170/1‑149 (MQ=255)

                                                        tGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAg                                                                                         >  1:3379085/1‑149 (MQ=255)

                                                           aCTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTgc                                                                                      <  2:2471237/149‑1 (MQ=255)

                                                              gTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                   <  1:2918047/149‑1 (MQ=255)

                                                              gTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCa                                                                                   <  1:44361/149‑1 (MQ=255)

                                                               tttCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAt                                                                                  >  2:678624/1‑149 (MQ=255)

                                                                    gACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCt                                                                             >  1:1570720/1‑149 (MQ=255)

                                                                      cAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGa                                                                           <  1:2225563/149‑1 (MQ=255)

                                                                       aGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGACCCGCCGTCCTTAAGTGCGCATTGCCTGaa                                                                          >  2:2824770/1‑149 (MQ=255)

                                                                        gTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaa                                                                         >  1:4290761/1‑149 (MQ=255)

                                                                         tGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGaaaa                                                                        <  2:1679263/149‑1 (MQ=255)

                                                                          gCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAAt                                                                       <  1:4077382/149‑1 (MQ=255)

                                                                          gCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAAt                                                                       <  1:4371709/149‑1 (MQ=255)

                                                                           cGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATg                                                                      <  2:2788479/149‑1 (MQ=255)

                                                                            ggAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATgg                                                                     >  2:1631083/1‑149 (MQ=255)

                                                                                ggACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGt                                                                 <  2:1780054/149‑1 (MQ=255)

                                                                                  aCTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGa                                                               <  2:440011/149‑1 (MQ=255)

                                                                                    tCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATg                                                             >  2:2015425/1‑149 (MQ=255)

                                                                                    tCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATg                                                             >  1:3844054/1‑149 (MQ=255)

                                                                                    tCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATg                                                             >  1:691100/1‑149 (MQ=255)

                                                                                      cTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGcc                                                           <  1:897567/149‑1 (MQ=255)

                                                                                       ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGccc                                                          <  1:2480077/149‑1 (MQ=255)

                                                                                       ttATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGccc                                                          <  1:3756474/149‑1 (MQ=255)

                                                                                        tATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCg                                                         >  1:2173932/1‑149 (MQ=255)

                                                                                          tGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTTCCTGAAAATGGCGGTGATGCCCgaa                                                       <  2:3721428/149‑1 (MQ=255)

                                                                                              tAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                   >  2:2110821/1‑149 (MQ=255)

                                                                                              tAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                   >  2:3704661/1‑149 (MQ=255)

                                                                                              tAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaag                                                   <  1:2656058/149‑1 (MQ=255)

                                                                                               aaaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaaga                                                  <  1:826638/149‑1 (MQ=255)

                                                                                                aaGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCgaagaagaa                                                 >  2:730206/1‑149 (MQ=255)

                                                                                                 aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                                <  2:3878264/149‑1 (MQ=255)

                                                                                                 aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                                <  2:4628901/149‑1 (MQ=255)

                                                                                                 aGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAt                                                >  2:3630385/1‑149 (MQ=255)

                                                                                                  gAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAAtt                                               >  1:3959854/1‑149 (MQ=255)

                                                                                                   aaaTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtg                                              <  1:512179/149‑1 (MQ=255)

                                                                                                     aTTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATtgtg                                            >  2:1536141/1‑149 (MQ=255)

                                                                                                      tttttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTgg                                           <  1:1917362/149‑1 (MQ=255)

                                                                                                      tttttACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTgg                                           <  1:601448/149‑1 (MQ=255)

                                                                                                          tACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaaaa                                       <  1:1399847/149‑1 (MQ=255)

                                                                                                          tACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGaaaa                                       >  2:3398079/1‑149 (MQ=255)

                                                                                                           aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAAc                                      <  1:4415080/149‑1 (MQ=255)

                                                                                                           aCCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAAc                                      <  2:1289453/149‑1 (MQ=255)

                                                                                                              cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGtt                                   <  1:3921664/149‑1 (MQ=255)

                                                                                                              cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGtt                                   <  1:4038934/149‑1 (MQ=255)

                                                                                                              cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGtt                                   >  2:4532257/1‑149 (MQ=255)

                                                                                                              cGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGtt                                   >  2:4573016/1‑149 (MQ=255)

                                                                                                                gCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCg                                 >  2:2555745/1‑149 (MQ=255)

                                                                                                                gCTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCg                                 >  2:3243688/1‑149 (MQ=255)

                                                                                                                 cTCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGa                                >  2:2946718/1‑149 (MQ=255)

                                                                                                                 cTCATTATCCATTTACATTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCCTTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGa                                >  2:3363224/1‑149 (MQ=255)

                                                                                                                  tCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGac                               <  1:2775095/149‑1 (MQ=255)

                                                                                                                  tCATTCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGac                               >  1:3131114/1‑149 (MQ=255)

                                                                                                                      tCGCCATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATGGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACAccc                           <  1:3657247/149‑1 (MQ=255)

                                                                                                                         ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATc                        >  1:3725754/1‑149 (MQ=255)

                                                                                                                         ccATTTACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATc                        <  1:3438716/149‑1 (MQ=255)

                                                                                                                             ttACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGa                    >  1:652001/1‑149 (MQ=255)

                                                                                                                             ttACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGa                    >  2:4202813/1‑149 (MQ=255)

                                                                                                                              tACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGaa                   >  1:2269855/1‑149 (MQ=255)

                                                                                                                              tACCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGaa                   <  2:1478379/149‑1 (MQ=255)

                                                                                                                               aCCTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAg                  >  1:3429810/1‑149 (MQ=255)

                                                                                                                                ccTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAgg                 <  1:1260367/149‑1 (MQ=255)

                                                                                                                                 cTTCCCGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAggg                >  2:3844164/1‑149 (MQ=255)

                                                                                                                                    cccGTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTa             >  1:1514494/1‑149 (MQ=255)

                                                                                                                                       gTCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAAt          >  2:3448556/1‑149 (MQ=255)

                                                                                                                                        tCTGGTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATc         <  2:1071964/149‑1 (MQ=255)

                                                                                                                                           ggTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg      <  1:755368/149‑1 (MQ=255)

                                                                                                                                           ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg      <  2:2438951/149‑1 (MQ=255)

                                                                                                                                           ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg      >  2:1332387/1‑149 (MQ=255)

                                                                                                                                           ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg      >  1:4019447/1‑149 (MQ=255)

                                                                                                                                           ggTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAgg      >  2:888756/1‑149 (MQ=255)

                                                                                                                                            gTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCAGTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGc     <  2:4124949/149‑1 (MQ=255)

                                                                                                                                            gTTCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGc     >  1:1328628/1‑149 (MQ=255)

                                                                                                                                             ttCATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCg    <  2:494363/149‑1 (MQ=255)

                                                                                                                                               cATAGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGa  <  2:2865465/149‑1 (MQ=255)

TCCCTGCTCAGCTTTATTACCTTCAGGTAAGGCTTCGCCACGCCTGAAGACATTTCTGTACTGTTTCAGACAGTGCGGAGGGACTCCTTATGTTTAAAGAAATTTTTACCCGGCTCATTCGCCATTTACCTTCCCGTCTGGTTCATCGTGATCCATTGCCTGGCGCGCAGCAGACAGTGAATACGGTGGTCCCGCCGTCCTTAAGTGCGCATTGCCTGAAAATGGCGGTGATGCCCGAAGAAGAATTGTGGAAAACGTTCGACACCCATCCGGAAGGGTTAAATCAGGCGGA  >  CP009273/4457353‑4457644

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