BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	4,541,648	A→G	R62R (CGT→CGC)	RPPX_20280 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	4,541,648	0	A	G	98.4%	497.2 / ‑3.0	188	R62R (CGT→CGC)	RPPX_20280	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (98/87); total (98/90)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 190 total aligned reads displayed.

CGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTGCAGGCA  >  CP009974/4541500‑4541795

cGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcgcg                                                                                                                                                   <  1:1084478/151‑1 (MQ=255)

 ggATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCgg                                                                                                                                                  >  2:3207641/1‑151 (MQ=255)

 ggATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCgg                                                                                                                                                  >  1:428543/1‑151 (MQ=255)

   ataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAg                                                                                                                                                >  1:4566780/3‑151 (MQ=255)

   ataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAg                                                                                                                                                >  2:4915681/3‑151 (MQ=19)

     aGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAggg                                                                                                                                              <  2:1429083/151‑1 (MQ=255)

      ggTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgt                                                                                                                                             <  2:4882237/151‑1 (MQ=19)

      ggTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgt                                                                                                                                             <  2:3680138/151‑1 (MQ=255)

       gTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCGGCTTTTCGCCCGGGGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtg                                                                                                                                            <  2:2667388/151‑1 (MQ=11)

       gTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtg                                                                                                                                            <  1:395107/151‑1 (MQ=255)

        tAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtgt                                                                                                                                           <  1:2225131/151‑1 (MQ=255)

         aGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTa                                                                                                                                          <  2:3055976/151‑1 (MQ=255)

          gcgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAg                                                                                                                                         <  1:2610492/151‑1 (MQ=255)

           cgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAgg                                                                                                                                        <  1:4181806/151‑1 (MQ=255)

             cTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGcc                                                                                                                                      >  1:3780405/1‑151 (MQ=25)

             cTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGcc                                                                                                                                      >  1:3086241/1‑151 (MQ=21)

             cTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGcc                                                                                                                                      <  1:1914547/151‑1 (MQ=255)

              tCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCg                                                                                                                                     <  2:5046628/151‑1 (MQ=255)

               cACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGa                                                                                                                                    >  1:5004719/1‑151 (MQ=25)

                 ggggggtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGt                                                                                                                                  <  2:3841225/147‑1 (MQ=255)

                  ggcggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGACCACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGCGTAGGCCGAgcg                                                                                                                                 >  2:3363260/1‑149 (MQ=255)

                  ggcggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtg                                                                                                                                 >  1:1946949/1‑151 (MQ=255)

                  ggcggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtg                                                                                                                                 >  1:4382657/1‑151 (MQ=25)

                    cggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgct                                                                                                                               <  1:2056927/151‑1 (MQ=16)

                      gtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctgc                                                                                                                             <  2:3924809/151‑1 (MQ=19)

                      gtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctgc                                                                                                                             <  1:3092537/151‑1 (MQ=255)

                       tggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCCAAGCGCGGAGGGTGTAGGCCGAGtgctgcc                                                                                                                            >  1:4092450/1‑150 (MQ=255)

                       tggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCa                                                                                                                            >  2:2594878/1‑151 (MQ=25)

                        ggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAg                                                                                                                           >  2:1717874/1‑151 (MQ=25)

                        ggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAg                                                                                                                           >  2:1978564/1‑151 (MQ=255)

                          tAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCg                                                                                                                         >  1:2641641/1‑151 (MQ=255)

                            gTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTg                                                                                                                       >  1:422862/1‑151 (MQ=255)

                             tcctccCGCCCTTGGTAGCCAGGGTCTGTGTAGAGGCCCGCGCCGACGTGCTCCTGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCCTGCACACGTAGACATAGCCCGGGCGCACAGCACGGAGGGTGTAGGCCGCGTGCTGCAGCGTgg                                                                                                                      <  2:736590/151‑1 (MQ=11)

                             tcctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTgg                                                                                                                      >  2:2429408/1‑151 (MQ=255)

                             taatccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTgg                                                                                                                      <  1:2417516/148‑1 (MQ=255)

                              cctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTggg                                                                                                                     <  2:904259/151‑1 (MQ=21)

                                tccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTggggg                                                                                                                   >  2:3478206/1‑150 (MQ=255)

                                 ccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGaa                                                                                                                  <  1:5211222/151‑1 (MQ=255)

                                  cAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAg                                                                                                                 <  1:2820018/151‑1 (MQ=255)

                                  cAGCCCTTGGTAGCTCAGCTCTTTGTAGAGCCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAg                                                                                                                 <  2:3303687/151‑1 (MQ=255)

                                      ccTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                             >  1:1421340/1‑151 (MQ=255)

                                       cTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTg                                                                                                            >  1:265377/1‑151 (MQ=19)

                                       cTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTg                                                                                                            >  2:790824/1‑151 (MQ=19)

                                         tGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGtt                                                                                                          >  1:1745028/1‑151 (MQ=19)

                                           gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                                        <  1:3641076/151‑1 (MQ=12)

                                            tAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCa                                                                                                       >  2:4785719/1‑151 (MQ=9)

                                             agctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAg                                                                                                      <  1:1106360/151‑1 (MQ=9)

                                              gctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAgg                                                                                                     <  2:2827051/151‑1 (MQ=9)

                                               ctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGt                                                                                                    >  1:2287634/1‑151 (MQ=14)

                                               ctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGt                                                                                                    >  1:626862/1‑151 (MQ=9)

                                                                                          ccAGCTTTTCCCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCGGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTAGGGAAGCCTTGTTCCAGGTTGAAGCCTGCAACGGCATAGCGGCAGGCGGGGGCGTCGCCGGt                                                         <  1:3796052/151‑1 (MQ=255)

                                                                                           cAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTg                                                        >  2:2886024/1‑151 (MQ=9)

                                                                                            ggCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGc                                                       <  2:1869481/150‑1 (MQ=9)

                                                                                             gCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCg                                                      <  1:4733793/151‑1 (MQ=14)

                                                                                             gCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCg                                                      >  1:3254351/1‑151 (MQ=9)

                                                                                              cTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgg                                                     >  1:1331547/1‑151 (MQ=9)

                                                                                                tttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgggg                                                   >  1:4812259/1‑151 (MQ=12)

                                                                                                tttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgggg                                                   <  2:1182136/151‑1 (MQ=17)

                                                                                                tttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgggg                                                   >  2:615869/1‑151 (MQ=12)

                                                                                                  tCGCCCTGAGGGGCGTGCAGGAACACGTAGACAGAGCCTGGGCGCAGAGCGCGGAGGGTGTAGGCCGAGTGCTGCCGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCa                                                 <  2:4695644/151‑1 (MQ=11)

                                                                                                  tCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCa                                                 >  1:1608176/1‑151 (MQ=9)

                                                                                                  tCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCa                                                 <  2:2216195/151‑1 (MQ=17)

                                                                                                   cGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCAc                                                <  1:1244713/151‑1 (MQ=17)

                                                                                                     cccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGa                                              >  1:2239443/1‑151 (MQ=12)

                                                                                                      ccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGAt                                             >  2:2374979/1‑151 (MQ=9)

                                                                                                       cTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATg                                            <  2:5132115/151‑1 (MQ=12)

                                                                                                        tGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCGTCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATgg                                           >  2:4578038/1‑151 (MQ=9)

                                                                                                        tGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATgg                                           <  2:1596233/151‑1 (MQ=12)

                                                                                                        tGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATgg                                           >  1:555878/1‑151 (MQ=12)

                                                                                                         gAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                                          >  1:1633694/1‑151 (MQ=9)

                                                                                                         gAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                                          <  1:1015801/151‑1 (MQ=12)

                                                                                                         aaGGGCCTTTCATGAACACGTAGACATAGCCTGGGCACAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                                          <  1:1015781/150‑1 (MQ=12)

                                                                                                          aGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCg                                         >  1:2656860/1‑151 (MQ=12)

                                                                                                           gggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGt                                        <  1:390661/151‑1 (MQ=12)

                                                                                                            ggCCTTTCGTGAACACGTAGACATAGCCTGGGGGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTa                                       <  2:2739986/151‑1 (MQ=255)

                                                                                                            ggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTa                                       <  2:4992614/151‑1 (MQ=17)

                                                                                                            ggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTa                                       >  1:1879764/1‑151 (MQ=9)

                                                                                                                tttCAGGAACACGTAGACATCGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgcgc                                   <  1:1074794/151‑1 (MQ=9)

                                                                                                                   cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAcc                                <  1:4461510/151‑1 (MQ=18)

                                                                                                                    aTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCg                               >  1:2294836/1‑151 (MQ=16)

                                                                                                                     tGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCTTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCgg                              >  2:3142515/1‑151 (MQ=255)

                                                                                                                     tGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCgg                              <  2:1268347/151‑1 (MQ=21)

                                                                                                                     tGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCgg                              <  2:1268708/151‑1 (MQ=21)

                                                                                                                      gAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGc                             <  1:3107791/151‑1 (MQ=19)

                                                                                                                         cacGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAgg                          >  1:1772390/1‑151 (MQ=19)

                                                                                                                           cGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGAt                        <  2:3599119/151‑1 (MQ=25)

                                                                                                                           cGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGAt                        >  2:2346662/1‑151 (MQ=19)

                                                                                                                            gTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATa                       >  1:298248/1‑151 (MQ=21)

                                                                                                                                 cATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAc                  >  1:1978816/1‑151 (MQ=255)

                                                                                                                                 cATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAc                  >  2:4872536/1‑151 (MQ=255)

                                                                                                                                  aTAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcg                 >  1:683435/1‑151 (MQ=255)

                                                                                                                                         gggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAATCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTgcgc          >  1:4041445/1‑151 (MQ=255)

                                                                                                                                          ggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCt         <  1:1662650/151‑1 (MQ=38)

                                                                                                                                            cgcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgc       >  1:581801/1‑151 (MQ=37)

                                                                                                                                            cgcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgc       <  1:790824/151‑1 (MQ=38)

                                                                                                                                            cgcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgc       <  1:3409895/151‑1 (MQ=255)

                                                                                                                                             gcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgca      >  1:3013279/1‑151 (MQ=255)

                                                                                                                                             gcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgca      >  2:3148749/1‑151 (MQ=37)

                                                                                                                                              cAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgcag     <  1:3310119/151‑1 (MQ=255)

                                                                                                                                                 aGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgcaggca  <  1:3578413/151‑1 (MQ=255)

CGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTGCAGGCA  >  CP009974/4541500‑4541795

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