BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,709,092	G→A	5.5%	G378S (GGC→AGC)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,709,092	0	G	A	5.5%	569.9 / 14.6	181	G378S (GGC→AGC)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (73/98); new base A (7/3); total (80/101)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 182 total aligned reads displayed.

TCGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGT  >  NC_002947/5708943‑5709239

tcGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTgg                                                                                                                                                     <  2:4912148/150‑1 (MQ=255)

  gACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGc                                                                                                                                                   >  1:13467/1‑150 (MQ=255)

      tGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATg                                                                                                                                               >  1:4211108/1‑150 (MQ=255)

       gCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCATACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGt                                                                                                                                              <  2:1456764/150‑1 (MQ=255)

       gCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGt                                                                                                                                              <  1:565337/150‑1 (MQ=255)

          ttCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGa                                                                                                                                           >  2:3460776/1‑150 (MQ=255)

          ttCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGa                                                                                                                                           <  2:2718997/150‑1 (MQ=255)

           tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  1:5349035/150‑1 (MQ=255)

           tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  2:2155625/150‑1 (MQ=255)

           tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  1:4739243/150‑1 (MQ=255)

           tCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGaa                                                                                                                                          <  1:4663171/150‑1 (MQ=255)

            cGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAg                                                                                                                                         <  2:1176579/150‑1 (MQ=255)

              cAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGat                                                                                                                                       >  1:3837394/1‑150 (MQ=255)

                 gaacgaacTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATc                                                                                                                                    >  2:2474892/1‑150 (MQ=255)

                    cgaacTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCGCACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGt                                                                                                                                 <  1:815562/150‑1 (MQ=255)

                        cTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGAtt                                                                                                                             <  1:3472219/150‑1 (MQ=255)

                         tGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTa                                                                                                                            <  1:203737/150‑1 (MQ=255)

                            gcCCGCAGTCAGCCCACCGCGAGGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgc                                                                                                                         <  2:5330055/150‑1 (MQ=255)

                              ccGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccg                                                                                                                       >  2:3502222/1‑150 (MQ=255)

                              ccGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccg                                                                                                                       >  2:1610526/1‑150 (MQ=255)

                              ccGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccg                                                                                                                       >  1:1204740/1‑150 (MQ=255)

                                    tCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccgctcgcc                                                                                                                 >  1:5225567/1‑150 (MQ=255)

                                       gCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTAcgccgctcgccgct                                                                                                              >  2:3935694/1‑150 (MQ=255)

                                        cccACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCtt                                                                                                             <  1:1649404/150‑1 (MQ=255)

                                        cccACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCtt                                                                                                             <  1:2989599/150‑1 (MQ=255)

                                           aCCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCg                                                                                                          >  1:4176380/1‑150 (MQ=255)

                                           aCCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCg                                                                                                          >  1:857680/1‑150 (MQ=255)

                                            ccGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGt                                                                                                         >  1:4614361/1‑150 (MQ=255)

                                            ccGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGt                                                                                                         <  1:2894435/150‑1 (MQ=255)

                                            ccGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGt                                                                                                         <  1:695545/150‑1 (MQ=255)

                                             cgcgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTg                                                                                                        <  1:1020546/150‑1 (MQ=255)

                                              gcgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTga                                                                                                       <  1:1009466/150‑1 (MQ=255)

                                              gcgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTga                                                                                                       <  1:2971500/150‑1 (MQ=255)

                                               cgATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTgag                                                                                                      >  2:1130784/1‑150 (MQ=255)

                                                   ggCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGAccc                                                                                                  <  1:3428089/150‑1 (MQ=255)

                                                     ctgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTg                                                                                                >  1:3380727/1‑150 (MQ=255)

                                                      tgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGc                                                                                               >  1:2700360/1‑150 (MQ=255)

                                                      tgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGc                                                                                               >  2:2402009/1‑150 (MQ=255)

                                                      tgctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGc                                                                                               >  1:2346149/1‑150 (MQ=255)

                                                       gctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATcc                                                                                                                                   <  2:4298806/113‑1 (MQ=255)

                                                       gctCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCa                                                                                              <  1:4760576/150‑1 (MQ=255)

                                                         tCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTgag                                                                                                      >  1:1765493/1‑140 (MQ=255)

                                                         tCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGa                                                                                            >  2:3576820/1‑150 (MQ=255)

                                                            aTGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCa                                                                                         <  2:1003735/150‑1 (MQ=255)

                                                            aTGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCa                                                                                         <  1:3283434/150‑1 (MQ=255)

                                                               ccAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCt                                                                                      >  2:4724652/1‑150 (MQ=255)

                                                                cAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTa                                                                                     >  1:1103688/1‑150 (MQ=255)

                                                                cAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTa                                                                                     >  1:3383633/1‑150 (MQ=255)

                                                                  gcgcTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATc                                                                                   <  2:3028116/150‑1 (MQ=255)

                                                                   cgcTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATcc                                                                                  <  1:2940661/150‑1 (MQ=255)

                                                                    gcTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCa                                                                                 <  2:611543/150‑1 (MQ=255)

                                                                            ccGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGc                                                                         <  2:5410990/150‑1 (MQ=255)

                                                                             cGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                        <  1:3778546/150‑1 (MQ=255)

                                                                             cGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                        <  2:4806053/150‑1 (MQ=255)

                                                                             cGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCa                                                                        <  2:947557/150‑1 (MQ=255)

                                                                                           cacGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGaa                                                          <  1:142478/150‑1 (MQ=255)

                                                                                                 caggaatgATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgt                                                    <  2:2892740/143‑1 (MQ=255)

                                                                                                  cTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                                   <  2:4176006/150‑1 (MQ=255)

                                                                                                  cTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAtgtg                                                   <  2:3687673/150‑1 (MQ=255)

                                                                                                     cTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGcac                                                >  2:321117/1‑150 (MQ=255)

                                                                                                      tCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGcaca                                               >  2:5386841/1‑150 (MQ=255)

                                                                                                         aTCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATcc                                            <  2:1774468/150‑1 (MQ=255)

                                                                                                         aTCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATcc                                            <  2:4616884/150‑1 (MQ=255)

                                                                                                          tCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCg                                           >  1:4644973/1‑150 (MQ=255)

                                                                                                            tGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCa                                         >  2:5259633/1‑150 (MQ=255)

                                                                                                            tGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCa                                         >  1:3644938/1‑150 (MQ=255)

                                                                                                             ggcggcGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCAt                                        <  2:1271292/150‑1 (MQ=255)

                                                                                                                 gcgATCCAAGTTCTGTAGAATGTCTGCAAATCTCTGGTCTTTGTCGAAGATATCCGTTATTACGAAGCTCGACTCTTGATTTAGACCCTGCAGAAAATCTATCCAGACAAGTCTGTACACGACCCTAACGATGTGCACATCCGCAtttaa                                    <  2:2803198/150‑3 (MQ=255)

                                                                                                                  cgAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTccc                                   <  1:5259633/150‑1 (MQ=255)

                                                                                                                   gAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCa                                  <  1:321117/150‑1 (MQ=255)

                                                                                                                    aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 <  1:2474892/150‑1 (MQ=255)

                                                                                                                    aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 <  1:1881292/150‑1 (MQ=255)

                                                                                                                    aaCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 >  1:3186742/1‑150 (MQ=255)

                                                                                                                     aCCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAg                                >  2:353997/1‑150 (MQ=255)

                                                                                                                        aaGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCca                             >  2:2065967/1‑150 (MQ=255)

                                                                                                                           gtcggtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCcatca                          >  2:1024305/1‑150 (MQ=255)

                                                                                                                             cggtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCGatcatc                        <  2:4177068/150‑1 (MQ=255)

                                                                                                                              ggtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCg                       >  2:2173843/1‑150 (MQ=255)

                                                                                                                              ggtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCg                       >  2:4944967/1‑150 (MQ=255)

                                                                                                                               gtcgAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCgg                      >  1:2189821/1‑150 (MQ=255)

                                                                                                                                  gAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCcat                            <  2:1690025/141‑1 (MQ=255)

                                                                                                                                    aGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGTACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCg                 <  1:3697276/150‑1 (MQ=255)

                                                                                                                                       gCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGca              >  2:2441756/1‑150 (MQ=255)

                                                                                                                                       cttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAt                                                      >  1:2110580/6‑110 (MQ=255)

                                                                                                                                       cttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGAt                                                      <  2:2110580/105‑1 (MQ=255)

                                                                                                                                        ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCaa                                 >  1:5175336/5‑130 (MQ=255)

                                                                                                                                        ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcac             >  2:3417256/5‑150 (MQ=255)

                                                                                                                                        ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcac             >  2:3182358/5‑150 (MQ=255)

                                                                                                                                        ttcgacaTCGCCCAGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcac             >  1:2724521/5‑150 (MQ=255)

                                                                                                                                         tGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcaca            >  2:4614806/1‑150 (MQ=255)

                                                                                                                                           cacaTCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcacaca          >  1:2151816/1‑150 (MQ=255)

                                                                                                                                              aTCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAgg       >  2:3944856/1‑150 (MQ=255)

                                                                                                                                               tCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCATGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCAACGGCGCAGGCACACAAggg      >  1:1472825/1‑150 (MQ=255)

                                                                                                                                                 gCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcacac           >  1:160747/1‑143 (MQ=255)

                                                                                                                                                 gCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGcacac           <  2:160747/143‑1 (MQ=255)

                                                                                                                                                 gCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCa    <  2:90230/150‑1 (MQ=255)

                                                                                                                                                  cTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAg   <  2:664313/150‑1 (MQ=255)

                                                                                                                                                   tGGGCGATGTCGAAGATATCCGTGTTTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGt  <  2:1735235/150‑1 (MQ=255)

                                                                                                                                                   tGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGt  <  1:2075979/150‑1 (MQ=255)

                                                                                                                                                   tGGGCGATGTCGAAGATATCCGTGATTACGCCGCACGCCGCTTGCGTGAGACCCTCCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACCATGTGCCCATCCGCATTTCCCAAGCCATCAACGCCGCCGCCACACAAGGGCAga  >  2:5190462/1‑149 (MQ=255)

TCGACCTGCCTTCGCAGGAACGAACTGCGCCCGCTGTCAGCCCACCGCGATGGCTGCTCAATGCCAGCGCTGAGGACCGCTACGCCTATCACACGGCACTGCTCGATCTGGCGGCGAACCAAGGTCGGTCGAAGGGCTGCACATCGCTGGGCGATGTCGAAGATATCCGTGATTACGCCGCTCGCCGCTTGCGTGAGACCCTGCAGACCAGCTATCCAGACAAGGCAGTACACGACCCGAACGATGTGCACATCCGCATTTCCCAAGCCATCATCGGCGCCGGCACACAAGGGCAGT  >  NC_002947/5708943‑5709239

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