BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,632,139	A→G	8.6%	E284G (GAG→GGG)	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,632,139	0	A	G	8.6%	482.6 / 24.4	164	E284G (GAG→GGG)	lepA	elongation factor 4
Reads supporting (aligned to +/- strand): ref base A (72/76); new base G (9/5); total (82/82)
Fisher's exact test for biased strand distribution p-value = 4.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 166 total aligned reads displayed.

ACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGC  >  NC_002947/1631992‑1632287

aCTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAgg                                                                                                                                                    >  1:807619/1‑150 (MQ=255)

  ttGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTc                                                                                                                                                  >  2:3858018/1‑150 (MQ=255)

    gggcGACAGCGTCGGGGTCGTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGa                                                                                                                                                <  2:588266/147‑1 (MQ=255)

    ggTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGa                                                                                                                                                <  2:588239/150‑1 (MQ=255)

     gTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAg                                                                                                                                               <  2:1229397/150‑1 (MQ=255)

     gTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAg                                                                                                                                               <  1:2084157/150‑1 (MQ=255)

       cGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgt                                                                                                                                             <  2:2964058/150‑1 (MQ=255)

        gACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgtg                                                                                                                                            >  2:4828195/1‑150 (MQ=255)

         aCAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAggaca                                                                                                                                           <  2:1576263/150‑4 (MQ=255)

           aGCGTCGTTGTCTTCACCCCGAAGCTCTCCCAGTCCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATGATCGCCAGCATCAAGGACATCCACGGTGCCCCGTTTGGTGTGACCCTGACCCTGTCCTCGACCTCAGGTGTCGAggacagg                                                                                                                                         <  2:4193573/150‑6 (MQ=255)

              gTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                      <  2:2803562/150‑1 (MQ=255)

              gTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                      <  1:688109/150‑1 (MQ=255)

               tCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAg                                                                                                                                     <  2:3614994/150‑1 (MQ=255)

                cGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAgg                                                                                                                                    <  1:3410615/150‑1 (MQ=255)

                  gtgtCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCt                                                                                                                                  >  2:2910854/1‑150 (MQ=255)

                  gtgtCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGAGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCt                                                                                                                                  >  1:2814076/1‑150 (MQ=255)

                  gtgtCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGAGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCt                                                                                                                                  >  1:2814059/1‑150 (MQ=255)

                     tCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCa                                                                                                                               >  2:1369532/1‑150 (MQ=255)

                      cTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                              <  1:2268263/150‑1 (MQ=255)

                       ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                             <  1:4498698/150‑1 (MQ=255)

                        tCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaaga                                                                                                                            >  2:3169433/1‑150 (MQ=255)

                            cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAt                                                                                                                        >  1:1659760/1‑150 (MQ=255)

                             ccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAggacagg                                                                                                                                         <  1:549596/132‑6 (MQ=255)

                             ccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAggacagg                                                                                                                                         >  2:549596/1‑127 (MQ=255)

                             ccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGAtt                                                                                                                       <  2:410399/150‑1 (MQ=255)

                              cGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGcc                                                                                                                                       <  1:1275618/133‑1 (MQ=255)

                               gAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCa                                                                                                                     <  1:1754956/150‑1 (MQ=255)

                                aaGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAg                                                                                                                    <  2:211850/150‑1 (MQ=255)

                                 aGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGc                                                                                                                   <  1:880347/150‑1 (MQ=255)

                                     cacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCCTCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGcc                                                                                                               <  2:4598191/150‑2 (MQ=255)

                                     cacCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCa                                                                                                               <  1:163668/150‑1 (MQ=255)

                                      acCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAg                                                                                                              >  2:1880033/1‑150 (MQ=255)

                                        ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                            >  1:108064/1‑150 (MQ=255)

                                        ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                            >  1:3577181/1‑150 (MQ=255)

                                        ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGt                                                                                                            >  2:786924/1‑150 (MQ=255)

                                          aGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCt                                                                                                          >  1:3934609/1‑150 (MQ=255)

                                               gCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGcc                                                                                                     >  2:514951/1‑150 (MQ=255)

                                                  gATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGc                                                                                                  <  1:3132464/150‑1 (MQ=255)

                                                  gATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGc                                                                                                  <  1:2414295/150‑1 (MQ=255)

                                                  gATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGc                                                                                                  <  1:4828195/150‑1 (MQ=255)

                                                      tGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGt                                                                                              >  2:5069564/1‑150 (MQ=255)

                                                       gAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGtt                                                                                             >  1:3341548/1‑150 (MQ=255)

                                                        aaaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTc                                                                                            <  1:2152740/150‑1 (MQ=255)

                                                           gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCg                                                                                         <  1:472439/150‑1 (MQ=255)

                                                           gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCg                                                                                         <  1:3844059/150‑1 (MQ=255)

                                                             cGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGt                                                                                       >  2:1570623/1‑150 (MQ=255)

                                                             cGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGt                                                                                       <  1:442597/150‑1 (MQ=255)

                                                               gCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCa                                                                                     <  2:744911/150‑1 (MQ=255)

                                                                   aGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTc                                                                                 <  1:230123/150‑1 (MQ=255)

                                                                    gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                                <  1:124315/150‑1 (MQ=255)

                                                                    gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                                <  2:1377535/150‑1 (MQ=255)

                                                                              tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCg                                                                      >  1:2572996/1‑150 (MQ=255)

                                                                               catcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                     <  2:5215702/150‑1 (MQ=255)

                                                                               catcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                     <  2:4088238/150‑1 (MQ=255)

                                                                                atcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAg                                                                    <  2:1462694/150‑1 (MQ=255)

                                                                                      aGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTc                                                              >  1:2174291/1‑150 (MQ=255)

                                                                                       gCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTcc                                                             >  1:5386563/1‑150 (MQ=255)

                                                                                        cATCAAGGACATCCTCGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCcg                                                            <  2:1006651/150‑1 (MQ=255)

                                                                                             aGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACg                                                       >  1:3262425/1‑150 (MQ=255)

                                                                                                   tCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                                 >  2:3544772/1‑150 (MQ=255)

                                                                                                    ccACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCa                                                <  2:2814076/150‑1 (MQ=255)

                                                                                                    ccACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCa                                                <  1:5059197/150‑1 (MQ=255)

                                                                                                      aCGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGa                                              >  2:4207837/1‑150 (MQ=255)

                                                                                                        ggTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAg                                            <  2:3552699/150‑1 (MQ=255)

                                                                                                        ggTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAg                                            >  2:2940762/1‑150 (MQ=255)

                                                                                                         gTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGc                                           <  2:4153634/150‑1 (MQ=255)

                                                                                                            ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                        <  2:2143748/150‑1 (MQ=255)

                                                                                                            ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                        <  2:2627406/150‑1 (MQ=255)

                                                                                                             cccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAc                                       >  2:5240327/1‑150 (MQ=255)

                                                                                                              ccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAcc                                      >  1:3793993/1‑150 (MQ=255)

                                                                                                                 gTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                                 <  1:4550816/136‑1 (MQ=255)

                                                                                                                    ggTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGTAGAAGCTGACCCTGAAc                                <  1:4591598/150‑1 (MQ=255)

                                                                                                                      tGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAg                                                                    <  1:1401624/112‑1 (MQ=255)

                                                                                                                      tGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGa                              >  2:4863604/1‑150 (MQ=255)

                                                                                                                          accctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcc                          >  2:5386075/1‑150 (MQ=255)

                                                                                                                            cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                        >  2:855354/1‑150 (MQ=255)

                                                                                                                            cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                        >  1:2340092/1‑150 (MQ=255)

                                                                                                                            cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                        >  1:3946199/1‑150 (MQ=255)

                                                                                                                             ctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCg                       <  2:1670372/150‑1 (MQ=255)

                                                                                                                             ctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCg                       <  2:4453803/150‑1 (MQ=255)

                                                                                                                               gaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                     <  2:4342679/102‑1 (MQ=255)

                                                                                                                               gaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                     >  1:4342679/1‑102 (MQ=255)

                                                                                                                                  cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa                  >  2:3983666/1‑150 (MQ=255)

                                                                                                                                  cctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCa                  >  1:2320800/1‑150 (MQ=255)

                                                                                                                                   aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                                <  1:1932353/83‑1 (MQ=38)

                                                                                                                                   aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTcc                                                                                >  2:1932353/5‑87 (MQ=38)

                                                                                                                                    tgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGt                >  1:5075820/1‑150 (MQ=255)

                                                                                                                                    tgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGt                >  1:3847908/1‑150 (MQ=255)

                                                                                                                                      accTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTAc              >  1:5318223/2‑150 (MQ=255)

                                                                                                                                      accTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTAc              >  2:507631/2‑150 (MQ=255)

                                                                                                                                        cTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACAt            >  2:3193627/1‑150 (MQ=255)

                                                                                                                                         tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGt                >  2:5035250/1‑145 (MQ=255)

                                                                                                                                          cGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGc          <  1:295851/150‑1 (MQ=255)

                                                                                                                                           gACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGcc         >  1:1084120/1‑150 (MQ=255)

                                                                                                                                            aCCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCg        <  1:3169433/150‑1 (MQ=255)

                                                                                                                                               ccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACGTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaa     >  2:2463540/1‑150 (MQ=255)

                                                                                                                                                cTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaaa    <  1:5069564/150‑1 (MQ=255)

                                                                                                                                                  gAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGc  >  1:242919/1‑150 (MQ=255)

                                                                                                                                                   aGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgacga                                                                           >  2:2040410/1‑76 (MQ=255)

                                                                                                                                                   aGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCcgacga                                                                           <  1:2040410/76‑1 (MQ=255)

ACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAGGTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAAGC  >  NC_002947/1631992‑1632287

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