BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,137,858	G→C	5.5%	L228V (CTC→GTC)	PP_4521 ←	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,137,858	0	G	C	5.5%	545.9 / 16.1	163	L228V (CTC→GTC)	PP_4521	aerotaxis receptor
Reads supporting (aligned to +/- strand): ref base G (59/95); new base C (5/4); total (64/99)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 164 total aligned reads displayed.

TCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATCAGCA  >  NC_002947/5137710‑5138007

 cgctgcgctgCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCt                                                                                                                                                     <  2:3619697/150‑1 (MQ=255)

 cgctgcgctgCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCt                                                                                                                                                     <  2:1085623/150‑1 (MQ=255)

  gctgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTg                                                                                                                                                    <  1:3255477/150‑1 (MQ=255)

  gctgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTg                                                                                                                                                    >  1:2200006/1‑150 (MQ=255)

   ctgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTgg                                                                                                                                                   <  2:606543/150‑1 (MQ=255)

    tgcgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGc                                                                                                                                                  >  1:2270382/1‑150 (MQ=255)

      cgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                                <  1:356574/150‑2 (MQ=255)

      cgctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGcca                                                                                                                                                <  1:1744327/150‑2 (MQ=255)

           ccGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAg                                                                                                                                           >  2:1799080/1‑150 (MQ=255)

             ggTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGc                                                                                                                                         >  2:372355/1‑150 (MQ=255)

             ggTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGc                                                                                                                                         >  2:1863288/1‑150 (MQ=255)

                cGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTg                                                                                                                                      <  2:5035675/150‑1 (MQ=255)

                cGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTg                                                                                                                                      <  1:1932817/150‑1 (MQ=255)

                cGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTg                                                                                                                                      <  1:1932663/150‑1 (MQ=255)

                 gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                                     <  1:1854092/150‑1 (MQ=255)

                 gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                                     <  1:501521/150‑1 (MQ=255)

                 gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTCGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                                     >  2:654125/1‑150 (MQ=255)

                    cATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGtt                                                                                                                                  >  1:3191810/1‑150 (MQ=255)

                       cggcgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGCCCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCGGATGTTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                               <  1:4819547/148‑1 (MQ=255)

                       cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                               <  1:2365018/150‑1 (MQ=255)

                       cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                               <  1:1387352/150‑1 (MQ=255)

                       cTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTc                                                                                                                               <  1:51063/150‑1 (MQ=255)

                          cgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCt                                                                                                                            <  2:745558/150‑1 (MQ=255)

                          cgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCt                                                                                                                            <  1:1355948/150‑1 (MQ=255)

                            ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGt                                                                                                                          <  1:4744321/150‑1 (MQ=255)

                            ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGt                                                                                                                          <  2:4116000/150‑1 (MQ=255)

                            ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGt                                                                                                                          <  2:1221645/150‑1 (MQ=255)

                            ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGt                                                                                                                          <  1:4014888/150‑1 (MQ=255)

                            ccTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGt                                                                                                                          <  2:4990282/150‑1 (MQ=255)

                             cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                         <  1:2926994/150‑1 (MQ=255)

                             cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                         <  2:2984080/150‑1 (MQ=255)

                             cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                         <  1:758386/150‑1 (MQ=255)

                             cTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTa                                                                                                                         >  2:4657629/1‑150 (MQ=255)

                              tGGGCATCGGACAAGTCGTAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCTTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAg                                                                                                                        <  1:2278754/150‑1 (MQ=255)

                              tGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAg                                                                                                                        <  1:1197594/150‑1 (MQ=255)

                              tGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAg                                                                                                                        <  1:3403468/150‑1 (MQ=255)

                              tGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAg                                                                                                                        <  2:1541668/150‑1 (MQ=255)

                                 gCATCGGACAAGTCGGAAGAGTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTc                                                                                                                     <  1:2804350/150‑1 (MQ=255)

                                   aTCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAg                                                                                                                   >  1:2654050/1‑150 (MQ=255)

                                    tCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAgc                                                                                                                  <  2:4353788/150‑1 (MQ=255)

                                    tCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGGTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAgc                                                                                                                  >  2:507720/1‑150 (MQ=255)

                                       gACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGcc                                                                                                               >  2:2407401/1‑150 (MQ=255)

                                       gACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCGTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCTCTGTAGGTCAGCGcc                                                                                                               >  1:3948700/1‑150 (MQ=255)

                                          aaGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccacc                                                                                                            >  1:1568309/1‑150 (MQ=255)

                                             tCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAAc                                                                                                         <  1:1310626/150‑1 (MQ=255)

                                              cGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACg                                                                                                        <  2:933886/150‑1 (MQ=255)

                                               ggAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACgg                                                                                                       >  1:4505300/1‑150 (MQ=255)

                                                 aaGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGAc                                                                                                     <  1:906283/150‑1 (MQ=255)

                                                        cgccgcCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTaaa                                                                                              <  2:222842/150‑1 (MQ=255)

                                                         gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGTGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAg                                                                                             >  1:356104/1‑150 (MQ=255)

                                                         gccgccCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGaaaca                                                                                             >  1:2169766/1‑148 (MQ=255)

                                                             cccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACGGGTCAAGGCCt                                                                                         >  2:3073924/1‑150 (MQ=255)

                                                               cgcgACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTg                                                                                       <  1:3003834/150‑1 (MQ=255)

                                                                gcgACTTTACCCCGACCCCGGCATCCACCAGGGGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTgg                                                                                      <  2:4492945/150‑1 (MQ=255)

                                                                  gACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGa                                                                                    <  2:4064712/150‑1 (MQ=255)

                                                                   aCTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGAt                                                                                   <  1:2369367/150‑1 (MQ=255)

                                                                     tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGc                                                                                 >  1:2907131/1‑150 (MQ=255)

                                                                     tttACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGc                                                                                 >  2:1946839/1‑150 (MQ=255)

                                                                        cccccgaccccgGCATCCACCAGGCGGGTCTGGGCCGTGTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTGGGGATGCTCa                                                                              <  2:2289568/149‑1 (MQ=255)

                                                                        cccccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCa                                                                              <  2:3790066/149‑1 (MQ=255)

                                                                        accccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATCGCCAGGTCGAGCTGGCCCATCAGGCCTGGGTTGTCGCTGTAGGTCAGAGCCACCAACGGACTGGTAAAGGCCATGGGATGCTCa                                                                              <  2:4304906/150‑1 (MQ=255)

                                                                         ccccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAg                                                                             >  2:1889035/1‑150 (MQ=255)

                                                                          cccgaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGc                                                                            >  2:145628/1‑150 (MQ=255)

                                                                             gaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCaa                                                                         <  1:1775908/150‑1 (MQ=255)

                                                                             gaccccgGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCaa                                                                         <  1:1586000/150‑1 (MQ=255)

                                                                                     cATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcgac                                                                 <  1:4657629/150‑1 (MQ=255)

                                                                                       tCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGcgacga                                                               >  2:1060443/1‑150 (MQ=255)

                                                                                       tCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCAGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTAGGGATGCTCAGCCAACGTGcgacga                                                               >  1:2730861/1‑150 (MQ=255)

                                                                                          accaGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATg                                                            >  2:3262412/1‑150 (MQ=255)

                                                                                            caGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACAGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGc                                                          >  1:2115924/1‑150 (MQ=255)

                                                                                               gCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTg                                                       <  2:2200006/150‑1 (MQ=255)

                                                                                               gCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTg                                                       >  1:3379406/1‑150 (MQ=255)

                                                                                                    tCAGGGCCGTTTGCAGCCGGGCTTCTTCGCTGATGCTGGCCCGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCt                                                  <  2:2279047/150‑1 (MQ=255)

                                                                                                    tCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCt                                                  >  2:3859752/1‑150 (MQ=255)

                                                                                                        ggCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCg                                              >  2:3587800/1‑150 (MQ=255)

                                                                                                           ggTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTg                                           <  1:2360959/149‑1 (MQ=255)

                                                                                                            gTTTGCTGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                          <  1:5019001/150‑1 (MQ=255)

                                                                                                            gTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                          <  2:3539248/150‑1 (MQ=255)

                                                                                                            gTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGc                                          <  2:3539581/150‑1 (MQ=255)

                                                                                                                gCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGa                                      <  2:1568309/150‑1 (MQ=255)

                                                                                                                  aGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGAcc                                    >  2:2363792/1‑150 (MQ=255)

                                                                                                                  aGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGAcc                                    >  2:4664700/1‑150 (MQ=255)

                                                                                                                      gggCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGt                                >  1:2718479/1‑150 (MQ=255)

                                                                                                                             ttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc                         >  1:732277/1‑150 (MQ=255)

                                                                                                                             ttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc                         >  1:3369824/1‑150 (MQ=255)

                                                                                                                             ttcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGc                         <  1:4178670/150‑1 (MQ=255)

                                                                                                                              tcGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGcc                        <  2:2052151/150‑1 (MQ=255)

                                                                                                                               cGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCca                       <  2:3449989/150‑1 (MQ=255)

                                                                                                                                  tgatgaTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagca                    <  2:4748168/150‑1 (MQ=255)

                                                                                                                                  tgatgaTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagca                    <  2:1822164/150‑1 (MQ=255)

                                                                                                                                  cctgaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagca                    >  1:1667421/8‑150 (MQ=255)

                                                                                                                                     tgaggGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcagca                 <  1:4295486/150‑1 (MQ=255)

                                                                                                                                     gaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcagca                 >  2:850673/5‑150 (MQ=255)

                                                                                                                                     gaagggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCcagcagca                 >  2:850703/5‑150 (MQ=255)

                                                                                                                                        gggCCAGCTCGACCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGAc              >  2:1552508/2‑150 (MQ=255)

                                                                                                                                          gCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACtt            <  1:1723865/150‑1 (MQ=255)

                                                                                                                                           ccAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTg           >  2:2692404/1‑150 (MQ=255)

                                                                                                                                           ccAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTg           >  2:2692228/1‑150 (MQ=255)

                                                                                                                                                tCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCTGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATc      <  1:372355/150‑1 (MQ=255)

                                                                                                                                                tCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATc      <  2:3191810/150‑1 (MQ=255)

TCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCATCAGCA  >  NC_002947/5137710‑5138007

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