BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,445,466	T→A	6.1%	F6Y (TTT→TAT) ‡	PP_2140 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,445,466	0	T	A	6.1%	585.3 / 19.2	174	F6Y (TTT→TAT) ‡	PP_2140	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (78/84); new base A (4/7); total (82/92)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 178 total aligned reads displayed.

GGTGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACAAAGT  >  NC_002947/2445317‑2445614

ggTGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCtt                                                                                                                                                       <  2:572432/150‑1 (MQ=255)

   gATCCGCTACAGCCGCAAGACCACGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGc                                                                                                                                                    <  2:2278775/150‑1 (MQ=255)

   gATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGc                                                                                                                                                    <  2:1890724/150‑1 (MQ=255)

     tCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCgg                                                                                                                                                  <  2:2337697/150‑1 (MQ=255)

     tCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCgg                                                                                                                                                  <  2:101473/150‑1 (MQ=255)

        gCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGccc                                                                                                                                               >  1:4258790/1‑150 (MQ=255)

        gCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGccc                                                                                                                                               <  2:2254737/150‑1 (MQ=255)

         cTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGcccc                                                                                                                                              <  1:3234843/150‑1 (MQ=255)

          tACAGCCGCAAGACCAAGGAGCGGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑a                                                                                                                                            >  1:3838104/1‑150 (MQ=255)

           aCAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCAAAGCGGCCCC‑aa                                                                                                                                           <  1:2984791/150‑1 (MQ=255)

           aCAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCAAAGCGGCCCC‑aa                                                                                                                                           <  1:354421/150‑1 (MQ=255)

            cAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGTGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCAAAGCGGCCCC‑aaa                                                                                                                                          <  2:1414103/150‑2 (MQ=255)

            cAGCCGCAAGACCAACGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCAGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AAt                                                                                                                                          <  2:3703002/150‑1 (MQ=255)

             aGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCAAAGCGGCCCC‑aaat                                                                                                                                         <  1:4385116/150‑3 (MQ=255)

              gCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCa                                                                                                                                        >  2:2336843/1‑150 (MQ=255)

               ccGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCAt                                                                                                                                       >  2:123198/1‑150 (MQ=255)

                 gCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCAttt                                                                                                                                     >  1:3486590/1‑150 (MQ=255)

                  cAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTa                                                                                                                                    <  2:443160/150‑1 (MQ=255)

                  cAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTa                                                                                                                                    <  2:945937/150‑1 (MQ=255)

                  cAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTa                                                                                                                                    <  1:3584881/150‑1 (MQ=255)

                        cAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGc                                                                                                                              >  2:1743175/1‑150 (MQ=255)

                          aGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAg                                                                                                                            <  2:4318069/150‑1 (MQ=255)

                            gAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGat                                                                                                                          <  1:2661483/150‑1 (MQ=255)

                            gAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGat                                                                                                                          <  1:1718428/150‑1 (MQ=255)

                                aGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATcc                                                                                                                      <  1:1592326/150‑1 (MQ=255)

                                 gTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCg                                                                                                                     >  2:150429/1‑150 (MQ=255)

                                 gTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCg                                                                                                                     >  1:3465071/1‑150 (MQ=255)

                                 gTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTAGCGGCCCC‑AATCATTTATGATGCAGATATCCg                                                                                                                     <  2:2476020/150‑1 (MQ=255)

                                    tgtgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCt                                                                                                                  <  1:3205974/150‑1 (MQ=255)

                                      tgCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTg                                                                                                                >  1:2831668/1‑150 (MQ=255)

                                       gCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGt                                                                                                               >  2:1668365/1‑150 (MQ=255)

                                         aGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAg                                                                                                             <  1:4187521/150‑1 (MQ=255)

                                          gTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGc                                                                                                            >  2:3291963/1‑150 (MQ=255)

                                               aGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTAt                                                                                                       <  1:2196532/150‑1 (MQ=255)

                                                 gTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATgg                                                                                                     >  1:2110392/1‑150 (MQ=255)

                                                  tCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATggg                                                                                                    >  1:4117067/1‑150 (MQ=255)

                                                   cGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGa                                                                                                   <  2:4144254/150‑1 (MQ=255)

                                                    gATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAg                                                                                                  <  2:2201777/150‑1 (MQ=255)

                                                     aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAgg                                                                                                 <  2:4011633/150‑1 (MQ=255)

                                                     aTGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAgg                                                                                                 <  1:1688177/150‑1 (MQ=255)

                                                      tGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGTGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAggg                                                                                                >  1:2953998/1‑150 (MQ=255)

                                                      tGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAggg                                                                                                <  1:229391/150‑1 (MQ=255)

                                                          aaGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCGGATATCCGCCTTGTAGCCTTATGGGAGGGCACt                                                                                            >  2:39064/1‑150 (MQ=255)

                                                           aGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGAGCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTg                                                                                           <  2:402779/150‑1 (MQ=255)

                                                             ccGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGaa                                                                                         <  2:4082908/150‑1 (MQ=255)

                                                               gACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAAt                                                                                       <  1:2155551/150‑1 (MQ=255)

                                                                 gcGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATgg                                                                                     >  1:2042535/2‑150 (MQ=255)

                                                                 ccGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATgg                                                                                     >  1:2042541/1‑150 (MQ=255)

                                                                  cGGCTGGAAGGCATTCTACGACGGTAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGc                                                                                    >  2:986552/1‑150 (MQ=255)

                                                                     cTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                                 >  1:4436687/1‑150 (MQ=255)

                                                                     cTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                                 >  1:4436634/1‑150 (MQ=255)

                                                                     cTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                                 <  1:1347126/150‑1 (MQ=255)

                                                                     cTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                                 <  1:2903269/150‑1 (MQ=255)

                                                                     cTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                                 <  1:3403206/150‑1 (MQ=255)

                                                                     cTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCa                                                                                 <  2:3272427/150‑1 (MQ=255)

                                                                      tGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAg                                                                                <  2:3663866/150‑1 (MQ=255)

                                                                        gAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGa                                                                              >  1:3060740/1‑150 (MQ=255)

                                                                              aTTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTa                                                                        <  2:74737/150‑1 (MQ=255)

                                                                              aTTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTa                                                                        <  2:1925033/150‑1 (MQ=255)

                                                                               ttCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTac                                                                       >  1:636352/1‑150 (MQ=255)

                                                                                    cgacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCg                                                                  <  1:1295333/150‑1 (MQ=255)

                                                                                     gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGt                                                                 >  2:2009657/1‑150 (MQ=255)

                                                                                     gacgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGt                                                                 >  2:2009662/1‑150 (MQ=255)

                                                                                      acgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTa                                                                >  1:2894120/1‑150 (MQ=255)

                                                                                      acgGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTa                                                                <  2:2634920/150‑1 (MQ=255)

                                                                                           aaTGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTaccaac                                                           <  1:1570412/150‑1 (MQ=255)

                                                                                            aTGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTaccaacc                                                          >  1:542246/1‑150 (MQ=255)

                                                                                            aTGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTaccaacc                                                          >  1:3999227/1‑150 (MQ=255)

                                                                                              gCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAg                                                        >  2:2666038/1‑150 (MQ=255)

                                                                                               ccTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGa                                                       <  2:3548053/150‑1 (MQ=255)

                                                                                               ccTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGa                                                       <  1:267295/150‑1 (MQ=255)

                                                                                                 tGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGaaa                                                     <  1:1038778/150‑1 (MQ=255)

                                                                                                 tGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGaaa                                                     >  1:2421629/1‑150 (MQ=255)

                                                                                                      ggTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGtt                                                <  2:4497900/150‑1 (MQ=255)

                                                                                                             gATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGcc                                         >  2:2560175/1‑150 (MQ=255)

                                                                                                                   cgcTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTc                                   <  1:102696/150‑1 (MQ=255)

                                                                                                                      tGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGcc                                <  2:4258790/150‑1 (MQ=255)

                                                                                                                        aTCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCt                              <  2:1092149/150‑1 (MQ=255)

                                                                                                                         tCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                             <  2:2371259/150‑1 (MQ=255)

                                                                                                                         tCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTc                             <  2:573133/150‑1 (MQ=255)

                                                                                                                            aTTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCCTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAg                          >  1:798859/1‑150 (MQ=255)

                                                                                                                             ttGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGt                         <  1:922270/150‑1 (MQ=255)

                                                                                                                                gTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTc                                   <  2:3761080/137‑1 (MQ=255)

                                                                                                                                gTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTc                                   >  1:3761080/1‑137 (MQ=255)

                                                                                                                                    tgtgATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGc                  >  2:4263317/1‑150 (MQ=255)

                                                                                                                                     gtgATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGcc                 <  2:2906420/150‑1 (MQ=255)

                                                                                                                                       gATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCaa               >  1:3718977/1‑150 (MQ=255)

                                                                                                                                        aTTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAg              >  2:1151433/1‑150 (MQ=255)

                                                                                                                                        aTTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACAGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAg              >  1:3894785/1‑150 (MQ=255)

                                                                                                                                         attGGGGCCGCAAAGCGGCCCCAAATCATTTATGATGCAGATCTCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAg              >  1:2271437/2‑150 (MQ=255)

                                                                                                                                         attGGGGCCGCAAAGCGGCCCCAAATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAg              >  1:521549/2‑150 (MQ=255)

                                                                                                                                            ggggCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTg          <  1:2506289/150‑1 (MQ=255)

                                                                                                                                             gggCCGCAAAGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGa         >  1:757065/1‑150 (MQ=255)

                                                                                                                                               gCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAAc       >  1:2641981/1‑150 (MQ=255)

                                                                                                                                               gCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAAc       >  1:385479/1‑150 (MQ=255)

                                                                                                                                                  gCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACaaa    <  2:3465071/150‑1 (MQ=255)

                                                                                                                                                  gCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACaaa    <  1:1743175/150‑1 (MQ=255)

                                                                                                                                                   cTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACAAAg   >  2:3835971/1‑150 (MQ=255)

                                                                                                                                                    tttGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACAAAGt  <  2:3104567/150‑1 (MQ=255)

                                                                                                                                                    tttGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACAAAGt  <  1:4086361/150‑1 (MQ=255)

GGTGATCCGCTACAGCCGCAAGACCAAGGAGCAGTATGTGCAGTCCGAGGTCGATGGCAAGCCGACCGGCTGGAAGGCATTCTACGACGGCAATGCCTGGAAGGTTGAAGATAAGCGCTGATCCATTGGTGCTGTGATTTGGGGCCGCTTTGCGGCCCC‑AATCATTTATGATGCAGATATCCGCCTTGTAGCCTTATGGGAGGGCACTGAAATGGCCCAGGAACTCTACACCCGTACCAACCAGAAACTGTTCTTCGCCGGCCTCGCCCTCGAGTCCATGGCCAAGGCTGAACAAAGT  >  NC_002947/2445317‑2445614

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