BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,378,768	G→T	10.6%	P256H (CCT→CAT) ‡	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,378,768	0	G	T	10.6%	553.0 / 44.1	179	P256H (CCT→CAT) ‡	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand): ref base G (73/87); new base T (7/12); total (80/99)
Fisher's exact test for biased strand distribution p-value = 6.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 179 total aligned reads displayed.

GCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCG  >  NC_002947/5378622‑5378911

 cTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATc                                                                                                                                             <  1:2384576/150‑1 (MQ=255)

 cTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATc                                                                                                                                             <  1:1257521/150‑1 (MQ=255)

 cTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATc                                                                                                                                             <  1:1565491/150‑1 (MQ=255)

  ttGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCg                                                                                                                                            >  2:3043039/1‑150 (MQ=255)

  ttGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCg                                                                                                                                            <  1:2870164/150‑1 (MQ=255)

    ggCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGtt                                                                                                                                          >  2:1364827/1‑150 (MQ=255)

       gcgGGAGATCTTCACGTTCAGGAGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttctt                                                                                                                                       <  2:1192160/150‑1 (MQ=255)

        cgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGTCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttcttc                                                                                                                                      <  2:828992/150‑1 (MQ=255)

        cgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttcttc                                                                                                                                      <  1:1168311/150‑1 (MQ=255)

        cgGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGttcttc                                                                                                                                      <  2:2206026/150‑1 (MQ=255)

          ggAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAg                                                                                                                                    <  2:1549872/150‑1 (MQ=255)

              aTCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                                                                <  2:3544080/150‑1 (MQ=255)

                 ttCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATg                                                                                                                             >  1:4906413/1‑150 (MQ=255)

                  tCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGc                                                                                                                            >  2:974030/1‑150 (MQ=255)

                   gattgcaaGGTGCTTAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTGCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGcc                                                                                                                           >  1:592795/8‑150 (MQ=255)

                   cACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGcc                                                                                                                           >  1:2415241/1‑150 (MQ=255)

                     ttgcaaGGTGCTTAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGccct                                                                                                                         <  1:2301631/145‑2 (MQ=255)

                          aGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGActct                                                                                                                    >  2:5302135/1‑150 (MQ=255)

                            gTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtcttt                                                                                                                  >  2:2835489/1‑150 (MQ=255)

                             tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                                                 <  1:4972818/150‑1 (MQ=255)

                             tGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttc                                                                                                                 <  1:355636/150‑1 (MQ=255)

                              gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                <  1:4838819/150‑1 (MQ=255)

                              gCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttct                                                                                                                <  2:603001/150‑1 (MQ=255)

                               cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCTGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                               <  1:369373/150‑1 (MQ=255)

                               cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                               <  1:1656585/150‑1 (MQ=255)

                               cTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACtctttctt                                                                                                               <  1:1364827/150‑1 (MQ=255)

                                ttAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                              >  1:2892962/1‑150 (MQ=255)

                                ttAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTg                                                                                                              >  2:3278734/1‑150 (MQ=255)

                                ttAGGACCGGTCTGGTCAGCCGTGATGAATGGCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTg                                                                                                              <  1:4533742/150‑1 (MQ=255)

                                 tAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGa                                                                                                             <  2:4906413/150‑1 (MQ=255)

                                 tAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGa                                                                                                             <  1:1591705/150‑1 (MQ=255)

                                    gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                                                                >  2:4300294/1‑128 (MQ=255)

                                    gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTcc                                                                                                                                <  1:4300294/128‑1 (MQ=255)

                                    gACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACt                                                                                                          >  1:2768720/1‑150 (MQ=255)

                                     aCCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtc                                                                                                         <  2:624619/150‑1 (MQ=255)

                                        ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                                                                      <  1:1772482/150‑1 (MQ=255)

                                        ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtcgtc                                                                                                      <  1:1772337/150‑1 (MQ=255)

                                        ggTGGCATCTGCAGTGATGTACGGCAGGTTGACGACGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGCACtcgtc                                                                                                      <  2:5405738/150‑1 (MQ=255)

                                          tGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGa                                                                                                    <  1:974030/150‑1 (MQ=255)

                                              aTCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAg                                                                                                <  2:3728169/150‑1 (MQ=255)

                                              aTCTGCAGTGATGTAAGGCAGGTTGACTTCGGTCCACTGAGCCGACGACAGCTCGATCATGGCTTTATCTGCGGCTACTTTCAGACGCTGCACGGCCAGAGTATCGTTCTTCAGGACCATGCACGACTCTTTCTTGAACTCGTCGACAAg                                                                                                >  2:3852006/1‑150 (MQ=255)

                                               tCTGCAGTGATGTACCGCAGGTTGACGTCGGTCGACTGTGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAgg                                                                                               >  2:3347300/1‑150 (MQ=255)

                                                   cAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAACTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtc                                                                                                         <  1:1750004/136‑1 (MQ=255)

                                                   cAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAACTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACtc                                                                                                         >  2:1750004/1‑136 (MQ=255)

                                                       gATGTACGGCAGGTTGACGTCGGTCGACTGAGAGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGAt                                                                                       >  1:4277558/1‑150 (MQ=255)

                                                         tgaatggCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCa                                                                                     <  1:3833449/145‑1 (MQ=255)

                                                         tgaatggCAGGTTGATTTCGGTCTGCTGAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCa                                                                                     <  2:2413450/145‑1 (MQ=255)

                                                          gTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAg                                                                                    >  2:4809859/1‑150 (MQ=255)

                                                                 aGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATg                                                                             >  1:1266681/1‑150 (MQ=255)

                                                                  ggTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGt                                                                            >  1:665210/1‑150 (MQ=255)

                                                                  ggTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGt                                                                            <  2:4085696/150‑1 (MQ=255)

                                                                   gTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTc                                                                           >  1:597710/1‑150 (MQ=255)

                                                                    ttGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCg                                                                          <  2:5128637/150‑1 (MQ=255)

                                                                     tGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGa                                                                         <  1:3278734/150‑1 (MQ=255)

                                                                      gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                        <  2:2955827/150‑1 (MQ=255)

                                                                      gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                        >  2:2634006/1‑150 (MQ=255)

                                                                      gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                        <  2:3136956/150‑1 (MQ=255)

                                                                      gACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGaa                                                                        <  2:1301451/150‑1 (MQ=255)

                                                                            ggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTc                                                                  <  2:3913436/150‑1 (MQ=255)

                                                                            ggtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTc                                                                  <  2:4658246/150‑1 (MQ=255)

                                                                             gtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCg                                                                 >  2:1773412/1‑150 (MQ=255)

                                                                             gtcgACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCg                                                                 <  1:2186037/150‑1 (MQ=255)

                                                                                    gAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccga                                                          <  2:1089387/150‑3 (MQ=255)

                                                                                    gAGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGccga                                                          <  1:386682/150‑3 (MQ=255)

                                                                                     aGCGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAg                                                         <  1:863700/150‑1 (MQ=255)

                                                                                       cGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGa                                                       >  1:4501987/1‑150 (MQ=255)

                                                                                       cGGACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGacg                                                                                 <  2:1731704/124‑1 (MQ=255)

                                                                                        ggACGACAGCTCGATCTTGGCTTTTTCTGCGGATTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGaa                                                      <  2:4477336/150‑1 (MQ=255)

                                                                                        ggACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTTCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGaa                                                      >  1:3282139/1‑150 (MQ=255)

                                                                                        ggACGAAAGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGaa                                                      <  2:3484172/150‑1 (MQ=255)

                                                                                         gacgacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAg                                                     <  2:4231532/150‑1 (MQ=255)

                                                                                          acgacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAgg                                                                                               >  1:3759775/1‑107 (MQ=255)

                                                                                            gacAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTc                                                                           >  1:1664400/1‑125 (MQ=255)

                                                                                            gaaaGCTCGATCTTGGCTTTTTCGGCAGCTTCCTTGAGGCGCTGCAGAGCAAGCTTGTCTTTCTTAAGATCGATGCCCTGTTCCTTCTTGAACTCGGCAACCAGATATTCGACCAGACGCATGTCAAAGTCTTCACCGCCGAGGAACgtg                                                  >  1:4805648/4‑150 (MQ=255)

                                                                                              cAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc                                                <  1:1051144/150‑1 (MQ=255)

                                                                                              cAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc                                                <  1:2317474/150‑1 (MQ=255)

                                                                                                  tCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCGGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCg                                            >  2:4035620/1‑150 (MQ=255)

                                                                                                   cGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGt                                           >  1:806735/1‑150 (MQ=255)

                                                                                                    gATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGtt                                          <  1:3056297/150‑1 (MQ=255)

                                                                                                        ttGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTa                                      >  1:1693997/1‑150 (MQ=255)

                                                                                                           gCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                                                                             >  1:642117/1‑92 (MQ=255)

                                                                                                           gCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTa                                                                                             <  2:642117/92‑1 (MQ=255)

                                                                                                            cTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCa                                  >  2:3833893/1‑150 (MQ=255)

                                                                                                               tttCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTa                               <  1:5134262/150‑1 (MQ=255)

                                                                                                               tttCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTa                               <  2:2498477/150‑1 (MQ=255)

                                                                                                                 tCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACt                             >  2:504798/1‑150 (MQ=255)

                                                                                                                 tCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACt                             >  2:4348891/1‑150 (MQ=255)

                                                                                                                   tGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTc                           <  2:1036615/150‑1 (MQ=255)

                                                                                                                   tGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTc                           >  1:1000104/1‑150 (MQ=255)

                                                                                                                   tGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTc                           >  1:1339230/1‑150 (MQ=255)

                                                                                                                          tcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTgg                    <  1:2640523/150‑1 (MQ=255)

                                                                                                                          tcttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTgg                    <  1:5273847/150‑1 (MQ=255)

                                                                                                                           cttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGt                   >  2:1430098/1‑150 (MQ=255)

                                                                                                                           cttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGt                   <  1:3208567/150‑1 (MQ=255)

                                                                                                                           cttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGt                   <  1:3347300/150‑1 (MQ=255)

                                                                                                                           cttTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGt                   <  2:3194387/150‑1 (MQ=255)

                                                                                                                                 gACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTc                                                >  1:1004507/1‑115 (MQ=255)

                                                                                                                                 gACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTggcta                 >  1:3132335/1‑143 (MQ=255)

                                                                                                                                            gCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCg  >  2:687009/1‑150 (MQ=255)

                                                                                                                                              cAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGc                                    <  1:4057684/114‑1 (MQ=255)

GCTTGGCGCGGGAGATCTTCACGTTCAGGTGCTTAGGACCGGTGGCATCTGCAGTGATGTACGGCAGGTTGACGTCGGTCGACTGAGCGGACGACAGCTCGATCTTGGCTTTTTCTGCGGCTTCTTTCAGACGCTGCAGGGCCAGAGGATCGTTCTTCAGGTCCATGCCCGACTCTTTCTTGAACTCGTCGACAAGGTAGTCGATCAGGCGCATGTCGAAGTCTTCGCCACCCAGGAAGGTGTCGCCGTTGGTAGCCAGTACTTCGAACTGGTGCTCACCGTCGACTTCG  >  NC_002947/5378622‑5378911

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