BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,282,608	T→C	5.7%	noncoding (85/171 nt) T9A (ACT→GCT)	rttR ← tpr ←	rtT sRNA, processed from tyrT transcript protamine‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,282,608	0	T	C	5.7%	344.1 / 6.0	120	noncoding (85/171 nt) T9A (ACT→GCT)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand): ref base T (46/64); new base C (2/5); total (51/69)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTG  >  CP009273/1282461‑1282756

ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGCAGTCGTCCACTCTCTCTCGCTGGCGCGCGACTc                                                                                                                                                     >  2:663001/1‑149 (MQ=14)

ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTAGTTGCCTGCGGCAACGCTCTCTAGCTGGCGCTCGAGta                                                                                                                                                     >  2:661170/1‑148 (MQ=255)

    aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGAGGCAAACCACTCTCGCTGGCGCGCGAGTCCaaa                                                                                                                                                 >  2:2300379/1‑148 (MQ=14)

     tGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGTCCGGCGGCAAACCCCTCTCGCGGGCGCGCGAGTCGAAcc                                                                                                                                                >  1:4295541/1‑149 (MQ=14)

        tcAGAATATTTGATTGCGGATTCGATTGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                                             <  1:429305/148‑1 (MQ=255)

        gCAGAAGTTTTGATTGCGGATTCGATTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                                             <  1:429280/149‑1 (MQ=255)

         cAGAATTTTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTgg                                                                                                                                            <  2:499507/149‑1 (MQ=255)

             atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGAGGCAAACCACAAACGCTAGCGCTCCAGGCGAAACATGGTCGa                                                                                                                                        >  2:3943743/1‑149 (MQ=14)

                tttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCACACTCGCTTGCGCTAGAGTCGAAACTTGGTAGaata                                                                                                                                     >  2:2063257/1‑147 (MQ=14)

                 ttGATTGCGGATTCGATTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCt                                                                                                                                    <  1:1357064/149‑1 (MQ=255)

                   gATTGCGGATTCGTTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                  <  2:3578688/149‑1 (MQ=255)

                   gATTGCGGATTCGGTTGAGAGTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                  <  2:3578614/149‑1 (MQ=255)

                   gATTGCGGATTCGCTTGTGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                  <  1:2778329/149‑1 (MQ=255)

                   gATTGCGGATTCGATTGAGAGTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                                                  <  2:3578690/149‑1 (MQ=255)

                     ttGCGGATTCGATTGAGATTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtctc                                                                                                                                <  1:1659774/149‑1 (MQ=255)

                          gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCGGAGGCCAACCACTCTCGCTGGCGCGAGAGGGGGAACGTGGTCGAAGCTTATCAAcca                                                                                                                           >  2:1493579/1‑148 (MQ=14)

                            ttCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAACCACACTCGCGGGCCCTCCCGTGGAACCATGGTGGAAGCATCTCaaacata                                                                                                                         >  2:2232423/1‑143 (MQ=14)

                             tCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTcc                                                                                                                        <  1:3094880/149‑1 (MQ=255)

                              cgattGAGATTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                                       <  2:3008641/146‑1 (MQ=255)

                              cgattGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                                       <  2:3008642/146‑1 (MQ=255)

                                 ttGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGGGGCAAAGGTCTCAAGCTTGCGCGAGAGTCGAACCTATGGCGAATCTTATAATACTTCCACGc                                                                                                                    >  2:3998458/1‑149 (MQ=14)

                                     gagTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                                <  2:4935553/149‑1 (MQ=255)

                                     gagTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                                                <  1:1125129/149‑1 (MQ=255)

                                       gTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGc                                                                                                              <  2:3158821/149‑1 (MQ=255)

                                         tCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCGCCCTCGCGGGCCCTCGAGTCCCACCGTGGGAGAAGCTAAACATCCTACCCCGCAACGGCaa                                                                                                            >  2:2881539/1‑148 (MQ=14)

                                            gggTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat                                                                                                         <  1:2085771/149‑1 (MQ=37)

                                            gggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTATCTCGCTTGCGCTCGAGTCGAACCATGGTCGAAGCATCTCAAACATCCCCCCAAGGGCCGaaa                                                                                                         >  2:4279867/1‑148 (MQ=14)

                                                        aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTAGCTGGCGCACGAGTCGAACCATGGTCGAAGCTTCTCAACCTTCCCCCCAAGGGCAGAATATTTAATTGCgg                                                                                             >  2:797746/1‑149 (MQ=14)

                                                        aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCTCACTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCACCCCAAGGGCAGAATATTTAATTGCgg                                                                                             >  2:797749/1‑149 (MQ=14)

                                                         gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGGCGCGCGAGTCGAACCTAGGTCGAAGCATAACATCCATACCCGCCAGGGCAGAATATTTAAtttccga                                                                                            >  2:5154729/1‑144 (MQ=11)

                                                         gTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCACTCACGCAGGCGCAAGAGTCGAACCCAGGTCCAATCTTCTCAACCTTACCCGCAAGGGCAGAAAAATTAATAGCggc                                                                                            >  2:3205761/1‑148 (MQ=11)

                                                                      gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCTCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGt                                                                               <  1:1014781/149‑1 (MQ=255)

                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTGGAAGCGTCTCATCCCTCCCCGCATGGGCAGAAGATTTAATTTCGGATTAGTTGGGAAGtt                                                                              >  2:2135431/1‑149 (MQ=14)

                                                                          aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTAATTGCGGGTTAGTTGGGGAGTTCAg                                                                           >  2:4235157/1‑149 (MQ=255)

                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                         <  1:4991883/149‑1 (MQ=255)

                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                         <  1:4991888/149‑1 (MQ=255)

                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAggg                                                                         <  2:3610673/149‑1 (MQ=255)

                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGa                                                                        <  1:3024697/149‑1 (MQ=255)

                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGa                                                                        <  1:3024741/149‑1 (MQ=255)

                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGa                                                                        <  1:4684084/149‑1 (MQ=255)

                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAc                                                                       <  1:4195623/149‑1 (MQ=255)

                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                   <  1:3895176/121‑1 (MQ=37)

                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                   >  2:3895176/1‑121 (MQ=37)

                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                                                      <  1:843959/149‑1 (MQ=255)

                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                                                      <  1:4715135/149‑1 (MQ=255)

                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                                                      <  1:1056089/149‑1 (MQ=255)

                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGa                                                                 <  2:1922800/149‑1 (MQ=255)

                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTGGGCAGCGGCCACCCTCAATCGCTGGGGCGCGAGTCGAACCGGGGTCGAAGCGTCACAACCGTGCCCGCCGGGGGAGAAAAGTTAATTTCGGCTTGGTTGGGGAGTTCAGGGCCGTTTgg                                                                 >  2:3102983/1‑148 (MQ=9)

                                                                                          gcgcgcTTCGCCCTTCGGGTCGGTGCCTGCGGCTACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCTAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGa                                                           <  1:3938057/149‑1 (MQ=255)

                                                                                              gcTTCGCCCTTCGGGTCGTTGCCTGCGGCCACGCTCTCTCGCTGGCGCTCGAGGCGTACCTTTGTCGAATCTTCGAATACATCCCCCCCAGGGGCGAATATTTAATTGGGGGTTAGGTGGGGAGTTAAGGGACGTTTTTAAGTGAttgt                                                       >  2:4285182/1‑149 (MQ=2)

                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtg                                                   <  2:2392611/149‑1 (MQ=255)

                                                                                                           tgTCGTTGCCTGCGGCACCGCTCACTCGCTGACGCTCTAGCCGAACCTAGGTCGAAGAGTCTCATCCGACCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGTAAGTTCAGGGCCTTTTGAAAGTGATGGTGGTGGGGGAAGGa                                          <  2:2409006/148‑1 (MQ=11)

                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAt                                         <  1:2319095/149‑1 (MQ=255)

                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGACCCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                        <  2:4322179/149‑1 (MQ=255)

                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt                                        <  2:326996/149‑1 (MQ=255)

                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                               <  1:2028377/74‑1 (MQ=39)

                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                               >  2:2028377/1‑74 (MQ=39)

                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc                                   >  2:4519569/1‑149 (MQ=255)

                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc                                   >  2:4521873/1‑149 (MQ=255)

                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAc                                   >  2:3854209/1‑149 (MQ=255)

                                                                                                                  gCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTGTCAACCATCCCCGCATGGGCAGAAAATTTAATTTCGGATTCGTTGGGAAGTTCAGGCACTTTTGAAATTGTTGGGGGTGGGGGGCGGATTAGaag                                   >  1:3577634/1‑148 (MQ=17)

                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTTGGAAGTTCAGGGACTTTTTAAAGTGATGGTGGTGGGGGAAGGATTAGAAcca                                 >  1:2316138/1‑148 (MQ=255)

                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATACTTCCCCGCAGGGGCAGAATATTTAATTGCGGATTAGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGGAGGATTCGaaaca                                 >  1:1047448/1‑146 (MQ=17)

                                                                                                                     tGCGGCTTCGCTCTCTCGCGGTCGCTCGAGTCGAACCGCGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt                                <  1:4481910/149‑1 (MQ=17)

                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTGGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCtt                                <  1:4481899/149‑1 (MQ=255)

                                                                                                                      gCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc                               <  1:530640/149‑1 (MQ=255)

                                                                                                                      gCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTc                               <  1:2306899/149‑1 (MQ=255)

                                                                                                                       cGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCg                              <  1:4025645/149‑1 (MQ=255)

                                                                                                                         gCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTATCATCCTTACCCGCAAGGGCAGAATATTTAATTGCGGATTGGTTGGGTAGGTCATGGAGTTTTGTAATGGATGGTGGGGGGGGAAGGATTAGCACCTTCGaa                            >  2:1732922/1‑149 (MQ=17)

                                                                                                                           aaCGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGt                                                                                       >  1:2584329/1‑88 (MQ=18)

                                                                                                                            aCGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTccc                                                                                                                       <  1:4402256/55‑1 (MQ=35)

                                                                                                                              gCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAATCTTCTCAAACTTCCCCGCATGGGCAGAATATTTAAATGCGGATGCGTTGGGACATTCAGGGCCGTTTGATAGTGATGGTGGGGGGGGTAGGATTCGAAACGTCGAAGACGa                       >  2:4830906/1‑149 (MQ=17)

                                                                                                                              gCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGa                       >  2:2869301/1‑149 (MQ=255)

                                                                                                                              gCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGa                       <  2:2333232/149‑1 (MQ=255)

                                                                                                                              gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                             >  1:2581772/1‑127 (MQ=255)

                                                                                                                              gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                                             <  2:2581772/127‑1 (MQ=255)

                                                                                                                               ctctctCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGAACCTTAGAAGTCGAt                      >  1:2973865/1‑149 (MQ=255)

                                                                                                                               ctctctCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAAACTAACCCGCCTGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTAAGGGACGTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTAAGAAGTCGAt                      >  1:958245/1‑149 (MQ=17)

                                                                                                                                tctctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATCCCCGCAGGGGCAGAATATTTAATTGGGGATTGGTTTGGGAGTTGAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGTTTCGAACCTTCGAAGTCGCTg                     >  1:1769914/1‑149 (MQ=17)

                                                                                                                                 ctctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                    <  1:3198880/149‑1 (MQ=255)

                                                                                                                                 ctctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGa                    <  1:3745858/149‑1 (MQ=255)

                                                                                                                                 ctctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAATATTTAATTGCGGAGTCGTGGGGAAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGCTTATATCCTTAGAAGACGaagt                    >  1:4302549/1‑146 (MQ=17)

                                                                                                                                  tctcGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGAc                   <  1:2159766/149‑1 (MQ=255)

                                                                                                                                    tcGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCa                                                                            >  1:1714477/1‑90 (MQ=255)

                                                                                                                                     cGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGc                >  2:1573972/1‑149 (MQ=255)

                                                                                                                                       cTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAg              >  1:4697299/1‑149 (MQ=255)

                                                                                                                                       cTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGGTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAg              >  2:4911853/1‑149 (MQ=255)

                                                                                                                                         ggagcTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGAt            <  1:4554150/146‑1 (MQ=255)

                                                                                                                                          gcgcTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGAtt           <  1:3406237/149‑1 (MQ=255)

                                                                                                                                          ccgcTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGAtt           <  1:608743/148‑1 (MQ=255)

                                                                                                                                            gcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                   >  1:4834455/1‑59 (MQ=16)

                                                                                                                                            gcTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTaa                                                                                                   <  2:4834455/59‑1 (MQ=16)

                                                                                                                                             cTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAc        <  2:2299049/149‑1 (MQ=255)

                                                                                                                                              tCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTAGAAGTCGATGACGGCAGATTTaaa       >  1:3382216/1‑147 (MQ=255)

                                                                                                                                                gAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGt     <  1:1298217/149‑1 (MQ=255)

                                                                                                                                                gAGTCGAACCTTGGTCGAAGCTTCACATCCGTCCCCGCATGGGCAGAAAATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGTTAGTGATGGTGGTGGGGGGTGGATTCGAACCTTCGAAGTCGATGTCGGCGGATTTATAGt     >  1:2634010/1‑149 (MQ=21)

                                                                                                                                                  gTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCCCCTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAAACATCGAAGTCGATGTCGGCAGATTTAAAGTCt   >  1:2388837/1‑149 (MQ=18)

                                                                                                                                                   tCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg  <  2:313418/149‑1 (MQ=255)

                                                                                                                                                   tCGAACCTTGGTCGAAGCTTCACATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg  <  1:2247872/149‑1 (MQ=255)

                                                                                                                                                   tCGAACCTTGCACGAAGCTTCTCATCCCTCCCCGCATGGTCATAATATTTAATTGCGGATTCGTTTGTAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTg  <  2:1330254/149‑1 (MQ=21)

CCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTG  >  CP009273/1282461‑1282756

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