BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,237,098	G→T	52.6%	H12N (CAT→AAT)	folE ←	GTP cyclohydrolase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,237,098	0	G	T	52.6%	1.2 / inf	281	H12N (CAT→AAT)	folE	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand): ref base G (51/81); new base T (69/78); total (122/159)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 282 total aligned reads displayed.

CATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCACAG  >  CP009273/2236949‑2237246

cATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATg                                                                                                                                                      >  1:1438032/1‑150 (MQ=255)

   cAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAc                                                                                                                                                   <  2:933638/150‑1 (MQ=255)

   cAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAc                                                                                                                                                   <  2:2948458/150‑1 (MQ=255)

    aaaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAAcc                                                                                                                                                  >  1:720279/1‑150 (MQ=255)

      aCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAg                                                                                                                                                >  1:666181/1‑150 (MQ=255)

          tcatcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGcc                                                                                                                                            >  1:3267154/1‑150 (MQ=255)

          tcatcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGcc                                                                                                                                            >  1:3049040/1‑150 (MQ=255)

           catcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCg                                                                                                                                           <  2:3640673/150‑1 (MQ=255)

               aGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttc                                                                                                                                       >  1:257793/1‑150 (MQ=255)

                  cAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCttt                                                                                                                                    >  1:4238898/1‑150 (MQ=255)

                   aGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTa                                                                                                                                   <  2:5199963/150‑1 (MQ=255)

                       cGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGa                                                                                                                               >  2:748940/1‑150 (MQ=255)

                        gagaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAg                                                                                                                              <  2:4615534/150‑1 (MQ=255)

                        gagaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAg                                                                                                                              <  2:572133/150‑1 (MQ=255)

                         agaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGt                                                                                                                             <  1:5145648/150‑1 (MQ=255)

                          gaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTg                                                                                                                            <  2:2554003/150‑1 (MQ=255)

                           aTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGa                                                                                                                           >  2:3854007/1‑150 (MQ=255)

                           aTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGa                                                                                                                           <  1:3644156/150‑1 (MQ=255)

                             tCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATg                                                                                                                         <  2:1576991/150‑1 (MQ=255)

                             tCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATg                                                                                                                         <  1:4802911/150‑1 (MQ=255)

                                 cagcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCAt                                                                                                                     >  2:2397404/1‑150 (MQ=255)

                                  agcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAtt                                                                                                                    <  2:2171389/150‑1 (MQ=255)

                                   gcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAttt                                                                                                                   <  1:3797453/150‑1 (MQ=255)

                                     tGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTAt                                                                                                                 <  1:1273526/150‑1 (MQ=255)

                                     tGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTAt                                                                                                                 <  1:5002090/150‑1 (MQ=255)

                                       cATGAGTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGa                                                                                                               <  2:1525749/150‑1 (MQ=255)

                                         tGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAtt                                                                                                             <  1:4781939/150‑1 (MQ=255)

                                             ttCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTtctc                                                                                                         >  1:2981380/1‑150 (MQ=255)

                                                     tatGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgt                                                                                                 <  2:3321087/150‑1 (MQ=255)

                                                     tatGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgt                                                                                                 <  1:257911/150‑1 (MQ=255)

                                                       tGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGATTGATGGCATTTATGATTTCTCCTGCAGgtgt                                                                                               <  1:2683215/150‑1 (MQ=255)

                                                        gACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATTGCATGTATGATTTCTCCTGCAGGtgtg                                                                                              >  1:3188993/1‑150 (MQ=255)

                                                         aCCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGa                                                                                             <  2:2197324/150‑1 (MQ=255)

                                                         aCCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGa                                                                                             <  1:3238457/150‑1 (MQ=255)

                                                          ccAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGAc                                                                                            >  1:1762629/1‑150 (MQ=255)

                                                          ccAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGAc                                                                                            >  2:1988129/1‑150 (MQ=255)

                                                          ccAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGAc                                                                                            <  1:873810/150‑1 (MQ=255)

                                                            aGCAATAAGTCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGc                                                                                          >  1:601784/1‑150 (MQ=255)

                                                                    ggCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGccc                                                                                  >  1:2731275/1‑150 (MQ=255)

                                                                     gCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCt                                                                                 >  2:1855801/1‑150 (MQ=255)

                                                                      cTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTg                                                                                >  1:5081286/1‑150 (MQ=255)

                                                                      cTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTg                                                                                >  2:5167085/1‑150 (MQ=255)

                                                                          tGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTg                                                                            <  2:1978642/150‑1 (MQ=255)

                                                                            cgcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggg                                                                          >  2:2608901/1‑150 (MQ=255)

                                                                            cgcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggg                                                                          >  1:3173581/1‑150 (MQ=255)

                                                                             gcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTgggg                                                                         <  1:2171281/150‑1 (MQ=255)

                                                                             gcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTgggg                                                                         <  2:5215243/150‑1 (MQ=255)

                                                                               gTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGc                                                                       >  1:5168884/1‑150 (MQ=255)

                                                                                tttCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCa                                                                      >  1:3540062/1‑150 (MQ=255)

                                                                                  tCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCaaa                                                                    >  2:2642217/1‑150 (MQ=255)

                                                                                   cGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTTATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAg                                                                   <  2:2874730/150‑1 (MQ=255)

                                                                                   cGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAg                                                                   <  2:421873/150‑1 (MQ=255)

                                                                                   cGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAg                                                                   <  2:4848414/150‑1 (MQ=255)

                                                                                    gTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt                                                                  >  2:1058099/1‑150 (MQ=255)

                                                                                    gTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt                                                                  >  2:3769462/1‑150 (MQ=255)

                                                                                     ttATCCATTTCAGGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGtt                                                                 <  2:5082789/150‑1 (MQ=255)

                                                                                      tATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCATCATGAACCAGGGCCGCTTCTTTACAGAGTGATGGCATTTATGATTTATCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGtaa                                                                >  2:4170217/1‑148 (MQ=255)

                                                                                         ccATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGttatta                                                             <  2:1210480/150‑1 (MQ=255)

                                                                                         ccATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGttatta                                                             <  1:94261/150‑1 (MQ=255)

                                                                                         ccATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGttatta                                                             <  2:4372826/150‑1 (MQ=255)

                                                                                           aTTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAg                                                           <  1:249785/150‑1 (MQ=255)

                                                                                            tttCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGc                                                          >  1:1972904/1‑150 (MQ=255)

                                                                                              tCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTg                                                        >  1:3017438/1‑150 (MQ=255)

                                                                                                aTGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAt                                                      >  2:2768808/1‑150 (MQ=255)

                                                                                                   cACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTa                                                   <  2:4722644/150‑1 (MQ=255)

                                                                                                     cgggcgggCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACa                                                 >  2:3723971/1‑150 (MQ=255)

                                                                                                      gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                                >  2:2844871/1‑150 (MQ=255)

                                                                                                      gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGGGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                                >  2:204804/1‑150 (MQ=255)

                                                                                                       ggcgggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGc                                               >  2:3549269/1‑150 (MQ=255)

                                                                                                         cgggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCt                                             <  1:3939156/150‑1 (MQ=255)

                                                                                                          gggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTg                                            >  1:1015541/1‑150 (MQ=255)

                                                                                                             cgcAGCGGTGTTTCCAGGCCTCGCGCAACTCACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatt                                         <  1:1800501/150‑1 (MQ=255)

                                                                                                              gcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatta                                        <  1:3604374/150‑1 (MQ=255)

                                                                                                               cAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCCAAGTTATTAAGCTGATTTACAGCCTGattat                                       >  1:4532053/1‑150 (MQ=255)

                                                                                                                 gCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATtg                                     <  1:4473120/150‑1 (MQ=255)

                                                                                                                 gCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATtg                                     <  1:762942/150‑1 (MQ=255)

                                                                                                                     tgtTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGATTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCg                                 <  2:1645649/150‑1 (MQ=255)

                                                                                                                     tgtTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCg                                 <  2:3125816/150‑1 (MQ=255)

                                                                                                                     tgtTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCg                                 >  1:2899053/1‑150 (MQ=255)

                                                                                                                       tttCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTg                               <  2:666181/150‑1 (MQ=255)

                                                                                                                         tCCAGTCCTCGCGCAACTAACGCTTCATTACCCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAg                             <  1:3516233/150‑1 (MQ=255)

                                                                                                                          ccAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAgg                            <  2:4282725/150‑1 (MQ=255)

                                                                                                                           cAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggc                           <  1:439669/150‑1 (MQ=255)

                                                                                                                           cAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggc                           <  1:265908/150‑1 (MQ=255)

                                                                                                                               ccTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGgcggcg                       >  2:3486999/1‑150 (MQ=255)

                                                                                                                                  cgcgCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGcac                    <  1:1078297/150‑1 (MQ=255)

                                                                                                                                    cgcAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACAt                  >  2:3977243/1‑150 (MQ=255)

                                                                                                                                    cgcAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACAt                  >  1:1643900/1‑150 (MQ=255)

                                                                                                                                     gcAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATa                 <  2:3473927/150‑1 (MQ=255)

                                                                                                                                      cAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATaa                >  1:3378761/1‑150 (MQ=255)

                                                                                                                                       aaCTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAAt               <  2:2481162/150‑1 (MQ=255)

                                                                                                                                       aaCTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAAt               >  1:4664924/1‑150 (MQ=255)

                                                                                                                                       aaCTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCAGCAGGTGTGACGCCTCCGCCCAGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAAt               >  1:2364830/1‑150 (MQ=255)

                                                                                                                                        aCTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATc              <  1:819317/150‑1 (MQ=255)

                                                                                                                                        aCTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATc              >  1:2101543/1‑150 (MQ=255)

                                                                                                                                           aaCGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTACGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAg           >  2:3050779/1‑150 (MQ=255)

                                                                                                                                           aaCGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAg           <  2:1463815/150‑1 (MQ=255)

                                                                                                                                               cTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAtt       >  1:4058382/1‑150 (MQ=255)

                                                                                                                                                  cATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTcac    <  1:5107100/150‑1 (MQ=255)

                                                                                                                                                  cATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTcac    <  2:4950953/150‑1 (MQ=255)

CATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCACAG  >  CP009273/2236949‑2237246

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