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	97.8%	484.5 / ‑2.2	186	P66P (CCG→CCA)	RPPX_20280	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/3); new base T (93/88); total (94/92)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

GGTGCTCTGGCCGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGT  >  CP009974/4541489‑4541783

ggTGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTggg                                                                                                                                                  >  2:928532/1‑151 (MQ=21)

  tGCTCTGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGgcg                                                                                                                                                >  2:3082555/1‑151 (MQ=255)

  ggCTCTGGCTGGATAGGTAGCGCGCACGGCGGGGGTAGTCGTCCAGCCCTGGGTAGCTCAGCTCTTGGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGCACACGTAGACATAGCCTGGgcg                                                                                                                                                <  2:4954583/150‑1 (MQ=11)

   gCTCGGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGgcgc                                                                                                                                               <  2:2513611/151‑1 (MQ=16)

    ggcgggCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCa                                                                                                                                              <  1:3714798/147‑1 (MQ=255)

    ctctGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCa                                                                                                                                              <  1:2261540/151‑1 (MQ=255)

      ctGGCCGGATAGGTAGCGCTCAGGGCGGGGGTGGTCCTCCAGCAGTTGGTAGCTCAGCTGTTTGTAGAGGCCCGCGCATACGTGCTCAGGAGTGACCAGCTTTTCGACATGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCaaa                                                                                                                                            <  2:5226046/151‑1 (MQ=11)

       tGGCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAg                                                                                                                                           <  2:2839068/151‑1 (MQ=255)

         gCCGGATAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGcg                                                                                                                                         <  2:2769075/151‑2 (MQ=255)

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

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

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

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

                aGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAggg                                                                                                                                  <  1:3568881/151‑1 (MQ=19)

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

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

                  gTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGgtg                                                                                                                                >  1:1884771/1‑151 (MQ=255)

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

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

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

                       gcgcACGGGGGTGGTAGTGCTCCAGCCCTTGGTAGCGCGGCTCTTTGTAGGGGCCCGCGCGGACATGCTCATGAGTGGCCAGCTTTTCGCGCTGAGGGCGTTTCATGAACACGGAGACATAGCCCGGGCGCAAAGCACGGGGGGGGTAGGc                                                                                                                           <  2:5008277/148‑1 (MQ=11)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                  tggtAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCa                                                                                                                >  2:2596534/1‑151 (MQ=19)

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

                                    gtAGTCCTGCAGCCCTGGGTAGCTTAGCTCTGTGTAGGGGCCCGCGCCGACGTGCTCATGGATGACGGGCTTTTCGCCCTGAGGGGCTTTCTTGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGc                                                                                                              <  2:2945163/151‑1 (MQ=11)

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

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

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

                                         cctccAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTggg                                                                                                         >  2:1244713/1‑151 (MQ=21)

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

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

                                                 ccctGGTAGCGCAGCTCTTTGGAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCACTGAGGGCCGTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCtt                                                                                                 <  2:4816464/148‑1 (MQ=1)

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

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

                                                  cgtGGTAGCTCAGCTCTTGGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGTGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCCTg                                                                                                <  2:241084/149‑1 (MQ=11)

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

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

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

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

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

                                                        agctcagctcTTGGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGTCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAg                                                                                          <  1:1922602/151‑1 (MQ=255)

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

                                                         gctcagctcTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAgg                                                                                         >  1:4221402/1‑151 (MQ=9)

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

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

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

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

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

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

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

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

                                                                                                          ttttCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCggg                                        >  2:1811152/1‑151 (MQ=9)

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

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

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

                                                                                                              cGCCCTGGGGGCCTTTGATGAACACGTAGACAGAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCAc                                    <  1:2176032/151‑1 (MQ=11)

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

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

                                                                                                                cccTGAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGa                                  >  1:2239418/1‑151 (MQ=9)

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

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

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

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

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

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

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

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

                                                                                                                    gAGGGCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGc                              <  2:1287125/151‑1 (MQ=17)

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

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

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

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

                                                                                                                        gCCTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCAGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAg                          <  1:1752902/151‑1 (MQ=12)

                                                                                                                         ccTTTCATGAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAgc                         >  1:1878609/1‑151 (MQ=12)

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

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

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

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

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

                                                                                                                                 gAACACGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGc                 <  1:3107907/151‑1 (MQ=18)

                                                                                                                                   acacGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAg               >  2:4772877/1‑151 (MQ=19)

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

                                                                                                                                      cGTAGACATAGCCTGGGCGCAAAGCGCGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGAt            >  2:2349928/1‑151 (MQ=16)

                                                                                                                                       gTAGACATCGCCTGGGCGCAAAGGGCGGAGGGTGTAGGCCGAGTGCGGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATa           <  2:3115942/151‑1 (MQ=255)

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

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

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

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

GGTGCTCTGGCCGGACAGGTAGCGCTCACGGCGGTGGTAGTCCTCCAGCCCTTGGTAGCTCAGCTCTTTGTAGAGGCCCGCGCCGACGTGCTCATGAATGACCAGCTTTTCGCCCTGAGGGCCTTTCATGAACACGTAGACATAGCCCGGGCGCAAAGCACGGAGGGTGTAGGCCGAGTGCTGCAGCGTGGGGAAGCCTTGTTCCAGGTTGAAGCCTGCATCGGCATAGCGGCAGGCGGGGGCGTCGCCGGTGCGGGGCACGATGGCGTAGCGCACCGGCAGGATAGGCACGCGT  >  CP009974/4541489‑4541783

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