BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,490,910	G→C	5.3%	P175R (CCG→CGG) ‡	clpV ←	protein ClpV1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,490,910	0	G	C	5.3%	370.2 / 9.2	113	P175R (CCG→CGG) ‡	clpV	protein ClpV1
Reads supporting (aligned to +/- strand): ref base G (57/50); new base C (4/2); total (61/52)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 113 total aligned reads displayed.

GTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGT  >  NC_002947/3490767‑3491057

gaagttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTgg                                                                                                                                               >  2:1530811/5‑150 (MQ=255)

gTACTTCGGACTCACGGCCCAGCACCGGGACGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGCCCGTGTACTGGTCCAGCGTCTGGGTCTTGCCGCTGGCCTTGCTAGACGTGCGCACTGGCGCCGCTGGGGGGGCCATCGGGGTgg                                                                                                                                               >  1:379732/1‑150 (MQ=255)

gTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTgg                                                                                                                                               <  2:2444468/150‑1 (MQ=255)

    ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTGGAGCACGCTTCGGCTGTTGTTGAGCGGGTACTGGGCCATCGCCGGGGCCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                                           >  1:57511/1‑150 (MQ=255)

    ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                                           >  2:2758463/1‑150 (MQ=255)

    ttCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctg                                                                                                                                           <  1:1405088/150‑1 (MQ=255)

       gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctgcct                                                                                                                                        <  2:3036965/150‑1 (MQ=255)

        aaCTCACGGACCAGCACGGGGTAGATGCGGGCTTGGGGGGTTGCGTGTGGAAGGGTGACCGTGGACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGcctgcctg                                                                                                                                       >  2:2339576/1‑150 (MQ=255)

           tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGAc                                                                                                                                    >  1:2270938/1‑150 (MQ=255)

           tCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGAc                                                                                                                                    >  1:499210/1‑150 (MQ=255)

             aCGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTc                                                                                                                                  >  1:2176094/1‑150 (MQ=255)

              cGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCg                                                                                                                                 >  2:1722351/1‑150 (MQ=255)

               ggCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCgg                                                                                                                                >  1:1107190/1‑150 (MQ=255)

                    aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTa                                                                                                                           <  1:997243/150‑1 (MQ=255)

                    aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCCCGCTGgccag                                                                                                                                           <  2:2667337/134‑3 (MQ=255)

                    aGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCCCGCTGgccag                                                                                                                                           >  1:2667337/1‑132 (MQ=255)

                      cACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGACCTGCCTGGACTCGGCGCTAcc                                                                                                                         >  1:237543/1‑150 (MQ=255)

                          gggTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTcc                                                                                                                     <  2:2622848/150‑1 (MQ=255)

                          gggTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTcc                                                                                                                     <  2:936110/150‑1 (MQ=255)

                                  gCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAgcgc                                                                                                             <  1:2714926/150‑1 (MQ=255)

                                   cGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGcgcg                                                                                                            <  1:2357483/150‑1 (MQ=255)

                                           gcgGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcct                                                                                                                                            <  2:1888208/110‑1 (MQ=255)

                                           gcgGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcct                                                                                                                                            >  1:1888208/1‑110 (MQ=255)

                                            cgGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGTTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTc                                                                                                   <  2:2687445/150‑1 (MQ=255)

                                              ggCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                                 <  2:1278786/150‑1 (MQ=255)

                                                   tttGGGTAAGGTTGACCGTGTACTGGTCTGTCTCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCa                                                                                            >  1:859614/1‑150 (MQ=255)

                                                   tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCa                                                                                            >  1:249691/1‑150 (MQ=255)

                                                    ttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAg                                                                                           <  2:2918087/150‑1 (MQ=255)

                                                      gggTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                                 <  2:3055547/142‑1 (MQ=255)

                                                      gggTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAg                                                                                                 >  1:3055547/1‑142 (MQ=255)

                                                        gtgaGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCGGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCa                                                                                       >  1:2923329/4‑150 (MQ=255)

                                                         tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCAt                                                                                      >  1:1920236/1‑150 (MQ=255)

                                                          aaGGTTGACCATGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATc                                                                                     <  1:110281/150‑1 (MQ=255)

                                                           aGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCg                                                                                    >  2:2591665/1‑150 (MQ=255)

                                                             gTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGc                                                                                  >  2:844859/1‑150 (MQ=255)

                                                             gTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGc                                                                                  >  1:351718/1‑150 (MQ=255)

                                                                   cGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAt                                                                            <  1:2758463/150‑1 (MQ=255)

                                                                     tgtACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTAGAttt                                                                          >  1:2358059/1‑150 (MQ=255)

                                                                       tACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGATCCCCGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                                        <  1:1722351/150‑1 (MQ=255)

                                                                        aCTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtttttt                                                                       <  2:729507/150‑1 (MQ=255)

                                                                         cTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGATACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                                      <  1:2380438/150‑1 (MQ=255)

                                                                         cTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGACCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTc                                                                      <  2:237543/150‑1 (MQ=255)

                                                                          tGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccagc                                                                                                                                          >  2:263858/1‑78 (MQ=38)

                                                                          tGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccagc                                                                                                                                          <  1:263858/81‑1 (MQ=38)

                                                                                aGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTcc                                                               >  2:801628/1‑150 (MQ=255)

                                                                                       gggTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtttt                                                                         <  1:2454814/133‑1 (MQ=255)

                                                                                       gggTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtttt                                                                         >  2:2454814/1‑133 (MQ=255)

                                                                                        ggTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGa                                                       >  2:2404757/1‑150 (MQ=255)

                                                                                             tGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGa                                                  <  2:2923329/150‑1 (MQ=255)

                                                                                              gCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGAt                                                 <  1:264703/150‑1 (MQ=255)

                                                                                                cGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCACTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGCTGATTTTTTCCATCTCCGGTGCCGAGGCGATca                                               <  2:3069731/150‑1 (MQ=255)

                                                                                                    cgcggttcTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATcaccac                                           <  2:1107190/143‑1 (MQ=255)

                                                                                                     gCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCAcg                                          >  2:2841702/1‑150 (MQ=255)

                                                                                                       cTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGATCCGGTGCCGAGGCGATCACCAcgcg                                        >  1:781619/1‑150 (MQ=255)

                                                                                                        ttGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTccac                                                                                                                   >  1:215492/1‑74 (MQ=255)

                                                                                                        ttGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTccac                                                                                                                   <  2:215492/74‑1 (MQ=255)

                                                                                                        ttGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGa                                       <  2:85916/150‑1 (MQ=255)

                                                                                                          gCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGcctgcctg                                                                                                                                       >  1:756311/1‑52 (MQ=255)

                                                                                                          gCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGcctgcctg                                                                                                                                       <  2:756311/52‑1 (MQ=255)

                                                                                                            tCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCa                                   <  2:57511/150‑1 (MQ=255)

                                                                                                            tCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCa                                   <  1:844859/150‑1 (MQ=255)

                                                                                                            tCGACGTGCTCACTGGCGCCGCTGAGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGTCGAGGCGATCACCACGCGACGCa                                   >  2:181215/1‑150 (MQ=255)

                                                                                                                 gTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCa                              >  1:1117702/1‑150 (MQ=255)

                                                                                                                   gCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGc                            >  2:2323286/1‑150 (MQ=255)

                                                                                                                    cTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTAGATTTTTTCCAGCTCCGGTGCCGAGGCGCTCACCACGCGACGCAGTTCAGCa                           <  2:2358059/150‑1 (MQ=255)

                                                                                                                     tCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAg                                                                                                       >  2:2353674/1‑73 (MQ=255)

                                                                                                                     tCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAg                                                                                                       <  1:2353674/73‑1 (MQ=255)

                                                                                                                     tCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCAt                          >  2:307569/1‑150 (MQ=255)

                                                                                                                      cACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATc                         >  1:1220065/1‑150 (MQ=255)

                                                                                                                           gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCaa                    <  2:1220065/150‑1 (MQ=255)

                                                                                                                           gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCaa                    <  2:2176094/150‑1 (MQ=255)

                                                                                                                           gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTACAGCTCCGGTGCCGAGGCGAACACCACGCGACGCAGTTCAGCATCGTCaa                    >  1:905393/1‑150 (MQ=255)

                                                                                                                           gcaccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGa                                                  <  2:1686447/117‑1 (MQ=255)

                                                                                                                             gccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                                        <  1:2696143/96‑1 (MQ=255)

                                                                                                                             gccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAttttt                                                                        >  2:2696143/1‑96 (MQ=255)

                                                                                                                             gccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGc                  >  2:2564972/1‑150 (MQ=255)

                                                                                                                              ccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCa                 <  1:1872072/150‑1 (MQ=255)

                                                                                                                               cgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCag                >  2:1863202/1‑150 (MQ=255)

                                                                                                                                gctgggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCaga               <  1:2742095/150‑1 (MQ=255)

                                                                                                                                   aggcaggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGcc            >  1:1712218/6‑150 (MQ=255)

                                                                                                                                    ggctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCa                              >  2:282763/1‑131 (MQ=255)

                                                                                                                                     gctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTg          <  1:1530811/150‑1 (MQ=255)

                                                                                                                                     gcaggCCAGCCCGCTGGCCAGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTg          >  2:358049/4‑150 (MQ=255)

                                                                                                                                      ctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCCCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGAAGAGCCTgc         <  1:2999680/150‑1 (MQ=255)

                                                                                                                                       tggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgca        >  2:2977754/1‑150 (MQ=255)

                                                                                                                                       tggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGCCTAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgca        <  1:2339576/150‑1 (MQ=255)

                                                                                                                                       aggCCAGCCCGCTGGCCAGCCTGGACTCGGCGCTACCCTCCACCAGCACCGCCAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgca        >  1:78531/2‑150 (MQ=255)

                                                                                                                                        ggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgcag       >  1:2920140/1‑150 (MQ=255)

                                                                                                                                         gCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATGAAAGCGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgcagc      >  1:28931/1‑150 (MQ=255)

                                                                                                                                         gCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgcagc      <  2:2920140/150‑1 (MQ=255)

                                                                                                                                         gCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgcagc      >  2:2843755/1‑150 (MQ=255)

                                                                                                                                         gCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTAATTCCAGCTCCTGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTGAAGCAGAGCCTgcagg      >  2:875745/1‑149 (MQ=255)

                                                                                                                                           cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATc                         >  2:3028883/1‑129 (MQ=255)

                                                                                                                                           cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATc                         <  1:3028883/129‑1 (MQ=255)

                                                                                                                                           cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGcagcag    <  2:499210/150‑1 (MQ=255)

                                                                                                                                           cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGcagcag    >  2:1855267/1‑150 (MQ=255)

                                                                                                                                           cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGcagcag    >  2:1507391/1‑150 (MQ=255)

                                                                                                                                            aGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGt                       >  1:1133290/1‑130 (MQ=255)

                                                                                                                                            aGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGt                       <  2:1133290/130‑1 (MQ=255)

                                                                                                                                             gCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGt  <  1:1507391/150‑1 (MQ=255)

GTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGGT  >  NC_002947/3490767‑3491057

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