BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,134,011	A→C	5.9%	N349H (AAC→CAC)	tolC →	agglutination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,134,011	0	A	C	5.9%	318.6 / 9.3	115	N349H (AAC→CAC)	tolC	agglutination protein
Reads supporting (aligned to +/- strand): ref base A (63/45); new base C (3/4); total (66/49)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 116 total aligned reads displayed.

CAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTACAAC  >  NC_002947/5133864‑5134160

cAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCCAGAC‑TTCGCAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAAc                                                                                                                                                       <  1:42296/150‑1 (MQ=255)

 aaGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACg                                                                                                                                                      >  2:172829/1‑150 (MQ=255)

     tGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACACGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCAt                                                                                                                                                  <  2:2044172/150‑1 (MQ=255)

      ctcgttctcactgtcgaggacgtccagcagggtgcgctggcccaggccgaactggtcctggtaggccgcacgcacacgcttggtggtttcggcgtattcacgcgcggtcggcagctgcttgctggcgttgttcaTGGCGTGGAACGCCAt                                                                                                                                                  >  1:1271915‑M1/135‑150 (MQ=255)

           ggAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGaacaa                                                                                                                                            <  1:1586120/150‑1 (MQ=255)

                tGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGCGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCa                                                                                                                                     >  1:337483/1‑150 (MQ=255)

                   aCTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagca                                                                                                                                  >  1:138865/1‑150 (MQ=255)

                   aCTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagca                                                                                                                                  >  1:140517/1‑150 (MQ=255)

                     tACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaag                                                                                                                                >  2:2339439/1‑150 (MQ=255)

                     tACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGAACGAAAACCTGAACCTGGCGTGGAACGCCATAAACAA‑C‑GCCagcaaa                                                                                                                                >  2:144143/1‑149 (MQ=255)

                       cAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaagca                                                                                                                              >  2:1073434/1‑150 (MQ=255)

                          ccTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCt                                                                                                                           >  1:2299769/1‑150 (MQ=255)

                           cTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTg                                                                                                                          <  1:504163/150‑1 (MQ=255)

                            tGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGc                                                                                                                         <  1:2339439/150‑1 (MQ=255)

                             gTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGcc                                                                                                                        <  2:652321/150‑1 (MQ=255)

                                  gcgGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACcg                                                                                                                   >  1:51288/1‑150 (MQ=255)

                                          aGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAATATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCa                                                                                                                                                   <  2:1380746/112‑1 (MQ=255)

                                          aGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCa                                                                                                                                                   >  1:1380746/1‑112 (MQ=255)

                                           gCGACAAGGCTCGCATACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAAt                                                                                                          <  2:1133549/150‑1 (MQ=255)

                                              aCAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACg                                                                                                       >  1:1640617/1‑150 (MQ=255)

                                               cAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCACGCGTGAATACGc                                                                                                      <  1:172829/150‑1 (MQ=255)

                                                 aGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCg                                                                                                    <  2:1640617/150‑1 (MQ=255)

                                                  ggCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGa                                                                                                   >  1:499415/1‑150 (MQ=255)

                                                   gCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaa                                                                                                  >  2:1375603/1‑150 (MQ=255)

                                                    cTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaaa                                                                                                 <  2:138865/150‑1 (MQ=255)

                                                    cTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaaa                                                                                                 <  2:140517/150‑1 (MQ=255)

                                                     tCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACCACGTACTGCGTTAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATGCGCCGAAac                                                                                                >  2:919694/1‑150 (MQ=255)

                                                     tCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAac                                                                                                <  2:1467985/150‑1 (MQ=255)

                                                        ccTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATAc                                                                                                        <  1:913423/139‑1 (MQ=255)

                                                        ccTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATAc                                                                                                        >  2:913423/1‑139 (MQ=255)

                                                         cTACAGTCC‑GACGCGCACAAGATCCACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAaccacc                                                                                            >  2:1622817/1‑150 (MQ=255)

                                                          tACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGTAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAAccacca                                                                                           >  2:1484252/1‑150 (MQ=255)

                                                          tACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCTTGAATACGCCGAAAccacca                                                                                           >  1:1062626/1‑150 (MQ=255)

                                                          tACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAAccacca                                                                                           >  2:1875767/1‑150 (MQ=255)

                                                            cAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGaaa                                                                                                 >  2:944307/1‑142 (MQ=255)

                                                             aGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGc                                                                                        <  1:1622817/150‑1 (MQ=255)

                                                                 c‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCTCGTGAATACGCCGAAACCACCAAGCgtgt                                                                                    >  1:2192340/1‑150 (MQ=255)

                                                                   gACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCGTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGtgtg                                                                                   <  2:1655536/150‑1 (MQ=255)

                                                                   gACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGtgtg                                                                                   >  2:1488823/1‑150 (MQ=255)

                                                                      gcgcACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTgtgcgt                                                                                >  1:1371799/1‑150 (MQ=255)

                                                                        gcACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTgtgcgtgc                                                                              <  1:144143/150‑1 (MQ=255)

                                                                          acaAGATCAACCAGGCACTGGACATTCGGAACAACGCCCTGCGTGAGTTGACCGAGGACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                            >  2:1255870/1‑150 (MQ=255)

                                                                          acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                            >  2:491432/1‑150 (MQ=255)

                                                                          acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                            >  2:1978671/1‑150 (MQ=255)

                                                                          acaAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCgg                                                                            >  2:1931310/1‑150 (MQ=255)

                                                                              gATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTa                                                                        >  1:1030462/1‑150 (MQ=255)

                                                                              gATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTa                                                                        >  1:1622746/1‑150 (MQ=255)

                                                                                tCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGCCGAATCTTACTGTATGCAGGTG‑‑GACCACGCGTTAAGACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTAcc                                                                      >  1:280203/1‑150 (MQ=255)

                                                                                tCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaa                                                                                                                                 >  2:1684237/1‑91 (MQ=255)

                                                                                tCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaa                                                                                                                                 <  1:1684237/91‑1 (MQ=255)

                                                                                  aaCCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCCTGTGCGTGCGGCATACCAg                                                                    >  2:228241/1‑150 (MQ=255)

                                                                                    ccAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAggc                                                                  >  2:1010572/1‑149 (MQ=255)

                                                                                     cAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTAGCACCACCTGTTTCC‑T‑GCTAGCAAGCAGCTGCTTTCCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAc                                                                 <  1:1484252/150‑1 (MQ=255)

                                                                                          aCTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGtt                                                            >  2:630527/1‑150 (MQ=255)

                                                                                           cTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTc                                                           >  1:996431/1‑150 (MQ=255)

                                                                                            tGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCg                                                          >  2:79721/1‑150 (MQ=255)

                                                                                            tGGACATTCGCAACAACGCCCGGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCg                                                          <  2:227995/150‑1 (MQ=255)

                                                                                              gACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGc                                                        <  2:499415/150‑1 (MQ=255)

                                                                                                 aTTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTg                                                     <  2:51288/150‑1 (MQ=255)

                                                                                                 aTTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTg                                                     >  2:2306714/1‑150 (MQ=255)

                                                                                                  ttCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                    <  1:79721/150‑1 (MQ=255)

                                                                                                  ttCGCAACAACACCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                    <  1:1073434/150‑1 (MQ=255)

                                                                                                   tCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaagca                                                                                                                              >  1:336468/1‑75 (MQ=255)

                                                                                                   tCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCagcaagca                                                                                                                              <  2:336468/75‑1 (MQ=255)

                                                                                                   tCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTggg                                                   >  2:913793/1‑150 (MQ=255)

                                                                                                      caacaaCGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCa                                                                     <  2:652528/129‑1 (MQ=255)

                                                                                                      caacaaCGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCa                                                                     >  1:652528/1‑129 (MQ=255)

                                                                                                      caacaaCGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCa                                                <  1:1875767/150‑1 (MQ=255)

                                                                                                       aacaacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAg                                               >  1:80729/1‑150 (MQ=255)

                                                                                                        acaacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCGGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAgc                                              <  1:1978671/150‑1 (MQ=255)

                                                                                                          aacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAc                                                                 >  1:847181/1‑129 (MQ=255)

                                                                                                          aacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAc                                                                 <  2:847181/129‑1 (MQ=255)

                                                                                                          aacGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAgcgc                                            >  1:948563/1‑150 (MQ=255)

                                                                                                            cGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCAc                                          >  2:1900660/1‑150 (MQ=255)

                                                                                                             gCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                    <  2:1570940/139‑1 (MQ=255)

                                                                                                             gCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCAcc                                         <  1:630527/150‑1 (MQ=255)

                                                                                                                 tGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGACCAC‑A‑GTCCGTGGTCAGCCGCTCACCCCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCctgc                                     <  1:2113566/150‑1 (MQ=255)

                                                                                                                     tGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGa                                 <  1:1900660/150‑1 (MQ=255)

                                                                                                                           gACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCt                           >  1:2348062/1‑150 (MQ=255)

                                                                                                                           gACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCt                           <  1:2306714/150‑1 (MQ=255)

                                                                                                                            ttgctggcgttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGa                                                                  >  2:2200573‑M1/17‑110 (MQ=255)

                                                                                                                            ttgctggcgttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGa                                                                  <  1:2200573‑M1/94‑1 (MQ=255)

                                                                                                                            aCCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTgg                                                    >  1:964302/1‑124 (MQ=255)

                                                                                                                              cGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGAcc                                                                >  2:2199416/1‑110 (MQ=255)

                                                                                                                                agaACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACa                      <  2:80729/150‑1 (MQ=255)

                                                                                                                                agaACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACa                      >  1:2224784/1‑150 (MQ=255)

                                                                                                                                agaACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACa                      <  2:2224784/150‑1 (MQ=255)

                                                                                                                                    gttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGt                                                                                                              <  1:1522919‑M1/50‑1 (MQ=255)

                                                                                                                                    gttgttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGt                                                                                                              >  2:1522919‑M1/9‑58 (MQ=255)

                                                                                                                                    ccTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTga                  >  2:2379812/1‑150 (MQ=255)

                                                                                                                                     cTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTgag                 >  2:1632117/1‑150 (MQ=255)

                                                                                                                                     cTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTgag                 >  2:1175378/1‑150 (MQ=255)

                                                                                                                                        ttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGt                    <  2:195413‑M1/140‑1 (MQ=255)

                                                                                                                                        ttcaTGGCGTTCCACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGt                    >  1:195413‑M1/5‑144 (MQ=255)

                                                                                                                                          ccTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCTTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGa            <  2:1062626/150‑1 (MQ=255)

                                                                                                                                           cTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAg           <  2:948563/150‑1 (MQ=255)

                                                                                                                                            TGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGc          <  2:639085‑M1/150‑1 (MQ=255)

                                                                                                                                             GGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCt         >  1:1255870‑M1/1‑150 (MQ=255)

                                                                                                                                                 TGGACCGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTAc     >  1:111261‑M1/1‑150 (MQ=255)

                                                                                                                                                    aaCGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTACAAc  >  2:2379500/1‑150 (MQ=255)

CAAGCTGGCGTGGAGATGAACTACAACCTGTTCCGCGGTGGCAGCGACAAGGCTCGCCTACAGTCC‑GACGCGCACAAGATCAACCAGGCACTGGACATTCGCAACAACGCCCTGCGTGAGTTGACCGAGAACCTGAGCCTGGCGTGGAACGCCATGAACAA‑C‑GCCAGCAAGCAGCTGCCGACCGCGCGTGAATACGCCGAAACCACCAAGCGTGTGCGTGCGGCCTACCAGGACCAGTTCGGCCTGGGCCAGCGCACCCTGCTGGACGTCCTCGACAGTGAGAACGAGCTGTACAAC  >  NC_002947/5133864‑5134160

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