BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,237,098	G→T	66.7%	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	66.7%	290.7 / inf	285	H12N (CAT→AAT)	folE	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand): ref base G (40/55); new base T (107/82); total (148/137)
Fisher's exact test for biased strand distribution p-value = 2.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

Only 100 of 285 total aligned reads displayed.

TCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCAC  >  CP009273/2236951‑2237244

tCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTaa                                                                                                                                                  <  1:4289704/150‑1 (MQ=255)

tCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGaa                                                                                                                                                  >  1:1393072/1‑150 (MQ=255)

   aaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCa                                                                                                                                               >  2:2098709/1‑150 (MQ=255)

   aaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCa                                                                                                                                               <  2:1601919/150‑1 (MQ=255)

   aaCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCa                                                                                                                                               >  2:3020547/1‑150 (MQ=255)

    aCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAg                                                                                                                                              <  1:1934491/150‑1 (MQ=255)

    aCTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAg                                                                                                                                              <  2:908345/150‑1 (MQ=255)

     cTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAgg                                                                                                                                             >  2:3071370/1‑150 (MQ=255)

     cTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAgg                                                                                                                                             <  2:3599713/150‑1 (MQ=255)

       gTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGc                                                                                                                                           >  2:2865411/1‑150 (MQ=255)

        tcatcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGcc                                                                                                                                          <  1:1879495/150‑1 (MQ=255)

         catcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCg                                                                                                                                         >  2:5419235/1‑150 (MQ=255)

         catcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCg                                                                                                                                         >  2:5177662/1‑150 (MQ=255)

         catcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCg                                                                                                                                         >  1:639670/1‑150 (MQ=255)

            caGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGctt                                                                                                                                      <  2:3289326/150‑1 (MQ=255)

            caGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGctt                                                                                                                                      >  2:246709/1‑150 (MQ=255)

             aGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttc                                                                                                                                     <  1:1227698/150‑1 (MQ=255)

                cAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCttt                                                                                                                                  >  1:3557241/1‑150 (MQ=255)

                  ggTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGGTCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTAc                                                                                                                                >  1:726968/1‑150 (MQ=255)

                   gTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACt                                                                                                                               >  1:3634629/1‑150 (MQ=255)

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

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

                        gaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTg                                                                                                                          <  1:1479899/150‑1 (MQ=255)

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

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

                         aTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGa                                                                                                                         >  2:3351609/1‑150 (MQ=255)

                           tCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATg                                                                                                                       >  2:2036228/1‑150 (MQ=255)

                           tCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATg                                                                                                                       >  2:1263123/1‑150 (MQ=255)

                             agcagcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGc                                                                                                                     >  1:692794/1‑150 (MQ=255)

                              gcagcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCa                                                                                                                    <  1:3069470/150‑1 (MQ=255)

                               cagcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAt                                                                                                                   >  1:5542032/1‑150 (MQ=255)

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

                                   tGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTAt                                                                                                               >  1:3482820/1‑150 (MQ=255)

                                    gCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATg                                                                                                              <  2:4167345/150‑1 (MQ=255)

                                    gCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATg                                                                                                              <  1:2875837/150‑1 (MQ=255)

                                       tGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAtt                                                                                                           >  2:2047996/1‑150 (MQ=255)

                                       tGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGAtt                                                                                                           <  2:3136345/150‑1 (MQ=255)

                                            tCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTcc                                                                                                      <  2:2596009/150‑1 (MQ=255)

                                              ggTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTg                                                                                                    >  1:607690/1‑150 (MQ=255)

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

                                                         cAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACg                                                                                         <  2:472015/150‑1 (MQ=255)

                                                            cAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCt                                                                                      >  2:3823060/1‑150 (MQ=255)

                                                             aaTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTc                                                                                     >  1:3023086/1‑150 (MQ=255)

                                                             aaTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTc                                                                                     >  2:2161542/1‑150 (MQ=255)

                                                             aaTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTc                                                                                     <  1:5443502/150‑1 (MQ=255)

                                                              aTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTcc                                                                                    >  2:1959916/1‑150 (MQ=255)

                                                              aTAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTcc                                                                                    <  1:3604659/150‑1 (MQ=255)

                                                                 aGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGcc                                                                                 >  2:3970605/1‑150 (MQ=255)

                                                                   gCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCt                                                                               <  2:3551347/150‑1 (MQ=255)

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

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

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

                                                                        tGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTg                                                                          >  1:4492470/1‑150 (MQ=255)

                                                                          cgcgTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTggg                                                                        <  1:2605636/150‑1 (MQ=255)

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

                                                                           gcgTTTCGTTATCCAGTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTgggg                                                                       <  2:5606379/150‑1 (MQ=255)

                                                                             gTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGc                                                                     >  1:4710342/1‑150 (MQ=255)

                                                                             gTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGc                                                                     <  2:5390387/150‑1 (MQ=255)

                                                                              tttCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCa                                                                    <  2:2031879/150‑1 (MQ=255)

                                                                                tCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCaaa                                                                  <  1:3861548/150‑1 (MQ=255)

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

                                                                                  gTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt                                                                <  1:219642/150‑1 (MQ=255)

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

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

                                                                                        cATTTCATGCACGGGCGGGCGAAGAGCTGTTTCCAGTCCTCGCGCAACTAAAGCATAATGAACAAGTGACGATTATTTACTGATAGATGGATTTAATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCAGCGTGCGGGCAAAGTTATaaa                                                          >  2:3919004/1‑150 (MQ=255)

                                                                                         aTTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAg                                                         >  2:444823/1‑150 (MQ=255)

                                                                                          tttCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGc                                                        <  2:726968/150‑1 (MQ=255)

                                                                                           ttCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCt                                                       <  1:1103441/150‑1 (MQ=255)

                                                                                            tcctGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTg                                                      <  2:4344746/147‑1 (MQ=255)

                                                                                            tCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTg                                                      >  1:5113626/1‑150 (MQ=255)

                                                                                              aTGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGACGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGGGTGACGCCTCCGCCCTGCGTGGGGGCCAAGTTATTAAGCTGAt                                                    >  2:1073449/1‑150 (MQ=255)

                                                                                               tGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAtt                                                   >  2:2494114/1‑150 (MQ=255)

                                                                                                gCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAttt                                                  >  1:3879035/1‑150 (MQ=255)

                                                                                                  aCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCACAGTTATTAAGCTGATTTAc                                                >  2:450520/1‑150 (MQ=255)

                                                                                                  aCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTAc                                                <  1:1787419/150‑1 (MQ=255)

                                                                                                  aCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTAc                                                <  2:4679775/150‑1 (MQ=255)

                                                                                                  aCGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTAc                                                >  1:748285/1‑150 (MQ=255)

                                                                                                   cgggcgggCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACa                                               <  2:1774901/150‑1 (MQ=255)

                                                                                                   cgggcgggCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACa                                               <  1:2036136/150‑1 (MQ=255)

                                                                                                   cgggcgggCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACa                                               <  2:1272986/150‑1 (MQ=255)

                                                                                                    gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                              >  2:889185/1‑150 (MQ=255)

                                                                                                    gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                              >  1:2719226/1‑150 (MQ=255)

                                                                                                    gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATGTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                              >  2:3737521/1‑150 (MQ=255)

                                                                                                           cgcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatt                                       >  1:2088129/1‑150 (MQ=255)

                                                                                                            gcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatta                                      <  1:1835250/150‑1 (MQ=255)

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

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

                                                                                                                    gtTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGt                              <  1:3073731/150‑1 (MQ=255)

                                                                                                                        ccAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAgg                          >  1:1148396/1‑150 (MQ=255)

                                                                                                                           gTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggcgg                       <  1:1034161/150‑1 (MQ=255)

                                                                                                                                       cTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATcc           <  2:1436568/150‑1 (MQ=255)

                                                                                                                                          aCGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGt        >  1:5063265/1‑150 (MQ=255)

                                                                                                                                            gCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAt      >  2:597072/1‑150 (MQ=255)

                                                                                                                                            gCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAt      <  1:4823024/150‑1 (MQ=255)

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

                                                                                                                                             cTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTAtt     <  2:1279987/150‑1 (MQ=255)

                                                                                                                                              ttCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTc    >  2:4922873/1‑150 (MQ=255)

                                                                                                                                               tCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTca   <  2:1229863/150‑1 (MQ=255)

                                                                                                                                               tCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTca   <  2:5069635/150‑1 (MQ=255)

                                                                                                                                                cATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTcac  >  1:4195199/1‑150 (MQ=255)

TCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCAC  >  CP009273/2236951‑2237244

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