BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,211,803	G→A	92.1%	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	92.1%	999.8 / 57.4	290	P205S (CCT→TCT)	arpA	ankyrin repeat protein
Reads supporting (aligned to +/- strand): ref base G (9/14); new base A (132/135); total (141/149)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 290 total aligned reads displayed.

TCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCG  >  CP009273/4211656‑4211950

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

     aTTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTg                                                                                                                                              >  2:1921378/1‑150 (MQ=255)

      ttCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTgg                                                                                                                                             >  2:667994/1‑150 (MQ=255)

       tCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGa                                                                                                                                            >  2:2168316/1‑150 (MQ=255)

         ttAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGaaa                                                                                                                                          >  1:4243903/1‑150 (MQ=255)

            gCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCa                                                                                                                                       >  1:4917824/1‑150 (MQ=255)

             ccGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCaa                                                                                                                                      >  1:4177362/1‑150 (MQ=255)

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

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

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

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

                 aCCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGcc                                                                                                                                  <  2:2741743/150‑1 (MQ=255)

                    gTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCaca                                                                                                                               >  2:2651242/1‑148 (MQ=255)

                    gTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCaaa                                                                                                                               <  1:3347413/150‑1 (MQ=255)

                        aaGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCt                                                                                                                           <  2:2043113/150‑1 (MQ=255)

                         aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTg                                                                                                                          >  1:5006863/1‑150 (MQ=255)

                           cATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGt                                                                                                                        <  1:2457226/150‑1 (MQ=255)

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

                             tGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGttt                                                                                                                      <  1:4862618/150‑1 (MQ=255)

                               gTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                                                                                                    >  1:2840243/1‑150 (MQ=255)

                               gTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                                                                                                    <  1:4471316/150‑1 (MQ=255)

                               gTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGttttt                                                                                                                    >  1:3934883/1‑150 (MQ=255)

                                 tATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGttttttc                                                                                                                  <  2:3085129/150‑3 (MQ=255)

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

                                       ccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttttt                                                                                                            >  1:3580103/1‑150 (MQ=255)

                                           tAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgtt                                                                                                        <  2:3780094/150‑1 (MQ=255)

                                            aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgttg                                                                                                       >  2:2960498/1‑150 (MQ=255)

                                            aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgttg                                                                                                       <  1:3256205/150‑1 (MQ=255)

                                                tGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGt                                                                                                   >  1:3609980/1‑150 (MQ=255)

                                                tGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGGTTTTGCATTTTTTGTTGAAGt                                                                                                   >  2:4466678/1‑150 (MQ=255)

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

                                                 ggCTGCCATAAAACTGATGCCATTTGAATCCATCTTGTAAAAATATGCATTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCAGGTATTTGCATTTTTTGTTGAAGtc                                                                                                  >  2:2505702/1‑149 (MQ=255)

                                                    tGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGt                                                                                               >  1:1774659/1‑150 (MQ=255)

                                                    tGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGt                                                                                               >  1:3081777/1‑150 (MQ=255)

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

                                                        aTAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgaga                                                                                           <  1:1462571/150‑1 (MQ=255)

                                                          aaCACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGaaa                                                                                         <  1:3790776/150‑1 (MQ=255)

                                                           acacTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAt                                                                                        <  2:2328062/150‑1 (MQ=255)

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

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

                                                                gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc                                                                                   >  2:4351991/1‑150 (MQ=255)

                                                                gATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTc                                                                                   <  1:3989174/150‑1 (MQ=255)

                                                                 aTGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCa                                                                                  <  2:5022094/150‑1 (MQ=255)

                                                                  tGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAg                                                                                 >  2:2465222/1‑150 (MQ=255)

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

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

                                                                            ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTa                                                                       <  1:2553775/150‑1 (MQ=255)

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

                                                                                      tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGa                                                             >  2:1254297/1‑150 (MQ=255)

                                                                                       cAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGAt                                                            <  1:656949/150‑1 (MQ=255)

                                                                                        aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATg                                                           >  2:1900269/1‑150 (MQ=255)

                                                                                            tatGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGt                                                       <  1:3428954/150‑1 (MQ=255)

                                                                                             atGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGtt                                                      >  2:749355/1‑150 (MQ=255)

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

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

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

                                                                                               gCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTg                                                    >  1:3189118/1‑150 (MQ=255)

                                                                                                cTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGt                                                   <  2:2147031/150‑1 (MQ=255)

                                                                                                 tttCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTa                                                  <  1:904171/150‑1 (MQ=255)

                                                                                                    cATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                               >  1:1291915/1‑150 (MQ=255)

                                                                                                     aTAACTAGCAAGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAAACAACGCCAAAGTCTGGATATTGCATTATTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATAATCTTGAGGGGGTTTGTAAtaa                                              >  2:4704768/1‑148 (MQ=255)

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

                                                                                                        aCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACa                                           <  1:3981635/150‑1 (MQ=255)

                                                                                                        aCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACa                                           <  1:1374461/150‑1 (MQ=255)

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

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

                                                                                                            gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCt                                       <  1:1108544/150‑1 (MQ=255)

                                                                                                            gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCt                                       <  1:4527920/150‑1 (MQ=255)

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

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

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

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

                                                                                                                 aaaTCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCACAGGCTGGTGTTTGCCTGTTGTGTTGAAGTGCGTGAGAACTTTTTCAGCTGACGAGTCCTCATCGTGAGGGGGTTTGTAATCAACCTTCTGAtca                                  >  2:213349/1‑148 (MQ=255)

                                                                                                                 aaaTCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATaa                                  >  2:239722/1‑150 (MQ=255)

                                                                                                                 aaaTCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATaa                                  >  2:3496969/1‑150 (MQ=255)

                                                                                                                    tCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGGTTTTGCATTGTGTGTTGAAGTGCGTGAGAAATTTTTCCGCTGAAGAGTAATCATCTTGAGGGGGTGTGGACTCAACATTCTGATaactc                               >  2:2964382/1‑147 (MQ=255)

                                                                                                                      aCCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGGTTGTAATCAACATTCTGATAAATATc                             >  2:4404135/1‑150 (MQ=255)

                                                                                                                       ccGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCa                            <  1:2020040/150‑1 (MQ=255)

                                                                                                                       ccGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCa                            <  1:549161/150‑1 (MQ=255)

                                                                                                                       ccGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCa                            <  1:3586110/150‑1 (MQ=255)

                                                                                                                        cGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCAt                           >  2:2580488/1‑150 (MQ=255)

                                                                                                                        cGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCAt                           <  2:2208039/150‑1 (MQ=255)

                                                                                                                         gcaggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCAtt                          >  2:1582067/1‑150 (MQ=255)

                                                                                                                         gcaggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCAtt                          <  1:5058011/150‑1 (MQ=255)

                                                                                                                         gcaggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCAtt                          >  2:3405626/1‑150 (MQ=255)

                                                                                                                             gcagcaACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAata                      <  1:2196489/150‑1 (MQ=255)

                                                                                                                               agcaACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAatata                    <  1:81780/150‑1 (MQ=255)

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

                                                                                                                                      caacgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgca             >  2:4215016/1‑150 (MQ=255)

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

                                                                                                                                        acgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcagc           <  2:4144419/150‑1 (MQ=255)

                                                                                                                                          gcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCaa         >  1:4454336/1‑150 (MQ=255)

                                                                                                                                           ccccTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAt        <  2:3674096/150‑1 (MQ=255)

                                                                                                                                            cccTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTGTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtg       >  2:4905390/1‑150 (MQ=255)

                                                                                                                                              cTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgtg     <  1:870885/150‑1 (MQ=255)

                                                                                                                                              cTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgtg     >  1:593884/1‑150 (MQ=255)

                                                                                                                                               tCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATgtgt    <  2:169535/150‑1 (MQ=255)

                                                                                                                                                cAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTc   <  1:2960498/150‑1 (MQ=255)

                                                                                                                                                cAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTc   <  1:3096586/150‑1 (MQ=255)

                                                                                                                                                 aGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCg  <  2:182727/150‑1 (MQ=255)

TCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTGTCG  >  CP009273/4211656‑4211950

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