BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	5,595,953	T→C	L225P (CTC→CCC)	RPPX_24915 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	5,595,953	0	T	C	99.6%	609.6 / ‑5.1	250	L225P (CTC→CCC)	RPPX_24915	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (135/114); total (135/115)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 253 total aligned reads displayed.

CAGCCCCATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCTCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACAAAACCCACCC  >  CP009974/5595805‑5596099

  gCCCCATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTgg                                                                                                                                                >  1:936228/1‑151 (MQ=17)

  gCCCCATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTgg                                                                                                                                                <  1:5063989/151‑1 (MQ=17)

    cccATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGAc                                                                                                                                              >  1:3056730/1‑151 (MQ=17)

     ccATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGAcc                                                                                                                                             <  1:5386075/151‑1 (MQ=17)

      cATGTGCTGCCAGTCAGTGCCCTGCGAACCGGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGAccc                                                                                                                                            <  1:659407/151‑1 (MQ=12)

      cATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGGCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGAccc                                                                                                                                            <  2:3489329/151‑1 (MQ=17)

       aTGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCa                                                                                                                                           >  1:1215570/1‑151 (MQ=17)

        tgtgCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTCGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCaa                                                                                                                                          >  2:828028/1‑151 (MQ=14)

             tgcCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCCCTGCCCTCGCCCGCCACTACCCCTGGACCCAACCgga                                                                                                                                     >  2:4629750/1‑149 (MQ=255)

              gcCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGaaa                                                                                                                                    <  1:2526701/151‑1 (MQ=17)

               ccAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAg                                                                                                                                   <  1:2346095/151‑1 (MQ=17)

               ccAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAg                                                                                                                                   >  1:1166670/1‑151 (MQ=17)

                 agtcagtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTg                                                                                                                                 >  2:4423926/1‑151 (MQ=14)

                  gtcagtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGc                                                                                                                                >  1:1296076/1‑151 (MQ=17)

                  gtcagtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGc                                                                                                                                >  1:3402684/1‑151 (MQ=17)

                    cagtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCg                                                                                                                              <  1:1655170/151‑1 (MQ=17)

                    cagtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCg                                                                                                                              <  2:1430147/151‑1 (MQ=17)

                     agtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCgg                                                                                                                             <  2:78528/151‑1 (MQ=14)

                     agtGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCgg                                                                                                                             <  2:3552319/151‑1 (MQ=17)

                       tGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCggtg                                                                                                                           <  2:3806592/151‑1 (MQ=14)

                          ccTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCggtggtg                                                                                                                        <  1:4385626/151‑1 (MQ=17)

                             gCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGgtggtggtg                                                                                                                     >  1:475118/1‑151 (MQ=17)

                               aaaCCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGgc                                                                                                                   <  2:3477702/151‑1 (MQ=14)

                                  ccTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTCCCGGTGGTGGTGGCGCt                                                                                                                <  1:2181448/151‑1 (MQ=17)

                                   cTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTc                                                                                                               >  2:3558558/1‑151 (MQ=17)

                                   cTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTc                                                                                                               <  2:1352949/151‑1 (MQ=17)

                                     gggTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGc                                                                                                             >  2:3990245/1‑151 (MQ=14)

                                     gggTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGc                                                                                                             <  1:579947/151‑1 (MQ=17)

                                      ggTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGcc                                                                                                            <  2:4058785/151‑1 (MQ=17)

                                          gAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATg                                                                                                        >  1:2805040/1‑151 (MQ=17)

                                           aaGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGc                                                                                                       <  1:5040312/151‑1 (MQ=14)

                                            aGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGcc                                                                                                      >  1:1272502/1‑151 (MQ=17)

                                              aCTTCAAGCCACCCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGa                                                                                                    <  2:2554964/151‑1 (MQ=14)

                                               cTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATg                                                                                                        <  2:2866441/146‑1 (MQ=17)

                                               cTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATg                                                                                                        >  1:2866441/1‑146 (MQ=17)

                                                ttCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATGACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAg                                                                                                  <  2:5221951/151‑1 (MQ=14)

                                                ttCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAg                                                                                                  <  2:188714/151‑1 (MQ=14)

                                                ttCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAg                                                                                                  >  1:2885959/1‑151 (MQ=17)

                                                 tCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAgg                                                                                                 <  2:5005634/151‑1 (MQ=14)

                                                   aaGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCa                                                                                               >  2:624407/1‑151 (MQ=18)

                                                   aaGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGCGGTGGCGCTCGCCCATGCCGAAGGCa                                                                                               >  2:3907495/1‑151 (MQ=16)

                                                    aGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCAt                                                                                              <  1:4564491/151‑1 (MQ=18)

                                                      ccAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATgg                                                                                            >  1:3027314/1‑151 (MQ=21)

                                                       cAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGc                                                                                           >  2:4265994/1‑151 (MQ=25)

                                                       cAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGc                                                                                           <  2:1272502/151‑1 (MQ=21)

                                                       cAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGc                                                                                           >  1:3751027/1‑151 (MQ=21)

                                                         aCCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGccc                                                                                         >  1:1991897/1‑151 (MQ=18)

                                                          ccGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCCCCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTCGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCt                                                                                        <  1:139494/151‑1 (MQ=16)

                                                           cGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTg                                                                                       >  1:3638663/1‑151 (MQ=25)

                                                           cGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTg                                                                                       <  1:4473591/151‑1 (MQ=21)

                                                            gAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTgg                                                                                      <  2:2885959/151‑1 (MQ=21)

                                                              gcggcgCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGAc                                                                                    <  1:1617995/151‑1 (MQ=25)

                                                                ggcgCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCt                                                                                  <  1:4257271/151‑1 (MQ=21)

                                                                 gcgcATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTg                                                                                 <  1:2282936/151‑1 (MQ=21)

                                                                  cgcATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCTATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCTCTGGACCTGa                                                                                >  1:2019865/1‑151 (MQ=255)

                                                                  cgcATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGa                                                                                >  1:3802694/1‑151 (MQ=21)

                                                                       gCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTg                                                                           <  2:711489/151‑1 (MQ=21)

                                                                         cGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTc                                                                         >  1:4805388/1‑151 (MQ=21)

                                                                           cTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCgg                                                                       >  2:4918083/1‑151 (MQ=18)

                                                                             aaCCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGtc                                                                     <  1:3990245/151‑1 (MQ=18)

                                                                              aCCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGtct                                                                    <  1:1413567/151‑1 (MQ=17)

                                                                               ccTGGAGCAGCGACCCGATCACCGAAGCCTTGCCCATCGCCACCCTCACTGCCGTCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGtctc                                                                   <  2:4402360/151‑1 (MQ=11)

                                                                               ccTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCTCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGtctc                                                                   <  2:3212593/151‑1 (MQ=32)

                                                                               ccTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGtctc                                                                   <  1:4896580/151‑1 (MQ=17)

                                                                                cTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCg                                                                  >  1:4942901/1‑151 (MQ=14)

                                                                                        gcGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCa                                                          >  1:168217/1‑151 (MQ=17)

                                                                                        gcGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCa                                                          <  2:4713215/151‑1 (MQ=17)

                                                                                          gACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGc                                                        <  2:1296076/151‑1 (MQ=17)

                                                                                          gACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGc                                                        >  2:2161844/1‑151 (MQ=17)

                                                                                          gACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGc                                                        >  1:3460256/1‑151 (MQ=14)

                                                                                             ccGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATc                                                     >  2:56574/1‑151 (MQ=17)

                                                                                               gATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCaa                                                   >  2:2983292/1‑151 (MQ=17)

                                                                                               gATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCaa                                                   <  1:3883767/151‑1 (MQ=17)

                                                                                                aTCACCGAACCCTTGCCCATCGGCAACCTCCCTGCCATCGCCCGGCACTACCCCTGGACCCACCCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAAc                                                  <  2:1142930/151‑1 (MQ=255)

                                                                                                aTCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAAc                                                  <  2:1143139/151‑1 (MQ=17)

                                                                                                aTCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAAc                                                  >  1:3971011/1‑151 (MQ=17)

                                                                                                      gAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGAACTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTg                                            <  1:827809/151‑1 (MQ=17)

                                                                                                       aaaCCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGa                                           >  2:1840426/1‑151 (MQ=17)

                                                                                                        aaCCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAGAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGAc                                          <  2:1627209/151‑1 (MQ=14)

                                                                                                        aaCCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGAc                                          >  2:4295303/1‑151 (MQ=14)

                                                                                                          ccTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCt                                        <  2:2215793/151‑1 (MQ=17)

                                                                                                             tGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGAt                                     <  2:5165552/151‑1 (MQ=17)

                                                                                                                  aTCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATgccgg                                >  1:2943720/1‑151 (MQ=17)

                                                                                                                  aTCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATgccgg                                >  1:4282160/1‑151 (MQ=14)

                                                                                                                   tCGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATgccggc                               <  2:2487917/151‑1 (MQ=14)

                                                                                                                   gcGGCAGCCTCGGTGCCAGCGCCCGGCACTACCCCTGGGGGCAGCCGAAAGTGGCGGTGGTGGTGGCGCTCGCCGGTGCCGAAGGCATGGCCCTGGACCGGAGCCTGTCGGTCGCGGCCTACCAGCATCAACTGCGTGACCTGATgccggc                               <  1:1977034/150‑1 (MQ=11)

                                                                                                                    cGGCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATgccggcc                              <  2:2666728/151‑1 (MQ=17)

                                                                                                                      gCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGa                            <  1:4154495/151‑1 (MQ=17)

                                                                                                                      gCAACCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGa                            >  1:3335687/1‑151 (MQ=17)

                                                                                                                        aaCCTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGc                          <  2:2202793/151‑1 (MQ=14)

                                                                                                                          ccTCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCaa                        >  2:4182019/1‑151 (MQ=14)

                                                                                                                            tCACTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACt                      >  1:907227/1‑151 (MQ=17)

                                                                                                                              aCTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTgg                    <  2:2011513/151‑1 (MQ=17)

                                                                                                                               cTGCCATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGa                   <  1:3196449/151‑1 (MQ=17)

                                                                                                                                   cATCGCCCGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAacac               <  1:5295232/151‑1 (MQ=21)

                                                                                                                                         ccGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACaaaa         >  1:957315/1‑151 (MQ=25)

                                                                                                                                          cGGCACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACAAAac        <  1:5315432/151‑1 (MQ=25)

                                                                                                                                             cACTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACAAAaccca     >  1:4660442/1‑151 (MQ=25)

                                                                                                                                               cTACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACAAAacccacc   >  1:192567/1‑151 (MQ=37)

                                                                                                                                                tACCCCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACAAAacccaccc  >  2:2794870/1‑151 (MQ=38)

CAGCCCCATGTGCTGCCAGTCAGTGCCCTGCAAACCTGGGTCGAAGACTTCAAGCCAACCGAGCGGCGCATGCCGCTAACCTGGAGCAGCGACCCGATCACCGAAACCTTGCCCATCGGCAACCTCACTGCCATCGCCCGGCACTACCTCTGGACCCAACCGAAAGTGCCGGTGGTGGTGGCGCTCGCCGATGCCGAAGGCATGGCCCTGGACCTGAGCCTGTCGGTCTCGGCCTACCAGCATCAACTGCGTGACCTGATGCCGGCCGAGCAACTGGAACACACAAAACCCACCC  >  CP009974/5595805‑5596099

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