BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,858,963	A→T	12.2%	V224E (GTA→GAA)	yidE ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,858,963	0	A	T	12.2%	822.2 / 89.7	279	V224E (GTA→GAA)	yidE	putative transporter
Reads supporting (aligned to +/- strand): ref base A (113/131); new base T (15/19); total (129/150)
Fisher's exact test for biased strand distribution p-value = 8.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 285 total aligned reads displayed.

CACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAT  >  CP009273/3858815‑3859111

cACCAGTTGCGCATTATGTACCTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTa                                                                                                                                                      <  1:2337161/149‑1 (MQ=255)

  ccAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACa                                                                                                                                                    >  2:249893/1‑149 (MQ=255)

  ccAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACa                                                                                                                                                    >  2:2770587/1‑149 (MQ=255)

     gTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCt                                                                                                                                                 <  1:3381866/149‑1 (MQ=255)

     gTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCt                                                                                                                                                 >  1:1342765/1‑149 (MQ=255)

     gTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCt                                                                                                                                                 >  1:2597155/1‑149 (MQ=255)

       tGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                                                                                                               <  1:1935107/149‑1 (MQ=255)

           cATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATg                                                                                                                                           <  2:4598015/149‑1 (MQ=255)

             ttATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGc                                                                                                                                         >  1:4164917/1‑149 (MQ=255)

              tATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCt                                                                                                                                        <  1:2013252/149‑1 (MQ=255)

                tGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaa                                                                                                                                      <  1:1137863/149‑1 (MQ=255)

                tGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaa                                                                                                                                      >  2:1094752/1‑149 (MQ=255)

                tGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaa                                                                                                                                      <  1:4671339/149‑1 (MQ=255)

                  tAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAAt                                                                                                                                    <  2:3189135/149‑1 (MQ=255)

                    aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCa                                                                                                                                  >  1:4274940/1‑149 (MQ=255)

                    aaTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGAAGATAATTTTGTCTCCGTTGAGAATCGGTACATCTGTAATTGCTAAATCa                                                                                                                                  >  1:3211759/1‑149 (MQ=255)

                     aTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                                                                                                                 <  2:1227753/149‑1 (MQ=255)

                          ctggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAgg                                                                                                                            <  1:4130730/149‑1 (MQ=255)

                           tggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGt                                                                                                                           >  2:881605/1‑149 (MQ=255)

                            ggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGtt                                                                                                                          >  2:3071751/1‑149 (MQ=255)

                            ggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGtt                                                                                                                          <  2:4138362/149‑1 (MQ=255)

                              ctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAg                                                                                                                        <  1:1587763/149‑1 (MQ=255)

                                gACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAggg                                                                                                                      >  1:1185194/1‑149 (MQ=255)

                                 aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGt                                                                                                                     <  2:1171327/149‑1 (MQ=255)

                                   ccaccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtc                                                                                                                   <  2:3534358/149‑1 (MQ=255)

                                   ccaccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtc                                                                                                                   <  2:3434165/149‑1 (MQ=255)

                                    caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                                                                  <  2:2221313/149‑1 (MQ=255)

                                    caccaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtct                                                                                                                  <  1:4361325/149‑1 (MQ=255)

                                     accaGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                                                                                                                 <  1:4744481/149‑1 (MQ=255)

                                        aGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                              <  2:2012214/149‑1 (MQ=255)

                                        aGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAc                                                                                                              <  1:74194/149‑1 (MQ=255)

                                         gATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCATTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcg                                                                                                             >  1:4553145/1‑149 (MQ=255)

                                               gcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                       <  2:743675/149‑1 (MQ=255)

                                               gcaAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAata                                                                                                       <  2:1439603/149‑1 (MQ=255)

                                                 aaaTCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATAtt                                                                                                     >  1:2101752/1‑149 (MQ=255)

                                                      gCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAg                                                                                                >  1:2875205/1‑149 (MQ=255)

                                                         cAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCt                                                                                             <  1:2183754/149‑1 (MQ=255)

                                                         cAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCt                                                                                             >  1:3944555/1‑149 (MQ=255)

                                                          cgtgggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt                                                                                            <  2:1739889/145‑1 (MQ=255)

                                                          aGTTGGATATTGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt                                                                                            <  1:2238440/149‑1 (MQ=255)

                                                              ggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATc                                                                                        <  1:3338439/149‑1 (MQ=255)

                                                              ggATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATc                                                                                        <  2:3303334/149‑1 (MQ=255)

                                                                     ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcgcgc                                                                                 <  1:2405725/149‑1 (MQ=255)

                                                                     ggTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAgcgcgc                                                                                 >  2:4654788/1‑149 (MQ=255)

                                                                      gTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                <  2:2251672/149‑1 (MQ=255)

                                                                      gTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgcc                                                                                <  2:3856614/149‑1 (MQ=255)

                                                                        aTCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGcgccgc                                                                              >  1:48322/1‑149 (MQ=255)

                                                                         tCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCgccgcc                                                                             <  2:4712302/149‑1 (MQ=255)

                                                                            ggCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCTGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCAtg                                                                          >  2:3792522/1‑148 (MQ=255)

                                                                               gAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGt                                                                       >  2:3877316/1‑149 (MQ=255)

                                                                               gAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGt                                                                       >  2:1333703/1‑149 (MQ=255)

                                                                                aaGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTg                                                                      <  1:3795889/149‑1 (MQ=255)

                                                                                aaGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTg                                                                      <  1:954337/149‑1 (MQ=255)

                                                                                aaGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTg                                                                      >  2:190903/1‑149 (MQ=255)

                                                                                 aGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGc                                                                     >  1:630160/1‑149 (MQ=255)

                                                                                   gAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGt                                                                   >  2:155924/1‑149 (MQ=255)

                                                                                       ttttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAg                                                               >  1:3312478/1‑149 (MQ=255)

                                                                                       ttttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAg                                                               >  2:469587/1‑149 (MQ=255)

                                                                                        tttAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGa                                                              >  1:2704874/1‑149 (MQ=255)

                                                                                             ggTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGt                                                         <  2:2529720/149‑1 (MQ=255)

                                                                                             ggTTTCGTCGCGTTTCAGACGCGCGCAGATCAGTTTGTCGCCGTTGAGAATCGGGACATCTGTCATGGCGAACTCATGCAGGTTAGGGTTCTCAACGCTAATATTGAGAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCgg                                                         <  2:1413559/149‑2 (MQ=255)

                                                                                                ttcttcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtgct                                                      <  1:3757038/149‑1 (MQ=255)

                                                                                                ttcgtcGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAACTCATGCAGGTTAGGGTTCTCAACGCGAATATTGAGAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGtgct                                                      <  1:3756778/149‑1 (MQ=255)

                                                                                                     cgcgTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAGGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGa                                                 <  2:2561166/149‑1 (MQ=255)

                                                                                                        gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                              >  2:2718329/1‑149 (MQ=255)

                                                                                                        gTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCt                                              >  2:4518151/1‑149 (MQ=255)

                                                                                                            cAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTg                                          <  1:3511233/149‑1 (MQ=255)

                                                                                                             aGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGt                                         >  1:835942/1‑149 (MQ=255)

                                                                                                             aGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGt                                         >  1:4769010/1‑149 (MQ=255)

                                                                                                                    aGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACa                                  <  2:4228157/149‑1 (MQ=255)

                                                                                                                       tgATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGAGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                               <  2:415276/148‑1 (MQ=255)

                                                                                                                       aGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                               <  1:4693724/149‑1 (MQ=255)

                                                                                                                       aGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTg                               >  1:1016827/1‑149 (MQ=255)

                                                                                                                        gATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGa                              >  2:2124878/1‑149 (MQ=255)

                                                                                                                        gATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGa                              <  1:2983376/149‑1 (MQ=255)

                                                                                                                        gATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGa                              >  2:3712624/1‑149 (MQ=255)

                                                                                                                            aTTTTGTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAcgcg                          >  2:881130/1‑149 (MQ=255)

                                                                                                                            aTTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAcgcg                          >  1:875406/1‑149 (MQ=255)

                                                                                                                            aTTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGAcgcg                          >  2:1373898/1‑149 (MQ=255)

                                                                                                                             ttttGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCgg                         >  2:3918852/1‑149 (MQ=255)

                                                                                                                             ttttGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCgg                         >  2:173187/1‑149 (MQ=255)

                                                                                                                                 gTCGCCGTTGAGAATCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGaaaa                     <  1:1102545/149‑1 (MQ=255)

                                                                                                                                 gTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGaaaa                     <  1:1510007/149‑1 (MQ=255)

                                                                                                                                  tCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAAt                    >  2:154645/1‑149 (MQ=255)

                                                                                                                                  tCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAAt                    <  2:4372157/149‑1 (MQ=255)

                                                                                                                                     ccGATGAGAATCGGTACATCTTTAATGGCTAACTCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAAc                 <  2:1503042/149‑1 (MQ=255)

                                                                                                                                        ttGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCg              >  2:1316636/1‑149 (MQ=255)

                                                                                                                                        ttGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACAGTGACGCGGAAAATAACCCg              >  2:3639125/1‑149 (MQ=255)

                                                                                                                                         tGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGc             <  2:4428602/149‑1 (MQ=255)

                                                                                                                                           agaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCaa           <  1:4627838/149‑1 (MQ=255)

                                                                                                                                            gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGACGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAAc          >  2:2993720/1‑149 (MQ=255)

                                                                                                                                            gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAAc          >  1:1765215/1‑149 (MQ=255)

                                                                                                                                            gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAAc          >  2:4001839/1‑149 (MQ=255)

                                                                                                                                            gaATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAAc          >  1:3048970/1‑149 (MQ=255)

                                                                                                                                             aaTCGGTTCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa         >  1:548427/1‑149 (MQ=255)

                                                                                                                                             aaTCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa         >  1:2922787/1‑149 (MQ=255)

                                                                                                                                             aaTCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACa         <  2:513412/149‑1 (MQ=255)

                                                                                                                                                    aCATCTTTAATGGCTAAATCATGTAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAt  <  1:4344042/149‑1 (MQ=255)

                                                                                                                                                    aCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAt  <  2:1831287/149‑1 (MQ=255)

                                                                                                                                                    aCATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAt  <  1:4484966/149‑1 (MQ=255)

CACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTTGATCAGCGCGCCGCCATTGGTGCGTGAAGACTCGTGCTGCTGAGCTTCTGTCTCGACATTGACGCGGAAAATAACCCGCAACATCCACAT  >  CP009273/3858815‑3859111

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