BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009974	4,541,636	C→T	P66P (CCG→CCA)	RPPX_20280 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009974	4,541,636	0	C	T	99.2%	685.7 / ‑5.5	250	P66P (CCG→CCA)	RPPX_20280	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (130/118); total (131/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

GGTGCTCTGGCCGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCG  >  CP009974/4541489‑4541782

 gTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGgc                                                                                                                                                <  2:3251931/151‑1 (MQ=21)

  tGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGgcg                                                                                                                                               <  1:3602945/151‑1 (MQ=255)

  tGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGgcg                                                                                                                                               >  2:5383080/1‑151 (MQ=25)

   gCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGgcgc                                                                                                                                              >  1:4560189/1‑151 (MQ=25)

     tctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaa                                                                                                                                            >  1:5114081/1‑151 (MQ=25)

     tctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaa                                                                                                                                            <  2:3992552/151‑1 (MQ=21)

       tGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAg                                                                                                                                          >  2:6624111/1‑151 (MQ=255)

       tGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAg                                                                                                                                          <  1:2227320/151‑1 (MQ=21)

        ggCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGc                                                                                                                                         >  2:1123107/1‑151 (MQ=255)

           gggATAGGTAGGGCTCACGGCGGTGGTGGGCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGGCGACGTGCTCATGAATGACCAGCGTGGCGGCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCCAAGcgcg                                                                                                                                      <  2:5364458/150‑1 (MQ=11)

            ggataGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCgg                                                                                                                                     <  1:2321055/147‑1 (MQ=255)

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

                aGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAggg                                                                                                                                 >  1:4340487/1‑151 (MQ=255)

                aGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAggg                                                                                                                                 >  2:4509013/1‑151 (MQ=19)

                 ggTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgt                                                                                                                                <  2:1201115/151‑1 (MQ=21)

                    ggcgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTa                                                                                                                             <  1:3503205/150‑1 (MQ=255)

                    aGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGCGGGTGTa                                                                                                                             >  1:542027/1‑151 (MQ=255)

                    aGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTa                                                                                                                             >  2:6231866/1‑151 (MQ=25)

                      cgcTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAgg                                                                                                                           >  2:6101895/1‑151 (MQ=25)

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

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

                        cTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGcc                                                                                                                         <  1:1076835/151‑1 (MQ=19)

                         tCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCg                                                                                                                        >  2:6398736/1‑151 (MQ=25)

                         tCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCg                                                                                                                        <  1:1228343/151‑1 (MQ=19)

                            cggcggTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGt                                                                                                                     >  2:21639/1‑151 (MQ=255)

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

                                ggtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctg                                                                                                                 <  1:4430814/151‑1 (MQ=19)

                                ggtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctg                                                                                                                 <  1:5196898/151‑1 (MQ=255)

                                 gtggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGtgctgc                                                                                                                >  2:6639328/1‑151 (MQ=255)

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

                                    gtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTGTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGc                                                                                                             <  2:4824675/151‑1 (MQ=255)

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

                                      aGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGt                                                                                                           >  1:1401271/1‑151 (MQ=255)

                                       gTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTGTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTg                                                                                                          >  1:3424807/1‑151 (MQ=16)

                                        tcctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTgg                                                                                                         <  1:5946009/151‑1 (MQ=255)

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

                                             cAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAg                                                                                                    >  2:5002087/1‑151 (MQ=19)

                                              aGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcgcg                                                                                                                                      >  2:2455925/1‑113 (MQ=255)

                                              aGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGc                                                                                                   >  2:7087291/1‑151 (MQ=25)

                                               gCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGcc                                                                                                  >  2:1845583/1‑151 (MQ=255)

                                                 ccTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                <  1:2056501/151‑1 (MQ=255)

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

                                                  cTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAACGCGCGGAGGGTGTAGGCCGAGCGCTGCAGCGTGGGGCAGCCtcg                                                                                               >  2:6478711/1‑149 (MQ=11)

                                                  cTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTg                                                                                               <  1:581856/151‑1 (MQ=19)

                                                   ttGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGt                                                                                              >  1:6045965/1‑151 (MQ=19)

                                                     ggTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTc                                                                                            <  1:6492942/151‑1 (MQ=16)

                                                      gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                           >  1:1736777/1‑151 (MQ=16)

                                                      gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                           >  2:5217881/1‑151 (MQ=16)

                                                      gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                           <  1:278590/151‑1 (MQ=17)

                                                      gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                           <  2:2773841/151‑1 (MQ=16)

                                                      gTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTcc                                                                                           >  1:939804/1‑151 (MQ=18)

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

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

                                                                                                       aGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGc                                          >  1:3773560/1‑151 (MQ=9)

                                                                                                        gCTTTTGGCCCTGAGGGCCTGTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCg                                         <  2:5748783/151‑1 (MQ=9)

                                                                                                            ttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGc                                     >  1:2219495/1‑151 (MQ=14)

                                                                                                              cGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCAc                                   <  1:4034275/151‑1 (MQ=12)

                                                                                                               gCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACg                                  <  2:2667158/151‑1 (MQ=12)

                                                                                                                cccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCCCCGCGCGGAGGGTGCAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTCCATCGGCATAGCGGCCGGCGGGCGCCTCGCCGGTCCGGGGCCCgc                                 >  2:4239787/1‑150 (MQ=11)

                                                                                                                cccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACgg                                 >  2:5499258/1‑150 (MQ=9)

                                                                                                                 ccTGAGGGCCTTTCCTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGAt                                <  2:712250/151‑1 (MQ=12)

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

                                                                                                                   tGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGTAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATgg                              <  1:3148867/151‑1 (MQ=9)

                                                                                                                    gAGGGCCTTTGATGAACACGTAGGCAGGGCCTGGGCGCAAGGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCTGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                             <  1:5479489/151‑1 (MQ=11)

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

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

                                                                                                                      gggCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGt                           <  1:6992548/151‑1 (MQ=17)

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

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

                                                                                                                        gCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCCCGATGGCGTAg                         >  2:4643250/1‑151 (MQ=9)

                                                                                                                         ccTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgc                        <  1:2425582/151‑1 (MQ=12)

                                                                                                                          cTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgcg                       >  1:4927064/1‑151 (MQ=9)

                                                                                                                           tttCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgcgc                      >  2:4764303/1‑151 (MQ=9)

                                                                                                                             tCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAc                    >  1:5281904/1‑151 (MQ=9)

                                                                                                                             tCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAc                    >  1:4493123/1‑151 (MQ=14)

                                                                                                                             tCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAc                    <  1:169536/151‑1 (MQ=12)

                                                                                                                              cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAcc                   >  2:1260943/1‑151 (MQ=12)

                                                                                                                              cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAcc                   >  2:405818/1‑151 (MQ=16)

                                                                                                                              cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCAcc                   <  1:6653680/151‑1 (MQ=16)

                                                                                                                              cATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGGGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTCGCGCAcc                   >  2:5818834/1‑151 (MQ=255)

                                                                                                                               aTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCg                  <  2:542027/151‑1 (MQ=18)

                                                                                                                                tGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGTGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCgg                 <  2:5248868/151‑1 (MQ=255)

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

                                                                                                                                  aaCACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGCGTGCTGCAGCGTGGCGAAGCATTGCACCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCa               >  2:5109190/1‑151 (MQ=11)

                                                                                                                                     acGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCTGTGCGGGGCACGATGGCGTAGCGCCCCGGCAGGa            >  1:2492713/1‑151 (MQ=255)

                                                                                                                                     acGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGa            >  1:4392103/1‑151 (MQ=16)

                                                                                                                                       gTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATa          >  1:6426958/1‑151 (MQ=255)

                                                                                                                                         aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATacg        >  2:2298904/1‑149 (MQ=255)

                                                                                                                                         aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg        <  1:717321/151‑1 (MQ=255)

                                                                                                                                         aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg        <  2:3347786/151‑1 (MQ=255)

                                                                                                                                         aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg        <  1:981171/151‑1 (MQ=255)

                                                                                                                                         aGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGCGCTGCAGCGTGGGGAAGCCTGGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAgg        <  1:6717209/151‑1 (MQ=255)

                                                                                                                                          gACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGc       >  1:776431/1‑151 (MQ=255)

                                                                                                                                           aCATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCa      <  1:4884838/151‑1 (MQ=255)

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

                                                                                                                                           aCATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCa      <  1:5030768/151‑1 (MQ=255)

                                                                                                                                            cATAGCCTGGGCGCCAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAc     <  1:3649482/151‑1 (MQ=255)

                                                                                                                                              tAGCCTGTGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGGAGCGTGGGGAAGCCTTGCTCCCGGTTGAGGCCTGCATCGCCATAGCGGCAGGCGGGGGCGACGCCGGTGGGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgc   <  2:6751855/151‑1 (MQ=11)

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

                                                                                                                                              tAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCAcgc   <  1:1011391/151‑1 (MQ=255)

GGTGCTCTGGCCGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCG  >  CP009974/4541489‑4541782

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