BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,688,819	C→A	28.4%	E504* (GAA→TAA)	uidA ←	beta‑D‑glucuronidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,688,819	0	C	A	28.4%	617.4 / 202.6	305	E504* (GAA→TAA)	uidA	beta‑D‑glucuronidase
Reads supporting (aligned to +/- strand): ref base C (95/123); new base A (43/43); total (138/167)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 310 total aligned reads displayed.

ATCGGCGAAATTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGATGGTAT  >  CP009273/1688670‑1688968

 tCGGCGAAATTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCg                                                                                                                                                      >  1:1044223/1‑150 (MQ=255)

  cGGCGAAATTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAgg                                                                                                                                                     <  2:2108523/150‑1 (MQ=255)

      gAAATTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGtgat                                                                                                                                                 <  1:1185068/150‑1 (MQ=255)

      gAAATTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGtgat                                                                                                                                                 <  1:4554222/150‑1 (MQ=255)

       aaaTTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGtgatg                                                                                                                                                >  1:155393/1‑150 (MQ=255)

         aTTCCATACCTGTTCACCGACGACGGCGCTGATGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTgatgat                                                                                                                                              <  1:2317941/150‑1 (MQ=255)

         aTTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTgatgat                                                                                                                                              <  2:4936557/150‑1 (MQ=255)

         aTTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTgatgat                                                                                                                                              <  2:104129/150‑1 (MQ=255)

          ttCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATa                                                                                                                                             >  1:907289/1‑150 (MQ=255)

          ttCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATa                                                                                                                                             <  2:117216/150‑1 (MQ=255)

           tCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATaa                                                                                                                                            >  2:2716933/1‑150 (MQ=255)

            ccATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAAt                                                                                                                                           <  1:834315/150‑1 (MQ=255)

                 ccTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGAGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCt                                                                                                                                      >  2:2762771/1‑150 (MQ=255)

                      tCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGc                                                                                                                                 >  1:3435163/1‑150 (MQ=255)

                           gacgacGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGttt                                                                                                                            >  2:4174014/1‑150 (MQ=255)

                            acgacgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTtc                                                                                                                           >  1:4415852/1‑150 (MQ=255)

                             cgacgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtct                                                                                                                          >  1:1836605/1‑150 (MQ=255)

                              gacgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtctc                                                                                                                         >  1:2419058/1‑150 (MQ=255)

                              gacgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtctc                                                                                                                         >  1:5309230/1‑150 (MQ=255)

                              gacgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtctc                                                                                                                         >  1:5223427/1‑150 (MQ=255)

                              gacgGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTtctc                                                                                                                         >  1:5066176/1‑150 (MQ=255)

                                 ggCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTg                                                                                                                      <  2:1112338/150‑1 (MQ=255)

                                   cgcTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTgcc                                                                                                                    <  1:2823331/150‑1 (MQ=255)

                                      tGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccagg                                                                                                                 >  1:1278228/1‑150 (MQ=255)

                                       gACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggc                                                                                                                >  1:2901487/1‑150 (MQ=255)

                                         cgcgATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTgccaggcca                                                                                                              <  1:5137925/150‑1 (MQ=255)

                                            gATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGaa                                                                                                           >  2:3023807/1‑150 (MQ=255)

                                            gATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGaa                                                                                                           >  2:3885594/1‑150 (MQ=255)

                                                 aaGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCt                                                                                                      <  2:2376203/150‑1 (MQ=255)

                                                 aaGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCt                                                                                                      <  2:3106025/150‑1 (MQ=255)

                                                 aaGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCt                                                                                                      <  2:1905382/150‑1 (MQ=255)

                                                  aGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCtt                                                                                                     <  1:4876869/150‑1 (MQ=255)

                                                  aGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCtt                                                                                                     >  1:2150668/1‑150 (MQ=255)

                                                   gACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCttt                                                                                                    <  2:3924188/150‑1 (MQ=255)

                                                      gcgGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTc                                                                                                 <  1:3672173/150‑1 (MQ=255)

                                                      gcgGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTc                                                                                                 <  2:4519276/150‑1 (MQ=255)

                                                       cgGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTcc                                                                                                <  2:1744869/150‑1 (MQ=255)

                                                        ggTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCa                                                                                               >  1:660205/1‑150 (MQ=255)

                                                        ggTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCa                                                                                               >  1:4975413/1‑150 (MQ=255)

                                                         gTGATACATCTCCAGCCATGCACACTGATACTCGTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCGAACGTATCCACGCCGGTTTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAg                                                                                              <  1:2627180/150‑1 (MQ=255)

                                                           gATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTa                                                                                            <  1:5178760/150‑1 (MQ=255)

                                                            aTACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTAc                                                                                           <  1:4496126/150‑1 (MQ=255)

                                                            aTACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTAc                                                                                           <  2:4842190/150‑1 (MQ=255)

                                                            aTACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTAc                                                                                           >  1:2726097/1‑150 (MQ=255)

                                                            aTACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTAc                                                                                           <  1:2491989/150‑1 (MQ=255)

                                                            aTACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTAc                                                                                           >  1:596550/1‑150 (MQ=255)

                                                                atatCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTtc                                                                                       <  1:1135687/150‑1 (MQ=255)

                                                                  atCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTtctc                                                                                     >  2:1851894/1‑150 (MQ=255)

                                                                   tCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTctct                                                                                    <  1:3993401/150‑1 (MQ=255)

                                                                   tCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTctct                                                                                    <  2:2611388/150‑1 (MQ=255)

                                                                   tCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTctct                                                                                    <  2:1659127/150‑1 (MQ=255)

                                                                    ccAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTg                                                                                   >  2:1414700/1‑150 (MQ=255)

                                                                     cAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGc                                                                                  <  1:963689/150‑1 (MQ=255)

                                                                      aGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGcc                                                                                 <  2:2039410/150‑1 (MQ=255)

                                                                       gCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGGATTAGGTGATGATAATCGGCTGAGGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCGTCTCTGCCg                                                                                >  2:5050177/1‑150 (MQ=255)

                                                                         cATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGtt                                                                              <  1:754355/150‑1 (MQ=255)

                                                                         cATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGtt                                                                              <  2:1908748/150‑1 (MQ=255)

                                                                         cATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGtt                                                                              <  1:3741415/150‑1 (MQ=255)

                                                                          aTGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGttt                                                                             <  1:1154503/150‑1 (MQ=255)

                                                                                acTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAAt                                                                       >  2:3791002/1‑150 (MQ=255)

                                                                                acTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGAGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGAACCTTCTCTGCCGTATCCAAAt                                                                       >  2:3979413/1‑150 (MQ=255)

                                                                                 cTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATc                                                                      <  2:3435163/150‑1 (MQ=255)

                                                                                  tGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATcg                                                                     >  1:2584424/1‑150 (MQ=255)

                                                                                   gATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATcgc                                                                    >  1:4304599/1‑150 (MQ=255)

                                                                                     tACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATcgccg                                                                  <  1:3111386/150‑1 (MQ=255)

                                                                                       ctctTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCt                                                                >  2:5100274/1‑150 (MQ=255)

                                                                                       ctctTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCt                                                                >  2:1299011/1‑150 (MQ=255)

                                                                                               tCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACAt                                                        <  2:1229499/150‑1 (MQ=255)

                                                                                               tCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACAt                                                        <  2:1577797/150‑1 (MQ=255)

                                                                                               tCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACAt                                                        >  1:1913537/1‑150 (MQ=255)

                                                                                                   caTGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCa                                                    >  1:3546645/1‑150 (MQ=255)

                                                                                                    aTGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCAt                                                   >  2:2265975/1‑150 (MQ=255)

                                                                                                     tGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATc                                                  >  1:1597221/1‑150 (MQ=255)

                                                                                                     tGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATc                                                  <  1:4586766/150‑1 (MQ=255)

                                                                                                         ggTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGta                                              >  2:3533772/1‑150 (MQ=255)

                                                                                                            gtACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGtaata                                           >  2:5150474/1‑150 (MQ=255)

                                                                                                            gtACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGtaata                                           <  1:1335275/150‑1 (MQ=255)

                                                                                                              aCATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAAc                                         <  1:315020/150‑1 (MQ=255)

                                                                                                                   gAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAATGGTTc                                    <  2:2194567/150‑1 (MQ=255)

                                                                                                                   gAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTc                                    <  2:193315/150‑1 (MQ=255)

                                                                                                                     gTGCGGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAg                                  <  1:4159695/150‑1 (MQ=255)

                                                                                                                      tGCAGCCCGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAgg                                 >  2:2812170/1‑150 (MQ=255)

                                                                                                                       gCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGgc                                <  1:2212910/150‑1 (MQ=255)

                                                                                                                             cGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGgcacagca                          >  2:2552139/1‑150 (MQ=255)

                                                                                                                             cGGCTAACGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGgcacagca                          >  1:3219299/1‑150 (MQ=255)

                                                                                                                              ggCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGgcacagcac                         <  1:1802386/150‑1 (MQ=255)

                                                                                                                                  aaCGTATCCACGCCGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCGTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCa                     <  2:4069152/150‑1 (MQ=255)

                                                                                                                                   aCGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCaa                    <  1:3791002/150‑1 (MQ=255)

                                                                                                                                       aTCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAagag                >  2:4878441/1‑150 (MQ=255)

                                                                                                                                        tCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAgaga               >  1:2031971/1‑150 (MQ=255)

                                                                                                                                         ccACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGAt              <  1:4789779/150‑1 (MQ=255)

                                                                                                                                         ccACGCCGTATTAGGAGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGAt              <  1:4032468/150‑1 (MQ=255)

                                                                                                                                              ccGTATTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTg         >  2:259898/1‑150 (MQ=255)

                                                                                                                                               cGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGa        <  2:3167302/150‑1 (MQ=255)

                                                                                                                                               cGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGa        >  2:127482/1‑150 (MQ=255)

                                                                                                                                               cGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTATTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGa        >  2:3644048/1‑150 (MQ=255)

                                                                                                                                                  aTTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGATgg     >  2:4380091/1‑150 (MQ=255)

                                                                                                                                                  aTTAGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGATgg     >  2:3346872/1‑150 (MQ=255)

                                                                                                                                                     cGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGATGGTAt  <  2:5223427/150‑1 (MQ=255)

                                                                                                                                                     cGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGATGGTAt  <  2:54839/150‑1 (MQ=255)

ATCGGCGAAATTCCATACCTGTTCACCGACGACGGCGCTGACGCGATCAAAGACGCGGTGATACATATCCAGCCATGCACACTGATACTCTTCACTCCACATGTCGGTGTACATTGAGTGCAGCCCGGCTAACGTATCCACGCCGTATTCGGTGATGATAATCGGCTGATGCAGTTTCTCCTGCCAGGCCAGAAGTTCTTTTTCCAGTACCTTCTCTGCCGTTTCCAAATCGCCGCTTTGGACATACCATCCGTAATAACGGTTCAGGCACAGCACATCAAAGAGATCGCTGATGGTAT  >  CP009273/1688670‑1688968

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