BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,137,853	G→C	5.5%	D229E (GAC→GAG)	PP_4521 ←	aerotaxis receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,137,853	0	G	C	5.5%	544.3 / 16.1	164	D229E (GAC→GAG)	PP_4521	aerotaxis receptor
Reads supporting (aligned to +/- strand): ref base G (61/93); new base C (5/4); total (66/98)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 167 total aligned reads displayed.

TCAGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCAT  >  NC_002947/5137704‑5138002

tCAGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAgg                                                                                                                                                       >  1:1544041/1‑150 (MQ=255)

tCAGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAgg                                                                                                                                                       >  2:2857429/1‑150 (MQ=255)

  aGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTc                                                                                                                                                     <  1:4367538/150‑1 (MQ=255)

  aGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTc                                                                                                                                                     >  1:2289947/1‑150 (MQ=255)

    tctcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGa                                                                                                                                                   <  2:1913134/150‑1 (MQ=255)

      tcGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGc                                                                                                                                                 >  2:2369367/1‑150 (MQ=255)

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

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

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

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

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

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

             gctgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGCTCGACCTGGccat                                                                                                                                          <  2:518290/150‑3 (MQ=255)

               tgcCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTc                                                                                                                                        <  1:271245/150‑1 (MQ=255)

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

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

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

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

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

                       gAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTgg                                                                                                                                >  1:44293/1‑150 (MQ=255)

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

                        aGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTggg                                                                                                                               <  1:2646030/150‑1 (MQ=255)

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

                              tGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCg                                                                                                                         >  1:4935766/1‑150 (MQ=255)

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

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

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

                                cgcCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCAGCCACCCGGCGGGTCAGGGCCGTGTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCGGGCCCTTCAGGCCTGGGTTGTCGCt                                                                                                                       <  2:1035301/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                      ggCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGt                                                                                                                 <  1:1472063/150‑1 (MQ=255)

                                       gCATCGGACCAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCTGGGCCTTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTc                                                                                                                <  1:2804332/150‑1 (MQ=255)

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

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

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

                                              aCAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGcca                                                                                                         <  1:3327913/150‑1 (MQ=255)

                                                   tCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGccaccacc                                                                                                    >  1:2834631/1‑148 (MQ=255)

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

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

                                                         gATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTg                                                                                              <  2:2072406/150‑1 (MQ=255)

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

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

                                                                 cgcccgcGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGc                                                                                      >  1:1904202/1‑150 (MQ=255)

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

                                                                      gcgACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTgg                                                                                 >  1:2976781/1‑150 (MQ=255)

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

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

                                                                          cTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATg                                                                             >  1:1702513/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                  ctggCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGc                                                     <  1:4744444/148‑1 (MQ=255)

                                                                                                     gCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGCTGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTg                                                  <  1:3541576/150‑1 (MQ=255)

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

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

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

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

                                                                                                          tCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCt                                             <  2:75709/150‑1 (MQ=255)

                                                                                                            aGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGc                                           <  2:1217716/150‑1 (MQ=255)

                                                                                                                 cGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTg                                      <  1:2361092/150‑1 (MQ=255)

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

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

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

                                                                                                                      gCAGACGGGCTTCTTCGCTGATGATGGCGAGTTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGa                                 <  2:4930056/150‑1 (MQ=255)

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

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

                                                                                                                        aGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTAGGGATGCTCAGCCAACGTGAGACGAATGGCCTGGTTCTGCCGGTGCTCGAcc                               >  1:1728132/1‑150 (MQ=255)

                                                                                                                           cGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAg                            <  1:1799080/150‑1 (MQ=255)

                                                                                                                              gCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTAc                         >  2:24120/1‑150 (MQ=255)

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

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

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

                                                                                                                                     ggcggttgatgGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTGAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCca                  <  2:2270382/144‑1 (MQ=255)

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

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

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

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

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

                                                                                                                                                 ccAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTg      <  1:2306796/150‑1 (MQ=255)

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

                                                                                                                                                     gTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCAt  >  1:5069939/1‑150 (MQ=255)

TCAGTCTCGCTGCGCTGCCGGTCGAGCATCTGCGCCTGGGCATCGGACAAGTCGGAAGATTGCGCCGCCCGCGACTTTACCCCGACCCCGGCATCCACCAGGCGGGTCAGGGCCGTTTGCAGACGGGCTTCTTCGCTGATGATGGCCAGGTCGAGCTGGCCCTTCAGGCCTGGGTTGTCGCTGTAGGTCAGCGCCACCAACGGACTGGTAAAGGCCTTGGGATGCTCAGCCAACGTGCGACGAATGGCCTGGTTCTGCCGGTGCTCGACCCAGTACCAAGCCAGCAGCAGACTTGCCAT  >  NC_002947/5137704‑5138002

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