BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,438,253	A→T	6.0%	K457* (AAG→TAG) ‡	PP_3932 →	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,438,253	0	A	T	6.0%	395.1 / 11.2	127	K457* (AAG→TAG) ‡	PP_3932	diguanylate cyclase
Reads supporting (aligned to +/- strand): ref base A (57/62); new base T (3/5); total (60/67)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

CCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGTA  >  NC_002947/4438105‑4438402

ccgcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCaa                                                                                                                                                       <  1:1643738/150‑1 (MQ=255)

 cgcgccATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAg                                                                                                                                                      <  1:1368494/150‑1 (MQ=255)

 cgcgccATACAGGACCGGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAg                                                                                                                                                      <  1:229100/150‑1 (MQ=255)

      cATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCg                                                                                                                                                 >  2:255060/1‑150 (MQ=255)

        tACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCggg                                                                                                                                               >  2:1126895/1‑150 (MQ=255)

         aCAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGTCCAAGGCCCgggg                                                                                                                                              <  1:1014152/150‑1 (MQ=255)

         aCAGGACCTGCACGTGACACAGGCCGACCAGCGAATGCCGATTAAAACGTGCAGAAGTTTGGCAGTCTGCAAGTGAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGGGCCGACCAGGCGCTGTACCGGGCCAAGGCCCgggg                                                                                                                                              >  1:1863182/1‑150 (MQ=255)

           aGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggtaca                                                                                                                                            <  1:628326/150‑5 (MQ=255)

             gACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGc                                                                                                                                          >  2:360585/1‑150 (MQ=255)

               ccTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCaa                                                                                                                                        <  1:529937/150‑1 (MQ=255)

                cTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTCTACCGGGCCAAGGCCCGGGGCCGCAAc                                                                                                                                       <  1:562386/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCg                                                                                                                                                 >  2:417216/1‑134 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCg                                                                                                                                                 <  1:417216/134‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                 <  2:1035156/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                 <  1:2239857/150‑1 (MQ=255)

                      gTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTc                                                                                                                                 <  1:360585/150‑1 (MQ=255)

                         acgcaGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGaa                                                                                                                              <  2:377650/147‑1 (MQ=255)

                         acacAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGaa                                                                                                                              >  1:1536975/1‑150 (MQ=255)

                            caGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGAGACGCTGACATGCCGGTCACCAGCCCTGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAg                                                                                                                           <  2:1147303/150‑1 (MQ=255)

                             aGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAgg                                                                                                                          <  1:799577/150‑1 (MQ=255)

                                  gACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggtac                                                                                                                                             >  1:592172/1‑123 (MQ=255)

                                   aCCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAg                                                                                                                    <  1:800622/150‑1 (MQ=255)

                                     cAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGca                                                                                                                  >  2:847935/1‑150 (MQ=255)

                                     cAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGca                                                                                                                  >  2:2290803/1‑150 (MQ=255)

                                      aGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGcac                                                                                                                 <  2:2316209/150‑1 (MQ=255)

                                         gcCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAgc                                                                                                              <  1:1680152/150‑1 (MQ=255)

                                             gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggt                                                                                                                                               >  2:442774/1‑112 (MQ=255)

                                             gCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑Ct                                                                                                         <  1:1711639/150‑1 (MQ=255)

                                                gATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGt                                                                                                      <  2:1649542/150‑1 (MQ=255)

                                                gATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGt                                                                                                      <  1:2048080/150‑1 (MQ=255)

                                                 aTTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCTTGGCCCggtac                                                                                                                                             >  1:2321324/1‑108 (MQ=255)

                                                   tACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGAt                                                                                                   <  2:1609584/150‑1 (MQ=255)

                                                   tACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGAt                                                                                                   <  2:1704552/150‑1 (MQ=255)

                                                        tGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGaa                                                                                              <  2:710923/150‑1 (MQ=255)

                                                        tGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGaa                                                                                              <  2:339082/150‑1 (MQ=255)

                                                            tATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAg                                                                                          >  1:779010/1‑150 (MQ=255)

                                                              tCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGcc                                                                                        >  2:474505/1‑150 (MQ=255)

                                                                 gtgtGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGc                                                                                     >  2:806944/1‑150 (MQ=255)

                                                                 gtgtGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGc                                                                                     >  2:1977823/1‑150 (MQ=255)

                                                                   gtgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCca                                                                                   <  2:2105374/150‑1 (MQ=255)

                                                                   gtgGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCca                                                                                   <  2:1755451/150‑1 (MQ=255)

                                                                       cAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACA‑CGCGCTGGTGATCTGAACTAGCCGGCcagcag                                                                               <  1:62675/150‑1 (MQ=255)

                                                                        aGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAgg                                                                              >  1:2259596/1‑150 (MQ=255)

                                                                            tGCGACGCAGACATGCCGGCCCCCAGCCTCTACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTccc                                                                          >  1:1147083/1‑150 (MQ=255)

                                                                            tGCGACGCAGACATGCAGGCACATCGCATTGACGAAGTGCTGGACATTTCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTccc                                                                          >  2:25642/1‑150 (MQ=255)

                                                                              cGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCTTATGCTCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCAc                                                                        <  2:790386/150‑1 (MQ=255)

                                                                                aCGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGa                                                                      >  1:2269542/1‑150 (MQ=255)

                                                                                  gCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACa                                                                    >  2:1457107/1‑150 (MQ=255)

                                                                                      aCATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCtt                                                                >  2:1558231/1‑150 (MQ=255)

                                                                                       cATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                               <  2:2176589/150‑1 (MQ=255)

                                                                                        aTGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGa                                                                                                    >  2:1685713/1‑112 (MQ=255)

                                                                                        aTGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGa                                                                                                    <  1:1685713/112‑1 (MQ=255)

                                                                                            cGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGAt                                                          <  1:1558231/150‑1 (MQ=255)

                                                                                             ggCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATc                                                         >  1:570093/1‑150 (MQ=255)

                                                                                              gCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCa                                                        <  1:1477751/150‑1 (MQ=255)

                                                                                               cccccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAg                                                       >  1:1507439/1‑150 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACa                                                                                                                            <  1:156715/79‑1 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACa                                                                                                                            >  2:156715/1‑79 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcac                                                    >  2:387213/1‑150 (MQ=255)

                                                                                                  ccAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcac                                                    >  2:640602/1‑150 (MQ=255)

                                                                                                     gCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGcacaca                                                 <  1:847935/150‑1 (MQ=255)

                                                                                                           acgaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAATCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGc                                                      <  1:2111828/139‑1 (MQ=255)

                                                                                                           acgaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAATCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGc                                                      >  2:2111828/1‑139 (MQ=255)

                                                                                                             gaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAg                                                                                                                    >  1:2446700/1‑76 (MQ=255)

                                                                                                             gaACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCa                                         >  1:338368/1‑150 (MQ=255)

                                                                                                                ccgctgGCCAGTGCCGAGCAGGCGACGTCCCGCGAAAACACTTGAGATGGTAACTGCTTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCa                                      <  2:2231421/148‑1 (MQ=255)

                                                                                                                 tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                     <  2:1368021/150‑1 (MQ=255)

                                                                                                                 tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                     >  1:128773/1‑150 (MQ=255)

                                                                                                                 tgctgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAg                                     >  1:377481/1‑150 (MQ=255)

                                                                                                                    tgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGaa                                  <  1:1457107/150‑1 (MQ=255)

                                                                                                                    tgGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGaa                                  <  2:2448041/150‑1 (MQ=255)

                                                                                                                     ggCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGAACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAAc                                 >  1:1280784/1‑150 (MQ=255)

                                                                                                                      gCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACa                                <  1:1231526/150‑1 (MQ=255)

                                                                                                                      gCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACa                                >  1:1187068/1‑150 (MQ=255)

                                                                                                                       ccAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAg                               >  1:525285/1‑150 (MQ=255)

                                                                                                                        cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGc                                                                                                                         >  1:104247/1‑60 (MQ=255)

                                                                                                                        cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGc                                                                                                                         <  2:104247/60‑1 (MQ=255)

                                                                                                                        cAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGc                              >  1:302574/1‑150 (MQ=255)

                                                                                                                           tGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcg                           >  1:892267/1‑150 (MQ=255)

                                                                                                                           tGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcg                           >  1:794893/1‑150 (MQ=255)

                                                                                                                            gCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgc                          <  2:779010/150‑1 (MQ=255)

                                                                                                                            gCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgc                          <  1:25642/150‑1 (MQ=255)

                                                                                                                               gACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGa                       <  2:1351099/150‑1 (MQ=255)

                                                                                                                                aCCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGTGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGAc                      >  1:988145/1‑150 (MQ=255)

                                                                                                                                      cgcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAAccc                <  2:570093/150‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAAccc                <  2:2390257/149‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAAccc                >  1:2390257/1‑149 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCg               <  2:226451/150‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCg               <  2:1536975/150‑1 (MQ=255)

                                                                                                                                       gcTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCg               <  2:525285/150‑1 (MQ=255)

                                                                                                                                           tACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAc           <  1:1977823/150‑1 (MQ=255)

                                                                                                                                           tACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAc           >  2:2382142/1‑150 (MQ=255)

                                                                                                                                            aCCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgcg                         <  1:242164/135‑1 (MQ=255)

                                                                                                                                            aCCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCcgcg                         >  2:242164/1‑135 (MQ=255)

                                                                                                                                            aCCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAcc          >  2:474567/1‑150 (MQ=255)

                                                                                                                                             ccGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCAccc         >  2:1941166/1‑150 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                               <  1:2350891/90‑1 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTg                                                               >  2:2350891/1‑90 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGt   >  2:409981/1‑150 (MQ=255)

                                                                                                                                                   cAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGt   >  1:2134737/1‑150 (MQ=255)

CCGCGCCATACAGGACCTGCACGTGACACAGGCCGACCAGCGCCTGCAGATTACCGTGAGTATCGGTGTGGCAGGCTGCGACGCAGACATGCCGGCCCCCAGCCTCGACGAACTGCTGGCCAGTGCCGACCAGGCGCTGTACCGGGCCAAGGCCCGGGGCCGCAACCGCGTCGAACAGGCCGAAGCACAGCGC‑CTGGTGATCTGAACTAGCCGGCCAGCAGGTCCCACGACAGCTTGGCGATCAGCACACACAACAGCACCAGGAACAGCCCGCGGACAAACCCGGCACCCTTGCGTA  >  NC_002947/4438105‑4438402

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