BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,211,803	G→A	86.4%	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	86.4%	692.3 / 87.6	243	P205S (CCT→TCT)	arpA	ankyrin repeat protein
Reads supporting (aligned to +/- strand): ref base G (15/18); new base A (116/94); total (131/112)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 243 total aligned reads displayed.

CTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTC  >  CP009273/4211657‑4211949

cTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACg                                                                                                                                                 <  2:4123509/150‑1 (MQ=255)

cTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACg                                                                                                                                                 <  1:4108550/150‑1 (MQ=255)

cTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACg                                                                                                                                                 <  2:3209522/150‑1 (MQ=255)

 tttATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACgg                                                                                                                                                >  1:472287/1‑150 (MQ=255)

    aTTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTg                                                                                                                                             <  1:516800/150‑1 (MQ=255)

       cTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGaa                                                                                                                                          >  1:1629321/1‑150 (MQ=255)

         tAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAAc                                                                                                                                        >  2:570991/1‑150 (MQ=255)

         tAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAAc                                                                                                                                        <  1:3881299/150‑1 (MQ=255)

             cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAAc                                                                                                                                    <  2:1024819/150‑1 (MQ=255)

              ggACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACg                                                                                                                                   >  1:2510978/1‑150 (MQ=255)

               gACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGc                                                                                                                                  >  2:530539/1‑150 (MQ=255)

                aCCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGcc                                                                                                                                 >  1:378092/1‑150 (MQ=255)

                  cGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCaa                                                                                                                               <  2:4490946/150‑1 (MQ=255)

                    tCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAg                                                                                                                             >  1:2186154/1‑150 (MQ=255)

                     cAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGt                                                                                                                            >  2:71580/1‑150 (MQ=255)

                      aaaGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTc                                                                                                                           >  1:3253688/1‑150 (MQ=255)

                       aaGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCt                                                                                                                          >  2:2573402/1‑150 (MQ=255)

                        aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTg                                                                                                                         >  2:2858603/1‑150 (MQ=255)

                        aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTg                                                                                                                         <  1:1701057/150‑1 (MQ=255)

                         gCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTgg                                                                                                                        <  1:2966368/150‑1 (MQ=255)

                           aTGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtt                                                                                                                      <  2:1597243/150‑1 (MQ=255)

                           aTGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtt                                                                                                                      <  2:184106/150‑1 (MQ=255)

                             ggTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGtttt                                                                                                                    <  1:746208/150‑1 (MQ=255)

                             ggTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtttt                                                                                                                    >  1:182899/1‑150 (MQ=255)

                               ttATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTATCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTg                                                                                                                  <  2:2078498/150‑1 (MQ=255)

                               ttATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTg                                                                                                                  >  2:1873977/1‑150 (MQ=255)

                                tATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGc                                                                                                                 >  1:2129057/1‑150 (MQ=255)

                                  tCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAt                                                                                                               >  2:2032877/1‑150 (MQ=255)

                                   ccaccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtt                                                                                                              <  1:547611/150‑1 (MQ=255)

                                    caccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttt                                                                                                             <  2:1724106/150‑1 (MQ=255)

                                     accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttt                                                                                                            >  1:2318522/1‑150 (MQ=255)

                                     accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttt                                                                                                            <  1:11996/150‑1 (MQ=255)

                                        aGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                         <  1:2627261/150‑1 (MQ=255)

                                        aGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                         <  1:984467/150‑1 (MQ=255)

                                        aGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                         <  1:2627530/150‑1 (MQ=255)

                                         gTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                                                                                                        <  1:3067573/150‑1 (MQ=255)

                                          tAATATGGCTGCCAGAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgtt                                                                                                       <  2:3902242/150‑1 (MQ=255)

                                             tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                                                                                                    >  2:334280/1‑150 (MQ=255)

                                             tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGaa                                                                                                    <  2:3733500/150‑1 (MQ=255)

                                              atGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                                                                                                   <  2:2136143/150‑1 (MQ=255)

                                                ggCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTg                                                                                                 >  1:3293905/1‑150 (MQ=255)

                                                  cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                                                                                               >  2:174807/1‑150 (MQ=255)

                                                   tGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGt                                                                                              <  1:4087723/150‑1 (MQ=255)

                                                   tGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGt                                                                                              <  1:3688072/150‑1 (MQ=255)

                                                     ccATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTga                                                                                            >  1:868995/1‑150 (MQ=255)

                                                      cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgca                                                                                           <  1:553793/150‑3 (MQ=255)

                                                      cATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgag                                                                                           <  1:887646/150‑1 (MQ=255)

                                                       aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgaga                                                                                          >  1:1237356/1‑150 (MQ=255)

                                                        tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGaa                                                                                         >  1:943424/1‑150 (MQ=255)

                                                           cacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtt                                                                                      >  2:791036/1‑150 (MQ=255)

                                                            acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                     <  1:2073267/150‑1 (MQ=255)

                                                            acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                     <  2:2201965/150‑1 (MQ=255)

                                                            acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                     >  2:3107698/1‑150 (MQ=255)

                                                            acTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGGTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                     >  2:3106451/1‑150 (MQ=255)

                                                             cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtttt                                                                                    >  1:309547/1‑150 (MQ=255)

                                                             cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtttt                                                                                    <  1:2466998/150‑1 (MQ=255)

                                                             cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtttt                                                                                    >  1:1090068/1‑150 (MQ=255)

                                                                aTGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCa                                                                                 >  2:2497611/1‑150 (MQ=255)

                                                                    cATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTg                                                                             >  2:1666721/1‑150 (MQ=255)

                                                                    cATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTg                                                                             >  1:2991745/1‑150 (MQ=255)

                                                                      tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGaa                                                                           <  2:1390766/150‑1 (MQ=255)

                                                                         gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAAGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAagag                                                                        >  1:1300595/1‑150 (MQ=255)

                                                                          aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGt                                                                       <  2:3253784/150‑1 (MQ=255)

                                                                              cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatc                                                                   >  2:2530076/1‑150 (MQ=255)

                                                                              cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatc                                                                   >  2:2530002/1‑150 (MQ=255)

                                                                                tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatcat                                                                 <  1:3789674/150‑1 (MQ=255)

                                                                                  ttGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCt                                                               <  2:1990910/150‑1 (MQ=255)

                                                                                    gTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTg                                                             <  1:3635006/150‑1 (MQ=255)

                                                                                     tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGa                                                            <  1:3430779/150‑1 (MQ=255)

                                                                                       aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATg                                                          <  2:4021102/150‑1 (MQ=255)

                                                                                       aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATg                                                          <  2:4019549/150‑1 (MQ=255)

                                                                                          atatGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATgggg                                                       <  2:2409999/150‑1 (MQ=255)

                                                                                             tGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTTAAGTGCGTGAGAAATTTTTAAGCTGAAGAGTAATAATCTTGATGGGGttt                                                    >  2:3499239/1‑150 (MQ=255)

                                                                                             tGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGttt                                                    <  2:189514/150‑1 (MQ=255)

                                                                                              gCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTg                                                   <  2:4331478/150‑1 (MQ=255)

                                                                                                   cATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                              <  1:3952388/150‑1 (MQ=255)

                                                                                                   cATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                              <  2:4246986/150‑1 (MQ=255)

                                                                                                    aTAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCa                                             <  1:2564297/150‑1 (MQ=255)

                                                                                                    aTAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCa                                             >  1:4409745/1‑150 (MQ=255)

                                                                                                    aTAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCa                                             <  2:4120468/150‑1 (MQ=255)

                                                                                                        cTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACAt                                         >  1:3405804/1‑150 (MQ=255)

                                                                                                          aGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTc                                       <  1:2881779/150‑1 (MQ=255)

                                                                                                            cccGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTg                                     >  1:669061/1‑150 (MQ=255)

                                                                                                             ccGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGa                                    >  2:2790483/1‑150 (MQ=255)

                                                                                                              cGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGAt                                   >  2:2730051/1‑150 (MQ=255)

                                                                                                              cGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGAt                                   <  1:1938120/150‑1 (MQ=255)

                                                                                                              cGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGAt                                   >  2:4286471/1‑150 (MQ=255)

                                                                                                               gAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATa                                  >  2:1805900/1‑150 (MQ=255)

                                                                                                                 aaTCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATGTTTTGTTGAAGGGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATaaa                                >  1:1155316/1‑150 (MQ=255)

                                                                                                                   tCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAata                              >  2:1438131/1‑150 (MQ=255)

                                                                                                                    cACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAatat                             >  2:889428/1‑150 (MQ=255)

                                                                                                                         caggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTa                        >  1:3139182/1‑150 (MQ=255)

                                                                                                                         caggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTa                        >  1:2761565/1‑150 (MQ=255)

                                                                                                                          aggcagCAACGCAACGCCCCTCAGGACGGCTGGACACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTaa                       <  2:1994019/150‑1 (MQ=255)

                                                                                                                          aggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTaa                       <  1:583282/150‑1 (MQ=255)

                                                                                                                          aggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTaa                       >  1:1763212/1‑150 (MQ=255)

                                                                                                                                   cgcaacgcCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATtgtg              >  2:47296/1‑150 (MQ=255)

                                                                                                                                      aacgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcag           >  2:4502638/1‑150 (MQ=255)

                                                                                                                                        cgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcagca         >  1:3596328/1‑150 (MQ=255)

                                                                                                                                           cccTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtg      >  1:3934389/1‑150 (MQ=255)

CTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTC  >  CP009273/4211657‑4211949

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