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	99.5%	517.3 / ‑4.6	188	R62R (CGT→CGC)	RPPX_20280	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (95/92); total (95/93)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 189 total aligned reads displayed.

CGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTGCAGGC  >  CP009974/4541500‑4541794

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

  gataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGa                                                                                                                                                >  2:2062224/4‑151 (MQ=255)

   ataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAACGCGCGGAg                                                                                                                                               >  2:4368569/3‑151 (MQ=255)

    taGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAgg                                                                                                                                              >  1:5715070/2‑151 (MQ=255)

    taGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAgg                                                                                                                                              <  1:824423/150‑1 (MQ=19)

        tAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtgt                                                                                                                                          <  1:904391/151‑1 (MQ=21)

        tAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtgt                                                                                                                                          >  1:5436969/1‑151 (MQ=255)

          gcgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAg                                                                                                                                        <  1:186683/151‑1 (MQ=25)

            gcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGc                                                                                                                                      >  2:5587764/1‑151 (MQ=25)

            gcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGc                                                                                                                                      <  1:1356643/151‑1 (MQ=19)

              tCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCg                                                                                                                                    <  2:3202687/151‑1 (MQ=25)

                    cggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgct                                                                                                                              <  1:2992956/151‑1 (MQ=25)

                       gggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCa                                                                                                                           <  2:5507920/150‑1 (MQ=255)

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

                         gtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGc                                                                                                                         <  1:2088164/151‑1 (MQ=255)

                          tAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCg                                                                                                                        <  1:3294438/151‑1 (MQ=19)

                          tAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCg                                                                                                                        <  1:4107569/151‑1 (MQ=255)

                           aGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGTGTGCTGCAGCGt                                                                                                                       >  2:2752573/1‑151 (MQ=255)

                           aGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGt                                                                                                                       <  2:5715070/151‑1 (MQ=255)

                            gTCCTCGAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCGTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTg                                                                                                                      <  1:5964162/151‑1 (MQ=255)

                               ctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTgggg                                                                                                                   >  1:432784/1‑151 (MQ=255)

                                 ccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGaa                                                                                                                 <  2:5901918/151‑1 (MQ=19)

                                    gCCCTTGGGAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGcc                                                                                                              <  1:2106294/151‑1 (MQ=255)

                                      ccTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                            <  1:711746/151‑1 (MQ=16)

                                      ccTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                            <  1:3132270/151‑1 (MQ=18)

                                      ccTTGGGAGCTCAGCTCTTTGTACAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATCGTCTGGGCGCAAAGCGAGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                            <  2:1371129/151‑1 (MQ=11)

                                       cTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTg                                                                                                           <  1:3464455/151‑1 (MQ=16)

                                        ttGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGt                                                                                                          >  1:869628/1‑151 (MQ=21)

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

                                          ggTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTc                                                                                                        >  2:2543191/1‑151 (MQ=16)

                                          ggTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTc                                                                                                        >  2:308960/1‑151 (MQ=16)

                                            tAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGa                                                                                                                                   <  2:4461935/122‑1 (MQ=16)

                                            tAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGGGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCa                                                                                                      >  1:5579398/1‑151 (MQ=9)

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

                                              gctccgctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAgg                                                                                                    <  2:2991763/151‑1 (MQ=9)

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

                                               ctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGt                                                                                                   <  2:6095450/151‑1 (MQ=14)

                                                tcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGtt                                                                                                  <  1:1238077/151‑1 (MQ=2)

                                                                     cggcGACGTGCTCATGAATGACCAGCTGTTCGCCCTGAGGGCCTTTGATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCg                                                                             <  2:6329414/148‑1 (MQ=0)

                                                                                           cAGCTTTTCGCCCTGAGGGCCTTTCATGCACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCCGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGGCGCCGGTg                                                       >  2:4111997/1‑151 (MQ=11)

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

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

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

                                                                                              cgtttCGCGCTGGGGGCCGTTCCTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGGGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgg                                                    <  1:3849453/149‑1 (MQ=11)

                                                                                              cTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgg                                                    >  2:683294/1‑151 (MQ=14)

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

                                                                                               ttttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCggg                                                   <  1:867437/151‑1 (MQ=17)

                                                                                                tttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCgggg                                                  >  2:2153441/1‑151 (MQ=9)

                                                                                                 ttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGc                                                 <  2:1904889/151‑1 (MQ=17)

                                                                                                  tCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCa                                                >  1:287500/1‑151 (MQ=12)

                                                                                                    gCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACg                                              >  1:2310005/1‑151 (MQ=9)

                                                                                                    gCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACg                                              <  2:6139305/151‑1 (MQ=17)

                                                                                                     cccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGa                                             <  2:3787884/151‑1 (MQ=17)

                                                                                                       cTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATg                                           >  1:3566452/1‑151 (MQ=12)

                                                                                                       cTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATg                                           >  1:3378223/1‑151 (MQ=9)

                                                                                                       cTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCAGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATg                                           >  1:4422204/1‑151 (MQ=9)

                                                                                                        tGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATgg                                          >  2:2205755/1‑151 (MQ=9)

                                                                                                          aGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCg                                        >  2:4187563/1‑151 (MQ=9)

                                                                                                          aGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCg                                        >  1:950031/1‑151 (MQ=14)

                                                                                                          aGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCg                                        <  1:466906/151‑1 (MQ=12)

                                                                                                            ggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTa                                      <  1:768600/151‑1 (MQ=12)

                                                                                                            ggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTa                                      <  1:4557611/151‑1 (MQ=12)

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

                                                                                                             gcgtttCATGAACGCGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGGTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAg                                     <  2:4608768/148‑1 (MQ=11)

                                                                                                             gCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAg                                     <  1:787084/151‑1 (MQ=12)

                                                                                                              ccTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgc                                    >  1:5662544/1‑151 (MQ=9)

                                                                                                              ccTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgc                                    <  2:3654223/151‑1 (MQ=17)

                                                                                                               cTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgcg                                   >  2:2891522/1‑151 (MQ=14)

                                                                                                                    gtGAACACGGAGACCTGGGCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGCGTGCTGCAGCGTGGGGAGGCCTTGTTACTGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCg                              <  2:2290640/150‑1 (MQ=11)

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

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

                                                                                                                      gAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGc                            >  2:3043698/1‑151 (MQ=16)

                                                                                                                       aaCACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTATGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGAGGGGCAAGATGGCGTAGCGCActcgaa                           >  2:5590875/1‑146 (MQ=11)

                                                                                                                        acacGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAg                          >  1:446184/1‑151 (MQ=255)

                                                                                                                          acGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGa                        >  1:1589899/1‑151 (MQ=19)

                                                                                                                          acGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGa                        >  1:3260525/1‑151 (MQ=255)

                                                                                                                           cGTAGACATAGCCTGGGCGCAAAGCGCGGTGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGAt                       <  1:189446/151‑1 (MQ=255)

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

                                                                                                                             tAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAg                     <  1:1108648/151‑1 (MQ=255)

                                                                                                                              aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg                    >  2:3282667/1‑151 (MQ=255)

                                                                                                                               gACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGAGAGGc                   >  2:1276891/1‑151 (MQ=255)

                                                                                                                                aCATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCa                  >  1:1130671/1‑151 (MQ=255)

                                                                                                                                 cATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAc                 <  1:5216161/151‑1 (MQ=255)

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

                                                                                                                                 cATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAc                 <  1:4504985/151‑1 (MQ=255)

                                                                                                                                   tAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgc               >  1:2505192/1‑151 (MQ=255)

                                                                                                                                   tAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgc               >  1:1052336/1‑151 (MQ=255)

                                                                                                                                    aGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgcg              <  1:6249794/151‑1 (MQ=255)

                                                                                                                                    aGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgcg              <  1:6250994/151‑1 (MQ=255)

                                                                                                                                     gCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGt             <  2:3199592/151‑1 (MQ=255)

                                                                                                                                      cgtgggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTg            <  2:2330525/148‑1 (MQ=255)

                                                                                                                                      cctgggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTg            >  2:4921517/4‑151 (MQ=255)

                                                                                                                                       ctgggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTgc           <  2:5253568/149‑1 (MQ=255)

                                                                                                                                        tgggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTgcg          <  1:4946226/150‑1 (MQ=255)

                                                                                                                                          ggCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCt        <  2:5060688/151‑1 (MQ=255)

                                                                                                                                            cgcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgc      >  1:4782152/1‑151 (MQ=255)

                                                                                                                                            cgcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgc      >  1:5193702/1‑151 (MQ=38)

                                                                                                                                             gcAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgcc     >  1:4820259/1‑150 (MQ=255)

                                                                                                                                               aaaGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgcagg   >  1:5728962/1‑151 (MQ=255)

                                                                                                                                                aaGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTgcaggc  <  2:5663255/151‑1 (MQ=255)

CGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGTGCGCTGCAGGC  >  CP009974/4541500‑4541794

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