BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,590,074	A→C	13.2%	V240G (GTA→GGA)	dnaC ←	DNA biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,590,074	0	A	C	13.2%	583.1 / 91.7	234	V240G (GTA→GGA)	dnaC	DNA biosynthesis protein
Reads supporting (aligned to +/- strand): ref base A (127/76); new base C (23/8); total (150/84)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 236 total aligned reads displayed.

AGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCG  >  CP009273/4589925‑4590222

aGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTa                                                                                                                                                      >  1:566016/1‑150 (MQ=255)

  gACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTAcc                                                                                                                                                    <  1:4276125/150‑1 (MQ=255)

   aCGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTAccc                                                                                                                                                   >  1:4080801/1‑150 (MQ=255)

        gcaTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTa                                                                                                                                              <  1:4126033/150‑1 (MQ=255)

          aTGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACg                                                                                                                                            >  2:459309/1‑150 (MQ=255)

             cACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTa                                                                                                                                         <  2:1699885/150‑1 (MQ=255)

                cccGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCt                                                                                                                                      >  1:439700/1‑150 (MQ=255)

                cccGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCt                                                                                                                                      >  1:2727586/1‑150 (MQ=255)

                  cGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAt                                                                                                                                    <  2:2694910/150‑1 (MQ=255)

                  cGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTAt                                                                                                                                    <  1:2106841/150‑1 (MQ=255)

                   ggTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATAACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATc                                                                                                                                   >  1:1675458/1‑150 (MQ=255)

                    gTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATcc                                                                                                                                  >  1:3968501/1‑150 (MQ=255)

                    gTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATcc                                                                                                                                  >  1:3968532/1‑150 (MQ=255)

                     tGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATccc                                                                                                                                 >  2:2107908/1‑150 (MQ=255)

                     tGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATccc                                                                                                                                 >  2:4353361/1‑150 (MQ=255)

                      ggAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCa                                                                                                                                >  1:205931/1‑150 (MQ=255)

                        aaaTTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGt                                                                                                                              <  2:254459/150‑1 (MQ=255)

                          aTTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTACATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGATg                                                                                                                            >  1:1465834/1‑150 (MQ=255)

                          aTTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGATTAACAGTATTCATCAACATTACTCAATGTGGAATTGTGTCGCAGAATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACTGTAGCTATCCCAGTTg                                                                                                                            >  2:1143820/1‑150 (MQ=255)

                            tAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGaa                                                                                                                          >  1:2600506/1‑150 (MQ=255)

                             aaTGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAg                                                                                                                         >  1:855772/1‑150 (MQ=255)

                             aaTGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAg                                                                                                                         >  2:429440/1‑150 (MQ=255)

                             aaTGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAg                                                                                                                         >  2:429492/1‑150 (MQ=255)

                              aTGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGa                                                                                                                        >  1:328801/1‑150 (MQ=255)

                                   cTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAcc                                                                                                                   <  2:3218458/150‑1 (MQ=255)

                                    tGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCAccc                                                                                                                  >  1:1197629/1‑150 (MQ=255)

                                     ggcggcAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCca                                                                                                                 >  1:1984819/1‑150 (MQ=255)

                                         gcAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACaaa                                                                                                             >  2:31423/1‑150 (MQ=255)

                                         gcAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACaaa                                                                                                             >  2:911114/1‑150 (MQ=255)

                                             tGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGt                                                                                                         >  2:3821710/1‑150 (MQ=255)

                                               gcacagcacaGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTa                                                                                                       <  1:3408938/150‑1 (MQ=255)

                                               gcacagcacaGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTa                                                                                                       >  1:3174700/1‑150 (MQ=255)

                                                   agcacagAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCa                                                                                                   >  2:3685334/1‑150 (MQ=255)

                                                    gcacagAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAg                                                                                                  <  2:3348482/150‑1 (MQ=255)

                                                     cacagAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAgg                                                                                                 <  1:2575368/150‑1 (MQ=255)

                                                        agaagaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAgg                                                                                                 >  2:2657242/1‑147 (MQ=255)

                                                          aagaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCAt                                                                                            <  1:576303/150‑1 (MQ=255)

                                                          aagaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCAt                                                                                            >  1:1941184/1‑150 (MQ=255)

                                                          aagaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCAt                                                                                            >  1:1941125/1‑150 (MQ=255)

                                                            gaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATAc                                                                                          >  2:135791/1‑150 (MQ=255)

                                                            gaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATAc                                                                                          >  2:2369895/1‑150 (MQ=255)

                                                            gaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATAc                                                                                          >  1:1182590/1‑150 (MQ=255)

                                                            gaTGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATAc                                                                                          <  1:116827/150‑1 (MQ=255)

                                                             atctttCACAGTTTTCATCCTCATGACTCAAGGGGGAATTGTGTCGCAGTATCCCGCCAGACACGGGGATTAATACTCTCTACCTGGTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACg                                                                                         <  2:2504190/147‑1 (MQ=255)

                                                                         tttCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAcgcg                                                                             >  1:3441524/1‑150 (MQ=255)

                                                                         tttCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCAcgcg                                                                             >  1:3441735/1‑150 (MQ=255)

                                                                           tcatcatcaTTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCTCATACGGTCCATCACGCGtt                                                                           <  2:2904544/150‑1 (MQ=255)

                                                                           tcatcatcaTTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGtt                                                                           <  2:4101049/150‑1 (MQ=255)

                                                                           tcatcatcaTTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGtt                                                                           <  2:4223012/150‑1 (MQ=255)

                                                                               catcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                   >  1:1027426/1‑138 (MQ=255)

                                                                               catcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCAt                                                                                   <  2:1027426/138‑1 (MQ=255)

                                                                               catcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGcc                                                                       <  2:125338/150‑1 (MQ=255)

                                                                               catcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGcc                                                                       >  1:176892/1‑150 (MQ=255)

                                                                               catcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGcc                                                                       >  2:1690877/1‑150 (MQ=255)

                                                                                atcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGccc                                                                      >  2:2898035/1‑150 (MQ=255)

                                                                                atcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGccc                                                                      >  2:3435049/1‑150 (MQ=255)

                                                                                 tcatTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCa                                                                     >  2:865028/1‑150 (MQ=255)

                                                                                   atTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGc                                                                   >  1:2152128/1‑150 (MQ=255)

                                                                                      aCTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAAc                                                                >  2:2737281/1‑150 (MQ=255)

                                                                                       cTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACt                                                               <  2:3341766/150‑1 (MQ=255)

                                                                                             gTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCa                                                         >  1:3325189/1‑150 (MQ=255)

                                                                                              tGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTTATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCTTAAGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAt                                                        <  2:3616339/150‑1 (MQ=255)

                                                                                              tGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAt                                                        <  2:3616437/150‑1 (MQ=255)

                                                                                                gAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtct                                                      <  2:4498259/150‑1 (MQ=255)

                                                                                                 aaTTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCAtctc                                                     >  1:791960/1‑150 (MQ=255)

                                                                                                  aTTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCACCTTGGTCATctct                                                    >  2:1596477/1‑150 (MQ=255)

                                                                                                  aTTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATctct                                                    <  1:2820708/150‑1 (MQ=255)

                                                                                                    tgtgTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTc                                                  <  2:1972329/150‑1 (MQ=255)

                                                                                                    tgtgTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTc                                                  <  2:525936/150‑1 (MQ=255)

                                                                                                      tgtCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCa                                                <  2:1898940/150‑1 (MQ=255)

                                                                                                          gCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATAtt                                            <  2:3898095/150‑1 (MQ=255)

                                                                                                            aGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTAc                                          >  1:3902588/1‑150 (MQ=255)

                                                                                                            aGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTAc                                          >  1:525569/1‑150 (MQ=255)

                                                                                                             gTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACt                                         >  2:1463324/1‑150 (MQ=255)

                                                                                                             gTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACt                                         >  1:3805823/1‑150 (MQ=255)

                                                                                                             gTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACt                                         >  1:1546383/1‑150 (MQ=255)

                                                                                                             gTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTAACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACt                                         >  2:1462990/1‑150 (MQ=255)

                                                                                                                   cGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGt                                   <  1:4353361/150‑1 (MQ=255)

                                                                                                                   cGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGt                                   <  1:2798171/150‑1 (MQ=255)

                                                                                                                      cAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAg                                >  2:2920006/1‑150 (MQ=255)

                                                                                                                      cAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAg                                >  2:1944814/1‑150 (MQ=255)

                                                                                                                         aaaCGGGGATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCAt                             >  2:2155589/1‑150 (MQ=255)

                                                                                                                            cGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATccc                          >  1:3143825/1‑150 (MQ=255)

                                                                                                                             ggggATTAATACTCTTTACCTGTTCCCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCg                         <  1:3849532/150‑1 (MQ=255)

                                                                                                                                 aTTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCTAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTg                     <  2:1866701/150‑1 (MQ=255)

                                                                                                                                 aTTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTg                     <  2:3042525/150‑1 (MQ=255)

                                                                                                                                  ttAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTgg                    >  1:979247/1‑150 (MQ=255)

                                                                                                                                  ttAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTgg                    <  1:1099621/150‑1 (MQ=255)

                                                                                                                                    aaTACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGGCCATCACGCGTTCGCCCAGCCACTGGGTCATCTCTTCGCTATTTCGGTTGGTCAGCGTCGCGGGTGGGc                  >  2:2624243/1‑150 (MQ=255)

                                                                                                                                     aTACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCg                 >  1:3978411/1‑150 (MQ=255)

                                                                                                                                      tACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGt                <  1:3646470/150‑1 (MQ=255)

                                                                                                                                        ctctTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGttt              <  1:4363785/150‑1 (MQ=255)

                                                                                                                                                cTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCCACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGaagaag      >  2:4359512/1‑150 (MQ=255)

                                                                                                                                                cTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGaagaag      >  1:3087207/1‑150 (MQ=255)

                                                                                                                                                 tGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAATCTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTGGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAagaaga     >  2:2496721/1‑150 (MQ=255)

                                                                                                                                                 tGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTAACCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAagaaga     >  2:197897/1‑150 (MQ=255)

                                                                                                                                                  gTTACCCGGCTACGGTCGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGATTGGAAGAagag    >  2:2797713/1‑150 (MQ=255)

                                                                                                                                                  gTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAagag    >  1:3213646/1‑150 (MQ=255)

                                                                                                                                                  gTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAagag    >  1:3213323/1‑150 (MQ=255)

                                                                                                                                                    tACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTTGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAgcg  <  2:1182630/150‑1 (MQ=255)

AGGACGCAGCATGCACCCCGGTGGAAATTAATGCGCTGGCGGCAATGGCACAGCACAGAAGATGTTTAACAGTTTTCATCATCATTACTCAAGGTGGAATTGTGTCGCAGTATACCGCCAGAAACGGGGATTAATACTCTTTACCTGTTACCCGGCTACGGTAGCTATCCCAGTTGAAGATCACCCACAAACTGTTACCCAGGCGCATACGGTCCATCACGCGTTCGCCCAGCAACTTGGTCATCTCTTCCATATTACTGTTGGTCAGCATCCCGGTTGGGCGTTTGGAAGAAGAGCG  >  CP009273/4589925‑4590222

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