BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,337,825	G→C	8.7%	G329A (GGG→GCG)	PP_2934 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,337,825	0	G	C	8.7%	383.5 / 22.2	126	G329A (GGG→GCG)	PP_2934	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (62/53); new base C (6/5); total (68/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 133 total aligned reads displayed.

TGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCC  >  NC_002947/3337676‑3337973

tGCCGGGCATTTTCTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTggg                                                                                                                                                       <  2:1472801/150‑1 (MQ=255)

tGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTggg                                                                                                                                                       >  1:2247552/1‑150 (MQ=255)

tGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTggg                                                                                                                                                       <  2:1833903/150‑1 (MQ=255)

 gCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTgggg                                                                                                                                                      <  1:379851/150‑1 (MQ=255)

   cGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTc                                                                                                                                                    <  2:306427/150‑1 (MQ=255)

    gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCg                                                                                                                                                   >  1:459767/1‑150 (MQ=255)

         tttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCC‑GGC‑GCTGGCGACCCCAgcg                                                                                                                                             >  1:971690/1‑150 (MQ=255)

          ttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGCAACGGCAGCGTTCTTCAGCCAGGC‑GCTggg                                                                                                                                                       <  1:2094694/140‑1 (MQ=255)

          ttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTggg                                                                                                                                                       >  2:2094694/1‑140 (MQ=255)

          ttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAgcg                                                                                                                                             <  2:291682/150‑1 (MQ=255)

             gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCC‑GGC‑GCTGGCGACCCCAGCGcctg                                                                                                                                         <  2:971690/150‑3 (MQ=255)

              tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑gg                                                                                                                                        >  2:1080938/1‑150 (MQ=255)

                      gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGt                                                                                                                                >  1:2897353/1‑150 (MQ=255)

                        ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTc                                                                                                                              >  2:2723758/1‑150 (MQ=255)

                        ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTc                                                                                                                              >  2:725479/1‑150 (MQ=255)

                        ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCC‑GGC‑GCTGGCGACCCCAGCGCCTGGCTGaa                                                                                                                                   <  2:2192739/145‑2 (MQ=38)

                         cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                             >  2:2579541/1‑150 (MQ=255)

                         cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                             >  2:2394618/1‑150 (MQ=255)

                          aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAg                                                                                                                            >  1:1475129/1‑150 (MQ=255)

                              cGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGc                                                                                                                        <  2:337927/150‑1 (MQ=255)

                              cGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGc                                                                                                                        <  2:1475129/150‑1 (MQ=255)

                                aTGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGa                                                                                                                      <  2:2529574/150‑1 (MQ=255)

                                 tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGAc                                                                                                                     <  1:632378/150‑1 (MQ=255)

                                 tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCCAc                                                                                                                     >  2:226635/1‑150 (MQ=255)

                                  gCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGcc                                                                                                                                                 >  1:1064035/1‑122 (MQ=255)

                                  gCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGcc                                                                                                                                                 <  2:1064035/122‑1 (MQ=255)

                                   ccGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACgg                                                                                                                   <  1:1080938/150‑1 (MQ=255)

                                    cGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                                  >  1:2759709/1‑150 (MQ=255)

                                           tACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCCTCAGCC‑GGC‑GCTGGCGACCCCAGCGcctggct                                                                                                                                      <  1:901913/123‑6 (MQ=37)

                                           tACGCGCGAATACCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                                  <  1:2095640/143‑1 (MQ=255)

                                              gcgcgAATGCCGGCGGTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcga                                                                                                        <  2:493716/150‑2 (MQ=255)

                                              gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTc                                                                                                                              <  1:255883/128‑1 (MQ=255)

                                              gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                        >  2:253321/1‑150 (MQ=255)

                                                gcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACgg                                                                                                                   <  2:1273751/137‑1 (MQ=255)

                                                gcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCAAGCGCC‑GGCTGAGCGTTCAGCTGCGACgg                                                                                                                   >  1:1273751/1‑137 (MQ=255)

                                                     tGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCgcgccgcgcccaggaa                                                                                                 <  2:976978/150‑7 (MQ=255)

                                                            gcTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTa                                                                                          <  1:2399442/150‑1 (MQ=255)

                                                              tgtgCAAGGGCGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCGGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                        <  1:2736443/150‑1 (MQ=255)

                                                              tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGc                                                                                                                                  >  1:1324179/1‑108 (MQ=255)

                                                              tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                        >  1:394400/1‑150 (MQ=255)

                                                              tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                        <  2:2787784/150‑1 (MQ=255)

                                                               gtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGc                                                                                       >  2:955778/1‑150 (MQ=255)

                                                                tgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGcc                                                                                      >  2:462858/1‑150 (MQ=255)

                                                                     ggACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgcgc                                                                                                               >  1:2498143/1‑120 (MQ=255)

                                                                           cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGTGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa                                                                           >  1:2037857/1‑150 (MQ=255)

                                                                           cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa                                                                           <  2:2044133/150‑1 (MQ=255)

                                                                               ggCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                               <  1:1838042/126‑1 (MQ=255)

                                                                               ggCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                               >  2:1838042/1‑126 (MQ=255)

                                                                                 cGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                     <  1:1009250/118‑1 (MQ=255)

                                                                                 cGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                     >  2:1009250/1‑118 (MQ=255)

                                                                                 cGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGa                                                                     >  2:1061894/1‑150 (MQ=255)

                                                                                  gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCGACAAAGTTGGGCACGCCCCGGGCAAGCCGAg                                                                    >  2:1426355/1‑150 (MQ=255)

                                                                                  gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                                    >  1:591600/1‑150 (MQ=255)

                                                                                      aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                                >  1:1053781/1‑150 (MQ=255)

                                                                                          ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAg                                                            >  1:2180285/1‑150 (MQ=255)

                                                                                                   ccATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagcca                                                                                                  <  1:270331/103‑1 (MQ=255)

                                                                                                   ccATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagcca                                                                                                  >  2:270331/1‑103 (MQ=255)

                                                                                                   ccATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCaga                                                   >  1:3002609/1‑150 (MQ=255)

                                                                                                   ccATCATTCACTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCaga                                                   <  2:1263540/150‑1 (MQ=255)

                                                                                                    catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCAcaca                                                     >  1:1487058/1‑144 (MQ=255)

                                                                                                    catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                                  >  1:787461/1‑150 (MQ=255)

                                                                                                    catcatTCGCTGCAGCGAGAAACGGCAGCGTCAAGCAG‑AAGGCGGCTTGGGTCGCCATCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                                  <  1:687015/150‑1 (MQ=255)

                                                                                                     atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                    <  2:2411732/99‑1 (MQ=255)

                                                                                                     atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                                    >  1:2411732/1‑99 (MQ=255)

                                                                                                       catTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGGCGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGAt                                               >  2:2400186/1‑150 (MQ=255)

                                                                                                             cTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGc                                         <  2:787461/150‑1 (MQ=255)

                                                                                                             cTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGc                                         >  1:1184309/1‑150 (MQ=255)

                                                                                                             cTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGc                                         <  1:955778/150‑1 (MQ=255)

                                                                                                               gCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCAt                                       >  1:2468304/1‑150 (MQ=255)

                                                                                                                 aGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAc                                     <  1:1188138/150‑1 (MQ=255)

                                                                                                                   cGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAccc                                   >  2:2530426/1‑150 (MQ=255)

                                                                                                                    tagaAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                                  >  1:893182/2‑150 (MQ=255)

                                                                                                                    gagaAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                                  >  2:2182687/1‑150 (MQ=255)

                                                                                                                     agaAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGACACGCAGAAGATCGAGGCATACCCCa                                 >  1:1946365/1‑150 (MQ=255)

                                                                                                                       aaaCGCCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAgg                               <  1:515022/150‑1 (MQ=255)

                                                                                                                        aaCGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGgt                              <  1:253321/150‑1 (MQ=255)

                                                                                                                          cGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCa                                                             <  1:1746735/117‑1 (MQ=255)

                                                                                                                          cGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCa                                                             >  2:1746735/1‑117 (MQ=255)

                                                                                                                           ggCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCTCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTa                           <  2:2882271/150‑1 (MQ=255)

                                                                                                                              aGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                                >  1:735127/1‑110 (MQ=255)

                                                                                                                              aGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                                <  2:735127/110‑1 (MQ=255)

                                                                                                                                 gTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGa                                            <  2:224612/127‑1 (MQ=255)

                                                                                                                                 gTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGa                                            >  1:224612/1‑127 (MQ=255)

                                                                                                                                  ttcttcAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTTAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCAc                    >  1:2845781/1‑150 (MQ=255)

                                                                                                                                  ttcttcAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCAc                    >  1:1896164/1‑150 (MQ=255)

                                                                                                                                   cgctcAGCC‑GGC‑GCTGGCGACCCCAGCGCC‑TGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAc                          <  1:1651681/139‑1 (MQ=255)

                                                                                                                                   cgctcAGCC‑GGC‑GCTGGCGACCCCAGCGCC‑TGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAc                          >  2:1651681/4‑142 (MQ=255)

                                                                                                                                    gctcAGCC‑GGC‑GCTGGCGACCCCAGCGCC‑TGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAgg                                                                        <  1:3041133/93‑1 (MQ=37)

                                                                                                                                    gctcAGCC‑GGC‑GCTGGCGACCCCAGCGCC‑TGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAgg                                                                        >  2:3041133/3‑95 (MQ=25)

                                                                                                                                     ttcAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgc                 >  2:1223027/1‑150 (MQ=255)

                                                                                                                                      tcAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgcg                >  2:2417644/1‑150 (MQ=255)

                                                                                                                                         gCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGccc             >  2:2469432/1‑150 (MQ=255)

                                                                                                                                         gCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGccc             <  2:1896164/150‑1 (MQ=255)

                                                                                                                                          ccAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                                    >  2:1675580/1‑94 (MQ=255)

                                                                                                                                          ccAGGC‑GCTGGGGTCGCCAGCCCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa            >  2:395121/1‑150 (MQ=255)

                                                                                                                                           cAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACACGGGACAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCAt           >  1:1136390/1‑150 (MQ=255)

                                                                                                                                                 gcTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGACCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCg      <  1:1825377/150‑1 (MQ=255)

                                                                                                                                                  cTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCgg     <  1:1426355/150‑1 (MQ=255)

                                                                                                                                                   tGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCggg    >  2:2037038/1‑150 (MQ=255)

                                                                                                                                                     gggTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGACTGCAGTGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACAGCCACGGCGCCCATCACCGGGcc  >  1:969297/1‑150 (MQ=255)

TGCCGGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGC‑GCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACCGGGCC  >  NC_002947/3337676‑3337973

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