BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	4,540,460	C→T	R458R (CGG→CGA)	RPPX_20280 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	4,540,460	0	C	T	99.1%	716.7 / ‑5.1	225	R458R (CGG→CGA)	RPPX_20280	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (107/116); total (108/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 225 total aligned reads displayed.

CCCCGCCAATCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAACGCCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAGGGGTTATCGATGTGTTG  >  CP009974/4540313‑4540609

 cccGCCAATCCCTCGATCATGTTGTCGCTGGCGCCCAACAAGGAGCCGATGGCGTCGGATTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGACGCGTCGGt                                                                                                                                                   <  1:525140/151‑1 (MQ=17)

    gCCAATCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCAt                                                                                                                                                >  2:2460511/1‑151 (MQ=255)

    gCCAATCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCAt                                                                                                                                                <  1:4630140/151‑1 (MQ=255)

     ccAATCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATg                                                                                                                                               <  2:2312316/151‑1 (MQ=255)

     ccAATCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATg                                                                                                                                               <  1:2398243/151‑1 (MQ=255)

       aaTCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGAt                                                                                                                                             >  2:4916324/1‑151 (MQ=255)

           ccTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTc                                                                                                                                         <  2:3609756/151‑1 (MQ=255)

               gATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTg                                                                                                                                     >  1:3292178/1‑151 (MQ=255)

                 tCACGTTGTCGCTGGCGCCGAACGAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGcc                                                                                                                                   <  2:608928/151‑1 (MQ=255)

                 ggaggTTGTCGGTGGGGGGGAGCAAGGGGCGGATGGCGTCGGCGTTGTCCTTGAGCGGGTTGAGCGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGcc                                                                                                                                   <  2:1318671/147‑1 (MQ=14)

                   aCGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCgg                                                                                                                                 >  2:1527056/1‑151 (MQ=255)

                       tGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTGTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGgtgt                                                                                                                             <  1:3612139/151‑1 (MQ=255)

                        gTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTc                                                                                                                            >  1:5021804/1‑151 (MQ=255)

                           gCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGc                                                                                                                         <  2:4396199/151‑1 (MQ=255)

                            cTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCt                                                                                                                        <  1:740643/151‑1 (MQ=255)

                             gggCGCGGAACAAGGGGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGGGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTg                                                                                                                       <  2:5072906/150‑1 (MQ=14)

                               gcgcCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTCAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGc                                                                                                                     <  1:260248/151‑1 (MQ=255)

                                  ccGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGcccc                                                                                                                  >  2:638319/1‑150 (MQ=255)

                                  ccGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCa                                                                                                                  >  2:3720670/1‑151 (MQ=255)

                                  ccGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCa                                                                                                                  <  1:3342065/151‑1 (MQ=255)

                                   cGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCAt                                                                                                                 >  1:2269546/1‑151 (MQ=255)

                                     gaCAAGGAGCCGAGGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTg                                                                                                               <  2:2115172/150‑1 (MQ=255)

                                     aaCAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTg                                                                                                               >  1:1503011/1‑151 (MQ=255)

                                        aaGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGt                                                                                                            <  2:2960906/151‑1 (MQ=255)

                                        aaGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGt                                                                                                            <  1:330213/151‑1 (MQ=255)

                                         aGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTg                                                                                                           <  1:5088521/151‑1 (MQ=255)

                                         aGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTg                                                                                                           <  1:4116088/151‑1 (MQ=255)

                                         aGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTg                                                                                                           >  1:965113/1‑151 (MQ=255)

                                          ggAGCGGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGa                                                                                                          <  2:2605002/151‑1 (MQ=255)

                                          ggAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGa                                                                                                          <  1:4487040/151‑1 (MQ=255)

                                          ggAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGa                                                                                                          <  2:2202149/151‑1 (MQ=255)

                                            aGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAAt                                                                                                        >  2:1860715/1‑151 (MQ=255)

                                             gCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATc                                                                                                       <  1:4692239/151‑1 (MQ=255)

                                                gATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGc                                                                                                    >  1:3604924/1‑151 (MQ=255)

                                                gATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGc                                                                                                    <  1:1720105/151‑1 (MQ=255)

                                                  tGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCaa                                                                                                  <  2:2363870/151‑1 (MQ=255)

                                                    gCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACg                                                                                                >  1:1377836/1‑151 (MQ=21)

                                                       gcGGCTTTGTCCTTGAACGGGTTGAGCGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCa                                                                                             <  2:3681768/150‑1 (MQ=255)

                                                        cGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTGGGCTGGCCCATTGGGTGAATCAGCCACGCGAg                                                                                            >  1:3701031/1‑151 (MQ=255)

                                                            tttGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTgag                                                                                        >  2:930020/1‑151 (MQ=12)

                                                            tttGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTgag                                                                                        >  1:442430/1‑151 (MQ=12)

                                                            tttGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTgag                                                                                        >  1:30775/1‑151 (MQ=12)

                                                             ttGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGCCGGTCATGATCTTCGGTGCCGGGAGTGGGCTGGCCGATTGGAGGAATAAGCCACGCTAGCgaga                                                                                       >  2:3153425/1‑151 (MQ=11)

                                                              tGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGTCAGTGAGAt                                                                                      <  1:3129875/151‑1 (MQ=9)

                                                               gTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGAtt                                                                                     <  2:3523965/151‑1 (MQ=12)

                                                               gTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATATTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGAtt                                                                                     <  1:2524842/151‑1 (MQ=9)

                                                                  cTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTcc                                                                                  <  2:3993594/151‑1 (MQ=12)

                                                                  cTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTcc                                                                                  >  1:4981144/1‑151 (MQ=12)

                                                                    tGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACTCCAGTGAGATTTCCAg                                                                                >  1:34152/1‑151 (MQ=12)

                                                                    tGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAg                                                                                <  1:2786960/151‑1 (MQ=12)

                                                                      aaCGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCAAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGAc                                                                              <  1:2978973/151‑1 (MQ=9)

                                                                       agtggTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATAATCTTCGGAGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGAcc                                                                             <  2:1650114/148‑1 (MQ=255)

                                                                       aCGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTGGTGCTGGGCCAGCCATGGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGAcc                                                                             <  2:1114619/151‑1 (MQ=255)

                                                                        cGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACcg                                                                            <  1:1911346/151‑1 (MQ=16)

                                                                           gTTGAGCGGCGCCAAGGCGTTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCgcgg                                                                         <  1:1931919/151‑1 (MQ=255)

                                                                           gTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGCCCgcgg                                                                         >  2:1613993/1‑151 (MQ=255)

                                                                           gTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCgcgg                                                                         <  1:2867210/151‑1 (MQ=19)

                                                                            ttGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCgcggg                                                                        >  1:3924479/1‑151 (MQ=19)

                                                                               aaCGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCgcgggcgg                                                                     <  2:1476133/151‑1 (MQ=255)

                                                                               aaCGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCgcgggcgg                                                                     <  1:3950162/151‑1 (MQ=25)

                                                                                  ggCGCCAGGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCGAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAgg                                                                  <  2:4961846/151‑1 (MQ=255)

                                                                                    cgcCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTg                                                                >  1:4344811/1‑151 (MQ=18)

                                                                                    cgcCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTg                                                                >  1:1199041/1‑151 (MQ=19)

                                                                                       cAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGt                                                             >  1:3864507/1‑151 (MQ=19)

                                                                                        aaGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGtt                                                            >  1:3231487/1‑151 (MQ=255)

                                                                                         aGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTc                                                           >  2:4685869/1‑151 (MQ=21)

                                                                                           gCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGt                                                         >  2:1174261/1‑151 (MQ=255)

                                                                                             ggTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCa                                                       <  1:628510/151‑1 (MQ=255)

                                                                                              gtgtAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCAc                                                      >  1:3653314/1‑151 (MQ=255)

                                                                                              gtgtAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCAc                                                      >  1:4429777/1‑151 (MQ=255)

                                                                                              gtgtAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCAc                                                      <  2:1356944/151‑1 (MQ=255)

                                                                                                   aaGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCg                                                 >  1:1615392/1‑151 (MQ=255)

                                                                                                     gCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTa                                               <  2:556249/151‑1 (MQ=255)

                                                                                                      cGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGCGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCCGAGGTGCGTTCGCCACggcccttc                                              >  2:2515835/1‑145 (MQ=11)

                                                                                                       gggCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCTTTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACTGCGGGCGGAGGTGCGTGCGTCACGGTCGGAAc                                             >  1:1273174/1‑151 (MQ=11)

                                                                                                       gggCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAAc                                             >  2:523182/1‑151 (MQ=255)

                                                                                                        ggCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACa                                            >  1:4187949/1‑151 (MQ=255)

                                                                                                        ggCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACa                                            <  2:3304424/151‑1 (MQ=255)

                                                                                                         gCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAg                                           <  1:1797813/151‑1 (MQ=255)

                                                                                                           tGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGc                                         >  1:3208952/1‑151 (MQ=25)

                                                                                                            gTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGcc                                        <  2:2158074/151‑1 (MQ=255)

                                                                                                            gTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGcc                                        >  1:5215872/1‑151 (MQ=255)

                                                                                                               gTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAg                                     >  2:1602133/1‑151 (MQ=255)

                                                                                                               gTGCTGGTCCAGCCATGGCTCGCGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAg                                     <  2:4488551/151‑1 (MQ=14)

                                                                                                                 gCTGGTCGAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAggc                                   <  2:4274221/151‑1 (MQ=255)

                                                                                                                  cTGGTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAggcg                                  <  1:3983073/151‑1 (MQ=255)

                                                                                                                     gTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCt                               <  1:107541/151‑1 (MQ=255)

                                                                                                                     gTCCAGCCATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCt                               <  1:1527056/151‑1 (MQ=255)

                                                                                                                           ccATTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAgg                         >  2:812797/1‑151 (MQ=255)

                                                                                                                             aTTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCt                       <  2:4836868/151‑1 (MQ=255)

                                                                                                                             aTTGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCt                       <  2:1400719/151‑1 (MQ=255)

                                                                                                                               tGCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGt                     >  1:157423/1‑151 (MQ=255)

                                                                                                                                gCTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTa                    <  2:2580482/151‑1 (MQ=255)

                                                                                                                                 cTCCAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAg                   >  2:1912362/1‑151 (MQ=255)

                                                                                                                                   ccAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAggg                 >  2:1973627/1‑151 (MQ=255)

                                                                                                                                    cAGGCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAgggg                >  1:4962876/1‑151 (MQ=255)

                                                                                                                                      ggCGTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCGATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAGGGGtt              <  2:5015321/151‑1 (MQ=255)

                                                                                                                                         gTTTGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAGGGGTTATc           >  1:4850858/1‑151 (MQ=255)

                                                                                                                                          tttGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAGGGGTTATCg          <  1:2720464/151‑1 (MQ=255)

                                                                                                                                            tGAAGCGTCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAGGGGTTATCGAt        <  2:3186852/151‑1 (MQ=255)

CCCCGCCAATCCCTCGATCACGTTGTCGCTGGCGCCGAACAAGGAGCCGATGGCGTCGGCTTTGTCCTTGAACGGGTTGAACGGCGCCAAGGCGGTGTAAAGCGGGCTGTCGTGCTGGTCCAGCCATTGCTCCAGGCGTTTGAAACGCCGGTCATGATCTTCGGTGCCGGGTGTCGGCTGGCCCATTGGGTGAATCAGCAACGCCAGTGAGATTTCCAGACCGCGGGCGGAGGTGCGTTCGTCACGGTCGTAACAGGCCAAGGCGGCTGCCAGGCTGTAGGGGTTATCGATGTGTTG  >  CP009974/4540313‑4540609

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