BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,237,098	G→T	90.2%	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	90.2%	594.5 / 52.7	201	H12N (CAT→AAT)	folE	GTP cyclohydrolase I
Reads supporting (aligned to +/- strand): ref base G (12/8); new base T (53/128); total (65/136)
Fisher's exact test for biased strand distribution p-value = 9.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 203 total aligned reads displayed.

CATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCAC  >  CP009273/2236949‑2237244

cATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATg                                                                                                                                                    <  2:3247091/150‑1 (MQ=255)

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

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

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

        tGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCGCGCGCAACTAACGCTTCATTAACCAggg                                                                                                                                            >  2:110685/1‑150 (MQ=255)

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

           catcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCg                                                                                                                                         <  1:3687704/150‑1 (MQ=255)

            atcaGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGc                                                                                                                                        >  2:2950623/1‑150 (MQ=255)

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

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

                 ccAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttctt                                                                                                                                   <  2:1244479/150‑1 (MQ=255)

                 ccAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGcttctt                                                                                                                                   <  2:3432931/150‑1 (MQ=255)

                  cAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCttt                                                                                                                                  <  1:1808715/150‑1 (MQ=255)

                   aGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTTCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTa                                                                                                                                 >  2:211153/1‑150 (MQ=255)

                         agaTTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGt                                                                                                                           >  1:3027321/1‑150 (MQ=255)

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

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

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

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

                               agcagcTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGc                                                                                                                     <  1:1372387/150‑1 (MQ=255)

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

                                     tGCATGATTTCGGTCATCTGACCAGCAATAAGGCTTTTGCGCTTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTAt                                                                                                               <  1:3101225/150‑1 (MQ=255)

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

                                            tttCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAGTTATGATTtct                                                                                                        >  1:4214571/1‑150 (MQ=255)

                                              tCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTcc                                                                                                      >  2:260431/1‑150 (MQ=255)

                                                  tCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCa                                                                                                  <  1:4561598/150‑1 (MQ=255)

                                                    atatGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCCTTTATGATTTCTCCTGCAgg                                                                                                >  1:4889554/1‑150 (MQ=255)

                                                    atatGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAgg                                                                                                <  2:1455696/150‑1 (MQ=255)

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

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

                                                      atGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgtg                                                                                              <  2:449171/150‑1 (MQ=255)

                                                      atGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGgtg                                                                                              <  2:1875341/150‑1 (MQ=255)

                                                       tGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCGTTTATGATTTCTCCTGCAGgtgt                                                                                             <  1:2128394/150‑1 (MQ=255)

                                                       tGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCGTTTATGATTTCTCCTGCAGgtgt                                                                                             <  1:2128225/150‑1 (MQ=255)

                                                          ccAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGAc                                                                                          <  1:4594120/150‑1 (MQ=255)

                                                           cAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACg                                                                                         <  2:1706690/150‑1 (MQ=255)

                                                           cAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACg                                                                                         <  2:2130502/150‑1 (MQ=255)

                                                                  aaGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTATTTACGGAGTGATTGCATTTATTATTTGTCCTGCAGGTGTGACCCCACCGc                                                                                  >  1:1836086/1‑150 (MQ=255)

                                                                     gCTTGTGCGCGTTTCGTTCTCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCt                                                                               <  1:3728258/150‑1 (MQ=255)

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

                                                                      cTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCACCGCCCTg                                                                              >  1:322423/1‑150 (MQ=255)

                                                                        tttGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTAAGGAGTGATTGCATTTATGATTTCTCCGGCCGGTGTGACGCCTCCGCCCAGCg                                                                            >  1:1558676/1‑150 (MQ=255)

                                                                         ttGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGt                                                                           <  1:1247185/150‑1 (MQ=255)

                                                                         ttGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGATTCATTAACCAGGGCCGCTTCTTTTCTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGt                                                                           <  1:2173727/150‑1 (MQ=255)

                                                                          tGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTg                                                                          <  2:4316338/150‑1 (MQ=255)

                                                                           gcgcGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTgg                                                                         >  2:3214119/1‑150 (MQ=255)

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

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

                                                                               gTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAGTTAAGATTTCTCCTGCAGGTGTGACGCCTCCGCCCAGCGTGGGGGc                                                                     >  1:779764/1‑150 (MQ=255)

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

                                                                                 ttCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCAGCGTGGGGGCaa                                                                   >  1:1177569/1‑150 (MQ=255)

                                                                                   cGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAg                                                                 >  1:3135818/1‑150 (MQ=255)

                                                                                    gTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGt                                                                <  1:4860473/150‑1 (MQ=255)

                                                                                      tATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGGGGGGGCAAAGtta                                                              >  2:3523494/1‑150 (MQ=255)

                                                                                       aTCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTATTTACTGAGTGATTGCATTTATGATTTCTCCAGCAGGTGTGAAGCCTCCGCCCGGCGTGGGGGCAAAGttat                                                             >  1:1124289/1‑150 (MQ=255)

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

                                                                                           aTTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCCGTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCAGCGTGGGGGCCCAGTTATTaat                                                         >  1:4181922/1‑149 (MQ=255)

                                                                                               cATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGa                                                     <  2:4609718/150‑1 (MQ=255)

                                                                                                 tGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGGTGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCACAGTTATTAAGCTGGtt                                                   >  2:4046627/1‑150 (MQ=255)

                                                                                                  gCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAttt                                                  <  2:1953019/150‑1 (MQ=255)

                                                                                                  gCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTTCTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGAttt                                                  >  2:2361021/1‑150 (MQ=255)

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

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

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

                                                                                                      gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                              <  2:1641875/150‑1 (MQ=255)

                                                                                                      gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                              <  2:2030897/150‑1 (MQ=255)

                                                                                                      gggcgggcGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAg                                              <  2:3027357/150‑1 (MQ=255)

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

                                                                                                       ggcgggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGc                                             <  1:2041852/150‑1 (MQ=255)

                                                                                                        gcgggcgCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGcc                                            <  2:1828056/150‑1 (MQ=255)

                                                                                                            gcgcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGat                                        <  1:3775558/150‑1 (MQ=255)

                                                                                                            gcgcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATGTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGat                                        <  2:37322/150‑1 (MQ=255)

                                                                                                             cgcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatt                                       <  1:3408859/150‑1 (MQ=255)

                                                                                                             cgcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatt                                       <  2:4626168/150‑1 (MQ=255)

                                                                                                             cgcAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGatt                                       <  1:3836510/150‑1 (MQ=255)

                                                                                                                aGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCGTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGattatt                                    >  1:3720638/1‑150 (MQ=255)

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

                                                                                                                   ggTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATtgtg                                 <  1:2482399/150‑1 (MQ=255)

                                                                                                                   ggTGTTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATtgtg                                 <  1:2373774/150‑1 (MQ=255)

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

                                                                                                                      gtTTCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCGTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGt                              <  2:694356/150‑1 (MQ=255)

                                                                                                                      gtTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCGTGTTTTCTGAGTGTTGGGCTGTATGATTTCTCCAGCAGGTGTGAAGCCGCCGCCCGGCGGGGGGGCAAAGTTATTAAGCTGATTTACAGCCAGATTAATGTGGGt                              >  1:1099089/1‑150 (MQ=255)

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

                                                                                                                        ttCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGa                            <  2:1585198/150‑1 (MQ=255)

                                                                                                                         tCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAg                           >  2:4313128/1‑150 (MQ=255)

                                                                                                                         tCCAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAg                           >  2:1275372/1‑150 (MQ=255)

                                                                                                                           cAGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggc                         >  2:358307/1‑150 (MQ=255)

                                                                                                                            aGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggcg                        <  2:3885645/150‑1 (MQ=255)

                                                                                                                            aGTCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggcg                        <  1:3417581/150‑1 (MQ=255)

                                                                                                                              tCCTCGCGCAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCAATTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAggcggc                      >  2:1137881/1‑150 (MQ=255)

                                                                                                                               ccTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCAGTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGgcggcg                     >  1:3805560/1‑150 (MQ=255)

                                                                                                                                   gcgcAACTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGcaca                 <  2:4748975/150‑1 (MQ=255)

                                                                                                                                       aaCTACCGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAAt             <  2:2222071/150‑1 (MQ=255)

                                                                                                                                       aaCTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAAt             <  2:3783967/150‑1 (MQ=255)

                                                                                                                                       aaCTAACGCTTCATGAACCAGGGCCGCTTCTTTAATGAGTGATGGCATTTATTATTCCTCCAGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTTTGCGTGAGGCGGCGCACATAAt             >  1:3386923/1‑150 (MQ=255)

                                                                                                                                         cTAACGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATcc           >  2:4210775/1‑150 (MQ=255)

                                                                                                                                            aCGCTTCATTAACCAGGGCCGCTTGTTTACTGAGTGATGGCATTTATGATTTATCCAGCAGGTGTGACGCCTCCGCCCAGCGTGGGGGCAAAGTTATTAAGCTGATTTAAAGCCTGATTATTTTTCGGGAGGCGGCGCACATAATACAGt        >  1:4984566/1‑150 (MQ=255)

                                                                                                                                            aCGCTTCATTAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGt        <  2:2776849/150‑1 (MQ=255)

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

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

CATCAAACTGTCATCAGCCAGGTCGAGATTCAGCAGCTGCATGATTTCGGTCATATGACCAGCAATAAGGCTTTTGCGCGTTTCGTTATCCATTTCATGCACGGGCGGGCGCAGCGGTGTTTCCAGTCCTCGCGCAACTAACGCTTCATGAACCAGGGCCGCTTCTTTACTGAGTGATGGCATTTATGATTTCTCCTGCAGGTGTGACGCCTCCGCCCTGCGTGGGGGCAAAGTTATTAAGCTGATTTACAGCCTGATTATTGTGCGTGAGGCGGCGCACATAATCCAGTATTCAC  >  CP009273/2236949‑2237244

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