BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,609,051	G→C	5.8%	P200A (CCT→GCT)	waaA ←	3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,609,051	0	G	C	5.8%	380.5 / 11.6	121	P200A (CCT→GCT)	waaA	3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand): ref base G (54/60); new base C (4/3); total (58/63)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 123 total aligned reads displayed.

GGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAGCGC  >  NC_002947/5608902‑5609198

ggCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTT‑TTCGTAGCCGGTGACTTGCACGCGCTCAgg                                                                                                                                                     <  1:64076/150‑1 (MQ=255)

 gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAggg                                                                                                                                                    >  2:409253/1‑150 (MQ=255)

  cATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGc                                                                                                                                                   <  2:777358/150‑1 (MQ=255)

    tCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                                 <  1:2360547/150‑1 (MQ=255)

    tCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                                 <  2:946933/150‑1 (MQ=255)

          ccGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCCCTGAGCGCGTGCaagt                                                                                                                                           <  2:1385651/150‑4 (MQ=255)

            gtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                         <  1:1236981/150‑1 (MQ=255)

            gtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                         <  2:424223/150‑1 (MQ=255)

            gtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                         <  2:1027896/150‑1 (MQ=255)

                 gggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCAc                                                                                                                                    <  2:1424144/150‑1 (MQ=255)

                       tGGNGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTAAC‑GTAAGAGCGT‑AGCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgc                                                                                                                              >  2:1520440/1‑150 (MQ=255)

                       tGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgc                                                                                                                              >  1:833049/1‑150 (MQ=255)

                       tGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgc                                                                                                                              >  1:801334/1‑150 (MQ=255)

                         gcggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTAAATGGTCAG‑GT‑AAAACCCGATCG‑AGCCTGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgcgc                                                                                                                            <  1:1488110/150‑1 (MQ=255)

                           ggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTg                                                                                                                          >  2:910025/1‑150 (MQ=255)

                               aTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGcc                                                                                                                      >  1:1063428/1‑150 (MQ=255)

                               aTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGGGGCAGCAGCTGTTCGTCGGTGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCCAGGAGACGGTGGGcc                                                                                                                      >  1:1993889/1‑150 (MQ=255)

                                  cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑TT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCt                                                                                                                     >  2:1657254/1‑148 (MQ=255)

                                  cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCt                                                                                                                     <  1:1657254/148‑1 (MQ=255)

                                  cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCg                                                                                                                   <  1:1470269/150‑1 (MQ=255)

                                  cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCg                                                                                                                   <  1:893157/150‑1 (MQ=255)

                                   acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGTTCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAGGCGCCGGGACTggg                                                                                                                  >  1:992199/1‑150 (MQ=255)

                                   acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTg                                                                                                                          <  2:245590/142‑1 (MQ=255)

                                   acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTg                                                                                                                          >  1:245590/1‑142 (MQ=255)

                                   acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCgg                                                                                                                  >  1:970004/1‑150 (MQ=255)

                                   acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCgg                                                                                                                  >  2:2307838/1‑150 (MQ=255)

                                     acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTgg                                                                                                                         <  1:581423/141‑1 (MQ=255)

                                     acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGtt                                                                                                                >  1:125333/1‑150 (MQ=255)

                                     acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGtt                                                                                                                <  2:125333/150‑1 (MQ=255)

                                        gggCGGTGGCCGCCCCCCAACTGTTCACGAAATGCCCTCGCGCGTGGGAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTGGATCG‑GGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                             <  2:207379/150‑1 (MQ=255)

                                        gggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGTGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                             <  1:1145558/150‑1 (MQ=255)

                                        gggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                             <  1:1647005/150‑1 (MQ=255)

                                          gcgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                           <  2:2270660/150‑1 (MQ=255)

                                            gcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                         <  2:1219711/150‑1 (MQ=255)

                                            gcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGCCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                         >  2:1598445/1‑150 (MQ=255)

                                               ggCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCa                                                                                                      >  1:828345/1‑150 (MQ=255)

                                                  cgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCCGCTGTTCGCCGATGGTCAG‑TT‑CGAACTTGATCT‑AGCCGGTGACTTGCACGCGCGCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACgg                                                                                                   >  2:1073653/1‑150 (MQ=255)

                                                   gccCCCCACCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGc                                                                                                  <  1:184937/150‑1 (MQ=255)

                                                    ccccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACGTGCACGCCCTGAGCGCGTGCaag                                                                                                                                            >  2:1768638/1‑105 (MQ=255)

                                                    ccccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACGTGCACGCCCTGAGCGCGTGCaag                                                                                                                                            <  1:1768638/107‑3 (MQ=255)

                                                     cccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGa                                                                                                <  2:208272/150‑1 (MQ=255)

                                                     cccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGa                                                                                                <  2:1826626/150‑1 (MQ=255)

                                                      ccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                               <  1:357032/150‑1 (MQ=255)

                                                      ccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                               >  1:1477548/1‑150 (MQ=255)

                                                       cccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATc                                                                                              <  1:1586838/150‑1 (MQ=255)

                                                          aaCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCTCCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                           <  1:662666/150‑1 (MQ=255)

                                                          aaCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                           <  1:2307838/150‑1 (MQ=255)

                                                          aaCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                           <  1:1733827/150‑1 (MQ=255)

                                                              gTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCa                                                                                       >  1:1158537/1‑150 (MQ=255)

                                                              gTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCa                                                                                       >  2:1476788/1‑150 (MQ=255)

                                                               ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAt                                                                                      >  1:2267500/1‑150 (MQ=255)

                                                                 cACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑AGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAtct                                                                                    <  1:1106434/150‑1 (MQ=255)

                                                                 cACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCCCAGATGCCGGGAGG‑GGAAGTCG‑GGG‑GCTCGGGGACGGGCGCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTGCGGTCTGCACGGCGATCAGGCTCAtct                                                                                    >  1:647042/1‑150 (MQ=255)

                                                                  aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGaa                                                                                                                                >  1:923894/1‑105 (MQ=255)

                                                                  aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                                    >  2:1929957/1‑117 (MQ=255)

                                                                  aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                                    <  1:1929957/117‑1 (MQ=255)

                                                                    gCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                                 <  2:1776608/150‑1 (MQ=255)

                                                                        tGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGTTTTCg                                                                                                             >  1:1947970/1‑118 (MQ=255)

                                                                        tGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                             <  2:1947970/118‑1 (MQ=255)

                                                                           ccTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCCG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCCCATCTCGGCCAGCAt                                                                          <  1:1716177/150‑1 (MQ=255)

                                                                             tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                         >  1:1447275/1‑85 (MQ=255)

                                                                             tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                         <  2:1447275/85‑1 (MQ=255)

                                                                             tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                                 >  2:974564/1‑141 (MQ=255)

                                                                                  cgTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                                   <  2:2422115/150‑1 (MQ=255)

                                                                                     ggCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCCCATCTCGGCCAGCATCGGCCGGGTc                                                                <  1:1077064/150‑1 (MQ=255)

                                                                                           agcTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTg                                                          <  1:1476788/150‑1 (MQ=255)

                                                                                            gcTGTTCGGCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTgg                                                         >  2:90965/1‑150 (MQ=255)

                                                                                              tGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCa                                                       >  2:1285785/1‑150 (MQ=255)

                                                                                               gTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGcca                                                      >  2:15909/1‑148 (MQ=255)

                                                                                                  cgtcgATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgc                                                   <  2:970004/150‑1 (MQ=255)

                                                                                                   gtcgATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                                   <  2:395896/133‑1 (MQ=255)

                                                                                                   gtcgATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                                   >  1:395896/1‑133 (MQ=255)

                                                                                                    tcgATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                                 <  2:1674906/150‑1 (MQ=255)

                                                                                                    tcgATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgcgc                                                 >  2:785448/1‑150 (MQ=255)

                                                                                                     cgATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGcc                                                <  2:833049/150‑1 (MQ=255)

                                                                                                          gTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCCATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAg                                           <  1:661961/150‑1 (MQ=255)

                                                                                                                gT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCCAAGCGCCCATAGACACg                                      >  1:733731/1‑150 (MQ=255)

                                                                                                                gT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACg                                      >  1:330265/1‑150 (MQ=255)

                                                                                                                gT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACg                                      >  2:1968564/1‑150 (MQ=255)

                                                                                                                       cTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                                    <  1:472520/114‑1 (MQ=255)

                                                                                                                       cTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                                    >  2:472520/1‑114 (MQ=255)

                                                                                                                       cTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAACGCGCCCATAGCCACGGGCCGa                                <  2:956562/150‑1 (MQ=255)

                                                                                                                       cTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGa                                <  2:2267500/150‑1 (MQ=255)

                                                                                                                        ttGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAg                               <  1:915693/150‑1 (MQ=255)

                                                                                                                          gATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcg                             <  2:1063428/150‑1 (MQ=255)

                                                                                                                           aTCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcgc                            <  2:647042/150‑1 (MQ=255)

                                                                                                                           aTCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcgc                            >  1:1579031/1‑150 (MQ=255)

                                                                                                                              g‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGACGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCa                         <  1:1968564/150‑1 (MQ=255)

                                                                                                                              g‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCa                         <  2:717866/150‑1 (MQ=255)

                                                                                                                                          ttGCACGCGCGCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTGCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcgc              >  2:1675695/1‑150 (MQ=255)

                                                                                                                                             cACGCGCTCAGGGCGTGCACCCAGGTCTCGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCAtt           >  2:728276/1‑150 (MQ=255)

                                                                                                                                              aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg          <  2:1993889/150‑1 (MQ=255)

                                                                                                                                              aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg          >  1:1342150/1‑150 (MQ=255)

                                                                                                                                                  gcTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCa      >  2:2385958/1‑150 (MQ=255)

                                                                                                                                                    tCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGa                                <  1:1785870/122‑1 (MQ=255)

                                                                                                                                                    tCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGa                                >  2:1785870/1‑122 (MQ=255)

                                                                                                                                                    tCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAgc    >  1:1215640/1‑150 (MQ=255)

                                                                                                                                                    tCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAgc    >  2:1099043/1‑150 (MQ=255)

                                                                                                                                                      aGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAgcgc  >  2:831775/1‑150 (MQ=255)

                                                                                                                                                      aGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAgcgc  >  2:236187/1‑150 (MQ=255)

GGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAG‑GT‑CGAACTTGATCG‑AGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCAGCGC  >  NC_002947/5608902‑5609198

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