BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,211,803	G→A	45.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	45.4%	142.5 / inf	351	P205S (CCT→TCT)	arpA	ankyrin repeat protein
Reads supporting (aligned to +/- strand): ref base G (101/90); new base A (72/87); total (174/177)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 355 total aligned reads displayed.

GTTCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCGGC  >  CP009273/4211654‑4211952

gtTCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGATTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAgg                                                                                                                                                       >  2:4410281/1‑150 (MQ=255)

  tCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGAc                                                                                                                                                     <  1:2848492/150‑1 (MQ=255)

  tCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAAc                                                                                                                                                     <  2:1427051/150‑1 (MQ=255)

  tCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAAc                                                                                                                                                     <  2:1696350/150‑1 (MQ=255)

  tCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCAGAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGAc                                                                                                                                                     <  2:1495644/150‑1 (MQ=255)

      tATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCt                                                                                                                                                 <  2:2949734/150‑1 (MQ=255)

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

            tAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGACTGGAAAc                                                                                                                                           >  1:1431565/1‑150 (MQ=255)

            tAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGACTGGAAAc                                                                                                                                           >  1:3491479/1‑150 (MQ=255)

             aGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAAcc                                                                                                                                          <  1:4310885/150‑1 (MQ=255)

              gCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCa                                                                                                                                         <  2:5013641/150‑1 (MQ=255)

                cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                                                                                                       >  2:722333/1‑150 (MQ=255)

                cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAAc                                                                                                                                       <  2:2561037/150‑1 (MQ=255)

                cGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAAc                                                                                                                                       >  1:2722845/1‑150 (MQ=255)

                     cGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCaa                                                                                                                                  >  1:4140859/1‑150 (MQ=255)

                       tCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAg                                                                                                                                <  1:2445165/150‑1 (MQ=255)

                        cAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACAGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGt                                                                                                                               >  2:3723765/1‑150 (MQ=255)

                         aaaGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTc                                                                                                                              >  2:2937775/1‑150 (MQ=255)

                         aaaGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTc                                                                                                                              <  1:2497357/150‑1 (MQ=255)

                           aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTg                                                                                                                            <  1:2534075/150‑1 (MQ=255)

                            gCAGGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTgg                                                                                                                           <  2:2797851/150‑1 (MQ=255)

                             cATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGt                                                                                                                          >  1:80836/1‑150 (MQ=255)

                              aTGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGtt                                                                                                                         >  1:1051276/1‑150 (MQ=255)

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

                                 gTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                                                                                                      >  1:4119457/1‑150 (MQ=255)

                                 gTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGAGTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                                                                                                      <  2:2345898/150‑1 (MQ=255)

                                    ttCCACCAGTTATATTGCTGACATAACACTGATGACATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCa                                                                                                                   <  2:3226885/149‑1 (MQ=255)

                                        accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttt                                                                                                               <  2:448897/150‑1 (MQ=255)

                                            gTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                                                                                                           <  1:4608130/150‑1 (MQ=255)

                                              aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAAATAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCATGACGGCTGGAAACCAACGCCAAGGTCTGGTTTTTGCATTTTttgttg                                                                                                         >  2:4386304/1‑150 (MQ=255)

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

                                                 atGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                                                                                                      >  1:1522667/1‑150 (MQ=255)

                                                 atGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                                                                                                      >  1:1642041/1‑150 (MQ=255)

                                                     cTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                                                                                                  <  2:5140107/150‑1 (MQ=255)

                                                       gCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTg                                                                                                <  1:2620188/150‑1 (MQ=255)

                                                          aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgaga                                                                                             >  2:1262444/1‑150 (MQ=255)

                                                           tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGACGTGCGTGAGaa                                                                                            >  2:4195434/1‑150 (MQ=255)

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

                                                           tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGaa                                                                                            <  2:5031608/150‑1 (MQ=255)

                                                              cacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtt                                                                                         <  1:4001121/150‑1 (MQ=255)

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

                                                               acTGAGGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGTTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                        <  1:1584374/150‑1 (MQ=255)

                                                                  gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc                                                                                     >  2:297227/1‑150 (MQ=255)

                                                                   aTGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGTCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCa                                                                                    <  1:4042257/150‑1 (MQ=255)

                                                                   aTGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCa                                                                                    <  2:1874032/150‑1 (MQ=255)

                                                                    tGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAg                                                                                   <  1:3737024/150‑1 (MQ=255)

                                                                      ccATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCt                                                                                 >  2:1102866/1‑150 (MQ=255)

                                                                      ccATTTGATTCCATCTTGTCAAAATATGATTTCAAAACTAGCCCGAAATCACCGCAGGCAGCAACGAAACTCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCTTGAGAAATTTGTCAGCt                                                                                 >  1:4829322/1‑150 (MQ=255)

                                                                         tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGaa                                                                              <  1:1262444/150‑1 (MQ=255)

                                                                         tttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGaa                                                                              >  1:2299228/1‑150 (MQ=255)

                                                                          ttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAag                                                                             >  2:1822575/1‑150 (MQ=255)

                                                                          ttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAag                                                                             >  1:1707874/1‑150 (MQ=255)

                                                                            gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAagag                                                                           <  1:516629/150‑1 (MQ=255)

                                                                            gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAagag                                                                           >  1:2373367/1‑150 (MQ=255)

                                                                              ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTa                                                                         <  2:921558/150‑1 (MQ=255)

                                                                              ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTa                                                                         >  1:1989259/1‑150 (MQ=255)

                                                                                 cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatc                                                                      >  1:3364235/1‑150 (MQ=255)

                                                                                  aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatca                                                                     <  2:4641349/150‑1 (MQ=255)

                                                                                      tGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCtt                                                                 <  1:193584/150‑1 (MQ=255)

                                                                                      tGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCtt                                                                 <  2:4067543/150‑1 (MQ=255)

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

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

                                                                                        tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGa                                                               >  2:390149/1‑150 (MQ=255)

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

                                                                                                tGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGttt                                                       >  2:2045635/1‑150 (MQ=255)

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

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

                                                                                                  cTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGt                                                     >  2:4595726/1‑150 (MQ=255)

                                                                                                  cTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGt                                                     >  2:3938625/1‑150 (MQ=255)

                                                                                                  cTTACATAACTAGCCCGAAATCACCGCAGGCAGCAACGAAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTATTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGt                                                     >  1:4343016/1‑150 (MQ=255)

                                                                                                   tttCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTa                                                    <  1:3643047/150‑1 (MQ=255)

                                                                                                     tCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAAt                                                  >  1:3401370/1‑150 (MQ=255)

                                                                                                       aTAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCa                                                <  2:1636517/150‑1 (MQ=255)

                                                                                                       aTAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCa                                                <  1:546639/150‑1 (MQ=255)

                                                                                                        tAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCaa                                               >  2:1521434/1‑150 (MQ=255)

                                                                                                         aaCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAAc                                              >  1:3107037/1‑150 (MQ=255)

                                                                                                         aaCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAAc                                              >  1:3483770/1‑150 (MQ=255)

                                                                                                            tAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACAtt                                           >  2:3434813/1‑150 (MQ=255)

                                                                                                            tAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACAtt                                           >  2:2081389/1‑150 (MQ=255)

                                                                                                               cccGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTg                                        >  2:1046982/1‑150 (MQ=255)

                                                                                                                ccGAAATCACCGCAGGCAGCAACGCAACGACCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGa                                       <  2:2275994/150‑1 (MQ=255)

                                                                                                                  gAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGACACCAACGCCCATGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATa                                     >  2:1540627/1‑150 (MQ=255)

                                                                                                                  gAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATa                                     <  1:898759/150‑1 (MQ=255)

                                                                                                                     aTCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAat                                  >  1:4565540/1‑150 (MQ=255)

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

                                                                                                                       cACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGACGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAatat                                >  1:1231268/1‑150 (MQ=255)

                                                                                                                       cACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAatat                                >  2:3339896/1‑150 (MQ=255)

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

                                                                                                                              ggcagCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAat                         <  2:574490/150‑1 (MQ=255)

                                                                                                                               gcagcaACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAata                        <  1:1822748/150‑1 (MQ=255)

                                                                                                                                 agcaACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAatata                      <  2:37408/150‑1 (MQ=255)

                                                                                                                                 agcaACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAatata                      >  2:3022612/1‑150 (MQ=255)

                                                                                                                                   caacgcaacgCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATAtt                    >  1:4956010/1‑150 (MQ=255)

                                                                                                                                      cgcaacgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATtgtg                 <  2:5204166/150‑1 (MQ=255)

                                                                                                                                         aacgcCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcag              <  2:1593572/150‑1 (MQ=255)

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

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

                                                                                                                                          acgcCCCTCAGGACGGCTGGAAAACAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATATTGATGGGGTTTGTAATCAACATTATGATATATATAATTAATATATTGTgcagc             >  2:2719897/1‑150 (MQ=255)

                                                                                                                                                 tCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATgtgt      <  1:3841746/150‑1 (MQ=255)

                                                                                                                                                   aGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCg    >  1:841273/1‑150 (MQ=255)

GTTCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCGGC  >  CP009273/4211654‑4211952

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