BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,211,803	G→A	33.0%	P205S (CCT→TCT)	arpA ←	ankyrin repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,211,803	0	G	A	33.0%	450.7 / inf	314	P205S (CCT→TCT)	arpA	ankyrin repeat protein
Reads supporting (aligned to +/- strand): ref base G (103/107); new base A (57/46); total (161/153)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 318 total aligned reads displayed.

TCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCG  >  CP009273/4211656‑4211950

    tATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGc                                                                                                                                                <  1:1958732/149‑1 (MQ=255)

      ttCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTg                                                                                                                                              >  2:144756/1‑149 (MQ=255)

      ttCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTg                                                                                                                                              >  2:2748042/1‑149 (MQ=255)

       tCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTgg                                                                                                                                             <  2:3202849/149‑1 (MQ=255)

         ttAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGaa                                                                                                                                           <  2:1293609/149‑1 (MQ=255)

           aGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAAc                                                                                                                                         <  2:1636632/149‑1 (MQ=255)

           aGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAAc                                                                                                                                         >  1:4134567/1‑149 (MQ=255)

              cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCaa                                                                                                                                      <  2:441632/149‑1 (MQ=255)

              cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCaa                                                                                                                                      <  1:2648281/149‑1 (MQ=255)

              cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCaa                                                                                                                                      <  1:2218029/149‑1 (MQ=255)

                 aCCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGc                                                                                                                                   <  1:1024587/149‑1 (MQ=255)

                 aCCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGc                                                                                                                                   <  1:3050566/149‑1 (MQ=255)

                 aCCGTCAAAGCATGGTTATCCACCAGTAAAATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGACAAAATATGCTTTCATAACTAGACCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGCACGGCTGGAAACCAACGc                                                                                                                                   >  2:4785868/1‑149 (MQ=255)

                     tCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaaa                                                                                                                               >  2:2937983/1‑149 (MQ=255)

                      cAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAg                                                                                                                              >  2:1236703/1‑149 (MQ=255)

                       aaaGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGt                                                                                                                             <  2:1904727/149‑1 (MQ=255)

                         aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCt                                                                                                                           >  2:2680396/1‑149 (MQ=255)

                         aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCt                                                                                                                           >  2:2527263/1‑149 (MQ=255)

                          gCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTg                                                                                                                          <  2:1088190/149‑1 (MQ=255)

                             tGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGtt                                                                                                                       >  1:1176805/1‑149 (MQ=255)

                             tGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTgggt                                                                                                                       >  1:5198533/1‑147 (MQ=255)

                               gTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtttt                                                                                                                     <  2:2772087/149‑1 (MQ=255)

                                 tATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTg                                                                                                                   >  2:5118859/1‑149 (MQ=255)

                                   tCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCa                                                                                                                 >  1:3904590/1‑149 (MQ=255)

                                   tCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCa                                                                                                                 >  2:4622554/1‑149 (MQ=255)

                                    ccaccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAt                                                                                                                <  2:2452592/149‑1 (MQ=255)

                                    ccaccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAt                                                                                                                >  1:5148827/1‑149 (MQ=255)

                                        caGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttttt                                                                                                            >  1:4024368/1‑149 (MQ=255)

                                         aGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTGTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttttt                                                                                                           <  1:3209213/149‑1 (MQ=255)

                                              tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGa                                                                                                      >  2:5029545/1‑149 (MQ=255)

                                                 ggCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGt                                                                                                   <  1:3443570/149‑1 (MQ=255)

                                                   cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGc                                                                                                 <  2:2635360/149‑1 (MQ=255)

                                                   cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGc                                                                                                 >  1:4474384/1‑149 (MQ=255)

                                                   cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGc                                                                                                 >  1:1072492/1‑149 (MQ=255)

                                                   cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGc                                                                                                 <  1:2376857/149‑1 (MQ=255)

                                                       cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga                                                                                             >  1:1850198/1‑149 (MQ=255)

                                                       cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga                                                                                             >  1:555699/1‑149 (MQ=255)

                                                         tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgaga                                                                                           >  2:4879742/1‑149 (MQ=255)

                                                           acacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGaaa                                                                                         >  1:3688027/1‑149 (MQ=255)

                                                            cacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAt                                                                                        <  1:348466/149‑1 (MQ=255)

                                                              cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                      <  2:2165682/149‑1 (MQ=255)

                                                              cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                      <  1:3254131/149‑1 (MQ=255)

                                                                gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttttt                                                                                    <  1:2341027/149‑1 (MQ=255)

                                                                gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAAATAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTTTTGAAGTGCGTGAGAAAttttt                                                                                    >  1:4655392/1‑149 (MQ=255)

                                                                      aTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTg                                                                              >  2:2673101/1‑149 (MQ=255)

                                                                        ttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAATTTTCAGCTGaa                                                                            <  2:1989373/149‑1 (MQ=255)

                                                                         tGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAag                                                                           <  2:3593108/149‑1 (MQ=255)

                                                                         tGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAag                                                                           >  1:3578225/1‑149 (MQ=255)

                                                                          gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAaga                                                                          >  1:4078329/1‑149 (MQ=255)

                                                                           aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAagag                                                                         <  2:3141288/149‑1 (MQ=255)

                                                                           aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAagag                                                                         <  1:498500/149‑1 (MQ=255)

                                                                            ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGt                                                                        >  2:429719/1‑149 (MQ=255)

                                                                                  cTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatcat                                                                  <  1:4262246/149‑1 (MQ=255)

                                                                                   ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatcatc                                                                 <  2:3730148/149‑1 (MQ=255)

                                                                                   ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatcatc                                                                 >  2:2433631/1‑149 (MQ=255)

                                                                                        aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGAt                                                            <  1:680834/149‑1 (MQ=255)

                                                                                          aaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATgg                                                          >  2:3894647/1‑149 (MQ=255)

                                                                                           atatGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATggg                                                         <  2:1860170/149‑1 (MQ=255)

                                                                                            tatGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATgggg                                                        <  2:1528182/149‑1 (MQ=255)

                                                                                             atGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGt                                                       >  1:2712492/1‑149 (MQ=255)

                                                                                             atGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGt                                                       >  2:4751715/1‑149 (MQ=255)

                                                                                               gCTTTCATAACTAGCCCGAAATCACCGTAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGttt                                                     >  2:2195290/1‑149 (MQ=255)

                                                                                                cTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTg                                                    >  2:2926401/1‑149 (MQ=255)

                                                                                                 tttCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGt                                                   <  1:4936335/149‑1 (MQ=255)

                                                                                                  ttCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTa                                                  <  1:2190693/149‑1 (MQ=255)

                                                                                                  ttCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAACTTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTa                                                  >  2:3365442/1‑149 (MQ=255)

                                                                                                   tCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTaa                                                 >  1:1365321/1‑149 (MQ=255)

                                                                                                    cATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAAt                                                >  2:898892/1‑149 (MQ=255)

                                                                                                    cATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAAt                                                <  1:3723789/149‑1 (MQ=255)

                                                                                                     aTAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                               >  2:5049119/1‑149 (MQ=255)

                                                                                                       aaCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCaa                                             <  2:4197994/149‑1 (MQ=255)

                                                                                                        aCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAAc                                            <  2:4653619/149‑1 (MQ=255)

                                                                                                            gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTc                                        <  1:4896790/149‑1 (MQ=255)

                                                                                                            gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTc                                        >  1:2515740/1‑149 (MQ=255)

                                                                                                             cccGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCACAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCt                                       <  1:2436571/149‑1 (MQ=255)

                                                                                                             cccGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCt                                       <  1:4353346/149‑1 (MQ=255)

                                                                                                              ccGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTg                                      >  2:3813192/1‑149 (MQ=255)

                                                                                                               cGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGa                                     <  1:2862628/149‑1 (MQ=255)

                                                                                                                gAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGAt                                    >  2:1339706/1‑149 (MQ=255)

                                                                                                                gAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGAt                                    <  1:4588244/149‑1 (MQ=255)

                                                                                                                 aaaTCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATa                                   <  2:4300212/149‑1 (MQ=255)

                                                                                                                   aTCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATaaa                                 <  1:5199968/149‑1 (MQ=255)

                                                                                                                    tCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAat                                <  1:1725837/149‑1 (MQ=255)

                                                                                                                    tCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAat                                >  2:4483415/1‑149 (MQ=255)

                                                                                                                    tCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAat                                >  2:3468074/1‑149 (MQ=255)

                                                                                                                     cACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAata                               <  1:4811371/149‑1 (MQ=255)

                                                                                                                           aggcagCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTa                         >  2:1423187/1‑149 (MQ=255)

                                                                                                                           aggcagCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTa                         >  1:2390991/1‑149 (MQ=255)

                                                                                                                            ggcagCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTaa                        >  1:1778906/1‑149 (MQ=255)

                                                                                                                             gcagcaACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAat                       >  2:1181803/1‑149 (MQ=255)

                                                                                                                                  aacgcaacgcCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATAtt                  <  1:2748042/149‑1 (MQ=255)

                                                                                                                                  aacgcaacgcCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATAtt                  <  1:144756/149‑1 (MQ=255)

                                                                                                                                        acgcCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTGTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcag            >  2:4033078/1‑149 (MQ=255)

                                                                                                                                            cccTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAt        >  1:266963/1‑149 (MQ=255)

                                                                                                                                            cccTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAt        >  1:1150853/1‑149 (MQ=255)

                                                                                                                                              cTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgt      <  1:127679/149‑1 (MQ=255)

                                                                                                                                              cTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgt      <  2:640835/149‑1 (MQ=255)

                                                                                                                                              cTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgt      <  1:5054842/149‑1 (MQ=255)

                                                                                                                                               tCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgtg     >  2:2334594/1‑149 (MQ=255)

                                                                                                                                                  ggACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCg  >  2:4783333/1‑149 (MQ=255)

TCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCG  >  CP009273/4211656‑4211950

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